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In  this  thesis  we  present  the  first  measurements  of  simultane- 
ously recorded  lightning  return  stroke  vertical  electric  field  and 
induced  voltage  on  a  power  line.  Data  are  given  for  about  100  first 
strokes  and  over  200  subsequent  strokes  with  the  lightning  ground 
strike  point  located  by  a  combination  of  a  TV  network,  thunder 
ranging,  and  observer  comments.  Voltages  were  measured  at  one  end  of 
a  460  m  distribution  line  which  was  specially  constructed  for  the 
experiment.  The  voltage  and  vertical  electric  field  time-domain 
signals  are  recorded  on  film  from  oscilloscope  faces  and  on  magnetic 
instrumentation  tape  with  an  effective  frequency  response  from  less 
than  1  Hz  to  over  1  MHz.  The  measured  voltage  waveforms  do  not 
resemble  the  "classical  induced  positive  surge"  described  in  the 

literature.   Both  the  magnitude  and  polarity  of  the  induced  voltages 

xix 


are  found  to  be  strong  functions  of  the  location  of  the  lightning 
ground  strike  point. 

Theory  is  presented  which  predicts  induced  surges  of  both 
polarities,  as  measured.  First, the  return  stroke  horizontal  electric 
field  associated  with  a  finite  ground  conductivity  is  derived  from 
the  measured  vertical  electric  field.  Then,  the  Telegraphers' 
Equations  are  solved  by  computer  using  the  vertical  and  horizontal 
electric  fields  as  forcing  functions  to  predict  the  line  voltages. 
Calculated  voltage  waveforms  are  presented  for  a  500  m  line  and  for  a 
5  km  line  for  earth  conductivities  ranging  from  10"^  mhos/m  to  10"^ 
mhos/m,  earth  permittivities  ranging  from  s  =  15  to  e  =  3,  and 
lightning  ground  strike  points  ranging  from  0.2  km  to  5.0  km  from  the 
line. 

Direct  comparisons  are  made  between  the  measured  voltages  on  the 
460  m  line  and  calculated  voltage  waveforms.  Calculated  waveshapes 
are  in  moderately  good  agreement  with  the  measurements.  However, 
there  is  a  consistent  discrepancy  in  the  magnitudes:  the  calculated 
values  are  about  a  factor  of  4  lower  than  those  measured.  Possible 
errors  in  both  theory  and  measurement  are  discussed.  Where  possible, 
comparisons  are  made  between  voltage  measurements  reported  by  other 
investigators  and  our  theory. 

In  addition  to  the  lightning  return  stroke  data,  we  present  the 
first  stepped  leader  induced  voltages  correlated  with  stepped  leader 
vertical  electric  fields.  The  voltages  induced  by  stepped  leaders  on 
the  450  m  test  line  are  an  order  of  magnitude  lower  than  the  voltage 
induced  by  the  following  return  stroke. 


CHAPTER  I 
REVIEW 

1.1  Introduction 

In  this  chapter  a  brief  discussion  of  the  lightning  discharge 
process  will  be  presented  and  the  terminology  used  in  this  thesis 
will  be  defined.  Attention  will  be  focussed  on  the  cloud-to-ground 
lightning  discharge  because  it  is  this  discharge  which  is  of 
relevance  to  power  line  problems.  The  second  half  of  the  chapter 
will  be  devoted  to  a  critical  discussion  of  the  work  which  has  been 
done  to  date  regarding  the  understanding  of  lightning  induced 
voltages  on  power  lines. 

1.2  Lightning 

1 .2.1  The  Thundercloud        ' 

The  present  discussion  of  the  lightning  discharge  process  is  not 
intended  to  be  complete  or  exhaustive;  it  is  included  here  for  the 
sake  of  completeness.  The  discussion  presented  here  is  more  or  less 
the  consensus  view  of  the  discharge  process.  For  a  more  detailed 
discussion,  the  reader  may  refer  to  several  standard  monographs, 
e.g.,  Malan  (1963),  Chalmers  (1967),  Uman  (1969),  and  Golde 
(1977a).   Lightning  has  been  defined  by  Uman  (1969)  as  "a  transient, 
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high-current  electric  discharge,  whose  path  length  is  generally 
measured  in  kilometers."  The  most  common  source  of  lightning  is  the 
electrical  charge  separation  in  the  cumulonimbus  cloud.  All 
lightning  discharges  can  be  roughly  divided  into  cloud  discharges  and 
ground  discharges.  Lightning  discharges  which  take  place  entirely 
within  a  cloud,  between  two  clouds,  or  between  a  cloud  and  the 
surrounding  air  are  called  cloud  discharges  or  cloud  lightning. 
Lightning  discharges  which  take  place  between  the  cloud  and  the 
ground  are  referred  to  as  cloud-to-ground  discharges  or  ground 
lightning.  The  ratio  of  cloud  lightning  to  ground  lightning  has  been 
reported  by  Pierce  (1970)  and  Prentice  and  Mackerras  (1977)  to  be  a 
function  of  the  geographical  latitude.  Livingston  and  Krider  (1978) 
have  found  that  for  individual  storms,  the  ratio  is  a  function  of 
storm  stage. 

The  mechanisms  for  charge  separation  in  a  thundercloud  are  not 
fully  understood.  Generally,  thunderclouds  are  formed  in  an 
atmosphere  with  layers  of  cold  dense  air  at  higher  levels  and  warm 
moist  air  at  lower  levels,  giving  rise  to  an  unstable  condition. 
This  instability  results  in  the  warm  air  rising  in  a  strong  updraft, 
which  in  the  presence  of  a  temperature  gradient  and  the  gravitational 
field,  is  responsible  for  the  formation  of  the  thundercloud.  Recent 
reviews  of  the  proposed  cloud  electrification  processes  have  been 
presented  by  Magono  (1980)  and  Latham  (1981). 

The  electrical  configuration  of  a  typical  thundercloud  is  shown 
in  Figure  1.1.  Typically,  the  upper  portion  of  the  thundercloud 
carries  a  net  positive  charge  (the  P  region),  whereas  there  is  a 
preponderance  of  negative  charge  in  the  lower  portion  (the  N 


Figure  1.1.     Electrical  configuration  of  a  typical  thundercloud; 

the  bold  dots  represent  the  locations  of  the 

effective  charge  centers.  (Adapted  from  Mai  an, 
1963;  and  Uman,  1969) 


region).  In  addition  to  these  main  charge  centers,  there  is  a  small 
pocket  of  positive  charge  (the  p  region)  at  the  base  of  the  cloud. 
The  basic  electrical  structure  of  the  cloud  was  first  determined  from 
in-cloud  measurements  by  Simpson  and  Scrase  (1937).  More  recent 
studies  reported  by  Huzita  and  Ogawa  (1976)  and  by  Winn  et  al .  (1981) 
have  confirmed  the  validity  of  this  basic  electrical  structure, 
although  it  is  now  recognized  that  there  may  be  large  horizontal 
displacements,  of  the  order  of  kilometers,  between  the  P  and  N 
regions. 

1.2.2  Ground  Lightning 

A  typical  cloud-to-ground  discharge  commences  in  the  cloud,  and 
eventually  results  in  the  neutralization  of  tens  of  coulombs  of 
negative  cloud  charge.  There  is  some  disagreement  on  the  exact 
nature  and  location  of  the  process  responsible  for  the  origin  of  the 
ground  discharge  (e.g.,  Clarence  and  Malan,  1957;  Rustan  et  al . , 
1980;  Beasley  et  al..  1982).  However,  it  is  believed  that  the 
discharge  is  initiated  by  a  preliminary  breakdown  within  the  cloud 
probably  between  either  the  N  and  p  regions  or  the  P  and  N  regions 
shown  in  Figure  1.1.  This  preliminary  breakdown  sets  the  stage  for 
electrons  from  the  N  region  to  be  funnel ed  away  from  the  cloud  and 
towards  the  ground  in  a  series  of  faintly  luminous  steps.  This  phase 
of  the  discharge  is  the  stepped  leader.  Photographic  observations  of 
stepped  leaders  give  the  typical  duration  as  1  us  for  steps  40  to  80 
m  in  length,  with  a  pause  time  between  steps  of  about  50  usee  The 
total  stepped  leader  process  typically  takes  about  20  ms  and  results 
in  the  distribution  of  about  5  C  of  negative  cloud  charge  along  5  km 


of  channel  length  from  cloud- to-ground;  the  average  downward  velocity 
IS  about  2.5  X  10  m/sec.  The  stepped  leader  electric  field  records 
indicate  a  smooth  transfer  of  charge  from  the  cloud  onto  the  channel, 
indicating  that  the  stepping  process  itself  is  not  associated  with 
any  appreciable  lowering  of  charge. 

When  the  leader  gets  within  a  few  hundred  to  a  few  tens  of 
meters  of  the  ground,  the  electric  field  beneath  it  is  high  enough  to 
cause  the  initiation  of  upward-moving  discharges  from  objects  on  the 
ground,  thereby  commencing  the  attachment  process.  The  attachment 
process  is  of  considerable  practical  interest,  especially  to  power 
engineers,  since  it  determines  the  point  at  which  the  lightning 
strikes  the  transmission  or  distribution  lines.  The  height  of  the 
leader  above  ground  at  the  time  when  the  attachment  process  begins  is 
defined  as  the  "striking  distance"  of  the  lightning,  a  parameter 
which  has  assumed  importance  in  the  study  of  direct  lightning  strikes 
to  power  lines.  The  concept  of  striking  distance  has  been  discussed 
in  detail  in  Chapter  17  of  Golde  (1977b). 

When  one  of  the  upward-moving  discharges  contacts  the  downward- 
moving  leader,  typically  some  tens  of  meters  above  ground,  it 
initiates  the  return  stroke  phase.  The  first  return  stroke 
propagates  up  the  previously  ionized  leader  path,  and  is  typically 
associated  with  a  peak  current  of  30  kA  at  ground  level.  The  return 
stroke  wavefront  travels  up  the  channel  at  about  one-third  the  speed 
of  light,  making  the  trip  from  the  ground  to  the  cloud  in 
about  100  ys,  neutralizing  some  or  all  of  the  charge  stored  on  the 
stepped  leader  channel.  The  rapid  release  of  the  return  stroke 
energy  heats  the  channel  to  about  30,000°K,  generating  a  high- 


pressure  channel  which  expands  against  the  atmospheric  pressure. 
This  rapid  expansion  creates  a  shock  wave  which  eventually  becomes 
thunder. 

After  the  return  stroke  current  has  ceased  to  flow,  the  flash 
may  end.  However,  if  any  additional  charge  is  made  available  at  the 
top  of  the  channel  through  J-  or  K-processes,  a  dart  leader  may 
propagate  down  the  residual  channel.  A  typical  dart  leader  travels 
continuously  down  the  channel  at  an  average  velocity  of  3  x  10^ 
m/sec,  and  lowers  about  1  C  of  charge  into  the  channel.  The  dart 
leader  then  initiates  a  subsequent  return  stroke,  the  leader  return 
stroke  combination  occurring  three  or  four  times  during  a  typical 
flash.  Data  reported  by  Weidman  and  Krider  (1980)  indicate  that  a 
typical  subsequent  stroke  current  has  a  smaller  peak  current,  about 
15  kA,  faster  zero-to-peak  risetime,  but  similar  maximum  rate-of- 
change  of  current  as  a  typical  first  return  stroke. 

A  dart  leader  may  not  always  be  able  to  make  its  complete  trip 
from  cloud  to  ground  down  the  defunct  stepped  leader  channel.  In 
fact,  there  is  a  significant  percentage  of  leaders  which  begin  as 
dart  leaders,  but  become  stepped  leaders  close  to  the  ground.  These 
leaders  are  known  as  dart-stepped  leaders.  The  average  return  stroke 
initiated  by  a  dart-stepped  leader  would  be  expected  to  have  a  peak 
current  which  lies  between  the  30  kA  value  for  an  average  first 
stroke  and  the  15  kA  value  for  an  average  subsequent  stroke,  as 
discussed  above.  It  is  therefore  not  very  clear  whether  subsequent 
strokes  initiated  by  dart-stepped  leaders  should  be  classified  as 
first  strokes  or  as  subsequent  strokes.  Historically,  the  term 
"first  stroke"  has  been  used  to  denote  only  the  ^Qry   first  stroke  in 


a  flash,  and  statistical  distributions  for  various  parameters 
associated  with  return  strokes  have  been  investigated  under  this 
assumption.  However,  if  the  conversion  of  the  dart  leader  to  a 
stepped  leader  takes  place  sufficiently  high  in  the  cloud  the 
location  of  the  ground  strike  point  due  to  the  new  channel  may  differ 
from  that  of  the  previous  channel  by  distances  of  the  order  of 
kilometers.  The  TV  pictures  for  data  presented  in  this  thesis  for 
such  cases  show  completely  different  channels  below  the  Florida  cloud 
base,  typically  about  1  km.  It  is  generally  believed  that  the 
initial  portions  of  the  vertical  electric  and  horizontal  magnetic 
fields  produced  by  the  return  stroke  are  due  to  the  bottom  1  km  or 
less  of  the  cloud  to  ground  channel.  Furthermore,  a  new  channel  to 
ground  involves  a  new  attachment  process  which  is  of  prime  importance 
in  studies  involving  overhead  power  lines.  For  these  reasons,  in 
this  thesis  all  "subsequent"  strokes  which  exhibit  completely 
different  channels  on  the  TV  pictures  are  treated  as  first  strokes. 
However,  to  facilitate  comparisons  with  extant  data,  these  "first 
strokes"  will  be  clearly  labeled  when  data  are  being  presented. 

The  time  interval  between  strokes  is  typically  40  ms.  However, 
the  interstroke  interval  may  be  tenths  of  a  second  if  a  continuing 
current  flows  in  the  channel  after  the  stroke.  Continuing  currents 
represent  direct  charge  transfer  between  cloud  and  ground  and  are  of 
the  order  of  100  A.  From  a  study  of  five  storms,  Livingston  and 
Krider  (1978)  report  that  29%  to  46%  of  all  ground  flashes  in  Florida 
have  a  continuing  current  component. 

We  have  discussed  the  cloud  to  ground  discharge  as  being 
initiated  by  a  downward-moving  negatively  charged  leader.   A  ground 


discharge  of  this  type  is  called  a  negative  ground  flash.  A  typical 
negative  ground  flash  lasts  about  0.5  seconds  and  has  three  or  four 
return  strokes.  However,  cloud-to-ground  lightning  may  also  be 
initiated  by  a  downward-moving  positive  discharge.  Positive  ground 
flashes  are  rare  in  Florida.  Positive  flashes  have  larger  peak 
currents  and  charge  transfers  than  negative  flashes,  and  rarely  have 
multiple  strokes  (Berger  et  al . ,  1975;  Brook  et  al . ,  1982).  In 
addition,  cloud-to-ground  discharges  may  be  initiated  by  upward-going 
positively-  or  negatively-charged  leaders.  In  this  thesis  the 
negative  ground  flash  is  the  only  type  of  cloud- to-ground  discharge 
that  we  will  consider.  Hence  the  phrase  "ground  flash"  will  be  used 
to  mean  a  negative  ground  flash  initiated  by  a  downward  moving 
leader. 

1.2.3     Cloud  Discharge 

Cloud  discharges  occur  between  negative  and  positive  charge 
centers  in  the  same  cloud,  or  between  two  clouds  or  between  cloud  and 
the  surrounding  air.  A  typical  cloud  flash  lasts  about  0.5  msec,  has 
a  path  length  between  5  and  10  km,  and  results  in  the  neutralization 
of  10  to  30  coulombs  of  charge.  The  discharge  is  believed  to  be 
initiated  by  a  continuously  propagating  leader  which  generates  5  to  5 
weak  return  strokes  known  as  recoil  streamers  or  K-changes  when  the 
leader  contacts  pockets  of  charge  opposite  its  own.  Cloud  discharges 
have  not  been  as  extensively  studied  as  cloud  to  ground  discharges. 
In  the  context  of  power  line  protection,  cloud  discharges  have  always 
been  assumed  to  produce  negligible  effects.  However,  it  is  important 
to    keep     in    mind     that    ground    lightning    has     significant    in-cloud 
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components  and  hence  may  share  many  common  features  with  cloud 
lightning. 

1.3  Induced  Overvoltages 

Lightning  is  responsible  for  a  significant  percentage  of  all 
disruptions  in  the  electrical  power  transmission  and  distribution 
network.  In  fact,  the  details  of  the  disruptive  mechanism  have  been 
the  subject  of  research  for  almost  eighty  years,  as  we  shall  discuss 
below. 

Lightning  overvoltages  on  power  lines  may  be  grouped  into  two 
categories:  (1)  voltages  caused  by  the  lightning  current  injection 
due  to  direct  strikes  to  phase  wires,  ground  wires,  or  towers;  and 
(2)  voltages  due  to  the  coupling  of  lightning  generated  electro- 
magnetic fields  emanating  from  nearby  lightning  which  does  not 
physically  contact  any  part  of  the  power  system,  often  referred  to  as 
an  indirect  strike  or  induced  effect.  In  the  present  thesis,  we 
consider  only  the  voltages  induced  on  the  line  by  nearby  lightning. 

Perhaps  the  earliest  attempt  at  studying  lightning  induced 
voltages  on  overhead  lines  was  due  to  K.  W.  Wagner  (1908).  According 
to  him  the  electrical  charge  was  induced  by  a  thundercloud  on  an 
overhead  line  which  was  electrostatically  held  at  ground  potential 
due  to  leakage  over  insulator  strings  and  grounding  of  transformer 
neutrals.  When  the  cloud  was  discharged  by  lightning,  the  bound 
charge  on  the  line  was  released,  giving  rise  to  traveling  waves  on 
the  line. 
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Adendorff  (1911)  presented  an  excellent  discussion  on  the 
various  power  line  problems  created  by  lightning.  Using  some  wery 
simple  models,  he  calculated  the  potential  difference  between  cloud 
and  ground  to  be  about  2.6  x  lO'^  V  for  a  cloud  height  of  5  miles.  He 
also  noted  that  overhead  ground  wires  solidly  earthed  to  ground  at 
the  poles  would  greatly  reduce  the  effects  due  to  a  direct  strike  on 
the  line,  the  effect  being  "entirely  preventative  and  not 
curative."  The  following  quotation  on  page  166  of  Adendorff  (1911) 
reflected  the  consensus  view  at  that  time  vis-a-vis  direct  lightning 
strikes  to  lines:  "If  the  discharge  is  very  heavy,  as  usually  is  the 
case,  the  probabilities  are  that  a  portion  of  the  section  struck  will 
completely  disappear  .  .  .  ." 

In  the  earlier  years  of  research,  it  was  believed  that  the  power 
line  should  be  protected  from  induced  surges,  the  implication  being 
that  not  much  could  be  done  in  the  case  of  direct  strikes.  However, 
as  power  transmission  expanded,  so  did  the  disruption  and  damage  due 
to  direct  lightning  strokes.  Thus,  it  is  not  surprising  that  during 
the  late  1920s  and  1930s  various  theoretical  and  experimental  studies 
were  undertaken  to  study  the  mechanism  and  effects  of  direct 
lightning  surges  on  power  transmission  lines,  e.g..  Cox  et  al . 
(1927),  Fortescue  et  al .  (1929),  Cox  and  Beck  (1930),  Lewis  and  Foust 
(1930),  Bewley  (1931),  Peek  (1931),  Rendell  and  Gaff  (1933),  Perry 
(1941),  and  Perry  et  al .  (1942). 

As  the  understanding  of  direct  strikes  to  lines  increased,  and 
lines  were  adequately  protected  against  them,  attention  was  again 
shifted  to  the  consideration  of  induced  voltages.  The  most  signif- 
icant of  the  early  research  on  induced  overvoltages  is  due  to  Wagner 
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and  McCann  (1942).  They  considered  most  of  the  essential  features  of 
the  discharge  process  including  the  effects  due  to  the  return  stroke 
channel.  Using  very  simple  models,  they  showed  that  induced  effects 
due  to  cloud  flashes  were  negligible,  whereas  induced  voltages  due  to 
close  ground  lightning  have  an  appreciable  magnitude. 

Szpor  (1948)  calculated  the  induced  voltage  on  the  line  caused 
by  a  nearby  lightning  stroke.  However,  according  to  his  own 
estimates,  the  quasi-static  nature  of  his  solution  restricted  the 
applicability  of  his  results  to  points  within  about  100  m  of  the 
lightning  discharge. 

A  comparison  of  the  frequencies  of  occurrence  of  surges  due  to 
direct  strikes  and  those  due  to  induced  effects  was  made  by  Golde 
(1954)  for  various  types  of  overhead  lines.  He  used  more  refined 
models  for  his  processes:  (a)  channel  charge  was  exponentially 
distributed  by  the  leader;  (b)  charge  neutralization  by  the  return 
stroke  did  not  take  place  instantaneously,  but  took  some  finite  time 
due  to  the  charge  stored  on  the  corona  sheath;  and  (c)  velocity  of 
the  return  stroke  was  assumed  to  decrease  exponentially  with 
height.  The  detailed  theory  used  in  his  calculations  was  presented 
in  an  earlier  paper  by  Bruce  and  Golde  (1941). 

Lundholm  (1957)  developed  explicit  expressions  for  induced 
voltages  based  on  the  simple  model  described  by  Wagner  and  McCann 
(1942).  He  used  the  concept  of  retarded  potentials,  thereby  avoiding 
the  quasi-static  solution  due  to  Szpor  (1948).  He  derived  independ- 
ent expressions  for  the  scalar  potential  in  terms  of  the  charge  on 
the  channel  and  the  vector  potential  in  terms  of  the  assumed  channel 
current.   However,  his  final  results  for  these  potentials  fail  to 
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satisfy  the  Lorentz  Condition  as  dictated  by  his  equations  and  are 
hence  in  error. 

Rusck  (1958)  developed  expressions  for  the  induced  voltage  on 
multi -conductor  lines.  He  found  that  the  induced  voltage  on  any  one 
conductor  was  not  influenced  by  the  presence  of  other  conductors--a 
surprising  result.  Further,  his  theoretical  analysis  followed  that 
due  to  Lundholm  (1957)  and  was  therefore  in  error. 

In  Wagner  and  McCann  (1942),  the  authors  argue  that  the  voltages 
induced  on  an  overhead  power  line  by  all  phases  of  ground  lightning, 
except  the  return  stroke,  are  insignificant,  because  only  "during  the 
progress  of  the  return  streamer,  the  rates  of  field  change  with  time 
are  sufficient  to  induce  a  voltage."   Owa  (1964)  calculated  the 
induced  voltage  on  a  power  line  due  to  a  stepped  leader.   Using  the 
stepped  leader  mechanism  proposed  by  Gri scorn  (1958)  and  Gri scorn  et 
al.  (1958),  which  gives  the  stepped  leader  current  as  100  kA,  and 
using  a  theoretical  method  of  solution  parallel  to  that  given  by 
Lundholm  (1957)  and  Rusck  (1958)  ,  he  found  peak  induced  voltages  of 
about  1  MV  produced  by  the  stepped  leader  for  a  line  height  of  about 
10  meters.   However,  currents  derived  from  measured  electric  field 
waveforms  by  Krider  and  Radda  (1975)  and  Krider  et  al .  (1977) 
indicate  that  the  peak  stepped  leader  current  is  of  the  order  of  1 
kA.    Further,  measurements  of  induced  voltages  on  power  lines 
reported  in  the  literature  have  never  approached  the  order  of 
magnitude  calculated  by  Owa  (1964),  indicating  that  the  theory  of 
Griscom  (1958)  is  in  error.    The  first  data  on  simultaneously 
recorded  stepped  leader  vertical  electric  field  and  induced  test-line 
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voltage  are  presented  in  Appendix  B  and  do  not  support  the 
theoretically  deduced  values  of  Owa  (1964). 

A  mathematical  determination  of  the  induced  voltage  on  a  power 
line  was  the  highlight  of  the  Ph.D.  dissertation  of  Chowdhuri 
(1966).  Contrary  to  the  findings  of  Rusck  (1958),  he  reported  that 
the  induced  voltage  on  a  conductor  of  a  multi -conductor  line  was 
influenced  by  the  presence  of  other  conductors.  Chowdhuri  and  Gross 
(1969)  reported  that  induced  voltage  was  higher  in  the  presence  of 
other  conductors  and  could  be  determined  from  a  consideration  of 
their  mutual  coupling.  However,  his  derivation  of  the  lightning 
fields  followed  that  due  to  Lundholm  (1957)  and  was  therefore  also 
incorrect. 

Another  theoretical  treatise  on  induced  voltages  on  overhead 
lines  was  presented  by  Singarajah  (1971).  He  used  the  old  model  of 
Wagner  and  McCann  (1942)  but  included  the  effect  of  the  upward-going 
streamer  on  the  return  stroke.  He  plotted  the  induced  voltage 
waveforms  from  his  equations,  using  various  parameters.  He  found 
that  the  voltage  on  any  one  conductor  was  not  affected  by  the 
presence  of  other  insulated  conductors  of  a  multi -conductor  line. 
However,  the  presence  of  a  grounded  earth  wire  was  found  to  reduce 
the  induced  voltage  on  the  line  by  about  ^0%   to  25%. 

An  experimental  program  initiated  by  the  South  African  National 
Electrical  Engineering  Research  Institute  (NEERI)  involves  the 
measurement  of  lightning  overvoltages  on  a  specially  constructed, 
unenergized,  9.9  km,  3-phase,  11  kV  distribution  line.  However,  no 
attempts  have  been  made  to  simultaneously  record  the  electric  and 
magnetic  fields  due  to  lightning.  A  description  of  the  test-line  and 


15 


some  of  the  preliminary  results  on  measured  overvoltages  were 
presented  by  Eriksson  and  Meal  (1980).  They  reported  measurements  of 
100  oscillograms  with  a  median  value  of  25  kV  for  voltages  induced  on 
the  line  due  to  nearby  lightning.  These  measurements  were  obtained 
as  triggered  records  with  a  threshold  of  12  kV.  In  more  recent 
reports,  Eriksson  et  al.  (1982)  and  Eriksson  and  Meal  (1982)  give  a 
distribution  of  300  overvoltages  measured  on  their  line.  For  32 
"good  quality  overvoltage  recordings"  the  mean  was  45  kV  with  a 
standard  deviation  of  35  kV.  The  maximum  measured  surge  is  300  kV 
which  is  also  the  B_asic  Jnsulation  Uvel  (BIL)  of  the  line.  Typical 
measured  overvoltages  are  shown  and  comparisons  are  made  with 
theoretically  calculated  waveshapes  based  on  the  incorrect  derivation 
due  to  Rusck  (1958). 

In  1973,  a  new  facility  to  study  artificially  triggered 
lightning  was  built  by  various  cooperating  agencies  in  France.  A 
description  of  their  main  facility  at  St.  Privat  d'Allier  has  been 
given  by  Fieux  et  al.  (1978).  As  part  of  the  program  two  overhead 
lines  were  built:  one  a  telecommunications  line,  2.1  km  in  length, 
and  the  other  a  260  m  medium  voltage  power  line.  Both  lines  were 
terminated  at  either  end  in  resistances  equal  to  their  respective 
characteristic  impedances.  Hamelin  et  al.  (1979)  reported  that  for 
triggered  lightning,  an  average  peak  voltage  of  about  1  kV  was 
measured  on  the  telecommunications  line  at  one  end,  the  lightning 
being  1.4  km  from  the  other  end.  Recordings  made  at  one  end  of  the 
medium  voltage  line  have  an  average  peak  of  74  kV,  the  distance  to 
the  lightning  being  50  m  from  the  other  end.  The  analytical  result 
derived  for  the  electric  and  magnetic  fields  produced  by  a  vertical 
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dipole  is  essentially  correct.  However,  their  theoretical  analysis 
suffers  because  they  use  a  very  simplified  return  stroke  model 
current  with  a  double  exponential  waveshape.  More  recently,  attempts 
have  been  made  to  incorporate  the  effects  of  an  imperfectly 
conducting  ground  plane  into  the  theory  by  Leteinturier  et  al .  (1980) 
and  Djebari  et  al.  (1981).  Though  the  theory  presented  appears  to  be 
correct,  the  lightning  fields  are  computed  from  the  same  simple  but 
incorrect  model.  Further,  due  to  computation  problems,  the 
calculations  are  performed  only  for  times  greater  than  2  ys,  when  the 
field  variation  is  relatively  slow.  Thus,  computed  results  for  the 
first  two  microseconds  are  not  available,  although  the  fast  changes 
occurring  within  this  period  may  be  of  prime  importance  in  power-line 
coupling. 

Prior  to  1980,  all  the  theories  put  forth  to  explain  power-line 
voltages  induced  by  nearby  lightning  were  based  on  the  premise  that 
the  induced  line  voltage  is  dependent  only  on  the  vertical  electric 
field  due  to  lightning.  Since  the  vertical  electric  field  at  the 
ground  due  to  a  negative  return  stroke  is  negative  (see  Appendix  A), 
it  follows  that  the  induced  line  voltage  would  have  a  positive 
polarity.  Thus,  we  have  the  classical  assumption  that  induced  surges 
on  lines  have  a  positive  polarity.  On  the  other  hand,  if  a  negative 
return  stroke  terminates  directly  on  the  line,  it  produces  a  negative 
current  injection  and  hence  a  negative  voltage  surge.  Based  on  this 
dichotomy  alone,  measured  surges  on  single  phase  lines  have  often 
been  classified  into  direct  and  induced  surges. 

However,  results  from  a  recent  Japanese  study  presented  by  Koga 
et  al.  (1981)  show  measurements  of  induced  surges  of  negative 
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polarity  as  well.  Simultaneously  measured  waveforms  for  induced 
voltage  on  a  1  km  line  indicate  a  surge  of  negative  polarity  at  one 
end  and  another  of  positive  polarity  at  the  other  end.  The  theory 
presented  by  Koga  et  al .  (1981)  uses  the  incorrect  model  presented  by 
Chowdhuri  (1966).  However,  their  theoretical  derivation  of  the 
horizontal  electric  field  due  to  a  finite  earth  conductivity  appears 
to  be  essentially  correct  and  is  similar  to  the  formulation  which 
will  be  adopted  in  this  thesis.  Further,  lack  of  any  simultaneous 
records  of  lightning  fields  and  induced  voltages  correlated  with  the 
ground  strike  point  prevented  them  from  making  any  direct  comparisons 
between  their  theory  and  data. 

The  experimental  data  recorded  during  the  summer  of  1979  and 
presented  in  this  thesis  show  that  the  induced  voltage  on  the  460- 
meter  unenergized  test-line  may  be  of  positive  or  of  negative 
polarity  depending  on  the  location  of  the  ground  strike  point. 
Further,  it  will  be  shown  theoretically  that  the  induced  voltage  on 
the  line  is  not  produced  by  the  vertical  electric  field  coupling 
alone:  the  horizontal  electric  field  resulting  from  a  finite  ground 
conductivity  must  also  be  taken  into  account,  since  the  voltage 
produced  by  it  may  dominate  that  due  to  the  vertical  electric  field. 

Coupling  of  incident  electromagnetic  fields  to  overhead  conduc- 
tors has  been  studied  by  various  groups  studying  the  effects  of  the 
Nuclear  Electromagnetic  Pulse  (NEMP)  on  overhead  conductors.  Formu- 
lations of  the  transmission  line  equations  involving  a  horizontal 
electric  field  have  been  discussed  by  Smith  (1977)  and  Vance  (197S) 
among  others.  A  complete  formulation  of  the  general  problem  in  the 
time  domain  has  been  presented  by  Agrawal  et  al.  (1980).   The 
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solution  scheme  outlined  by  Agrawal  et  al .  (1980)  will  be  used  in 
this  thesis.  The  complexity  involved  in  applying  this  type  of 
formulation  to  lightning  will  be  indicated. 


CHAPTER  II 
EXPERIMENTAL  RESULTS 

2.1  Introduction 

In  the  first  part  of  this  chapter,  a  complete  description  of  the 
experimental  setup  will  be  given.  The  systems  used  to  record  the 
vertical  electric  and  the  horizontal  magnetic  fields  produced  by 
lightning  will  be  described.  The  specifications  to  which  the 
unenergized  test-line  was  built,  the  apparatus  used  to  measure  the 
voltage,  and  the  various  systems  used  to  locate  the  ground  strike 
point,  or  points,  of  a  lightning  ground  flash  will  be  discussed. 

In  the  second  half  of  this  chapter,  data  on  induced  test-line 
voltages  will  be  presented  together  with  the  simulaneously  measured 
vertical  electric  field  due  to  lightning  return  strokes.  From  the 
analysis  of  the  vertical  electric  field  records  it  will  be 
established  that  the  lightning  return  stroke  sample  being  considered 
is  not  unusual  and  is  similar  to  other  observations  of  lightning, 
both  in  Florida  and  around  the  world.  Finally,  the  new  findings  with 
regard  to  the  measured  line  voltages  will  be  presented.  Data 
obtained  from  simultaneous  records  of  line  voltage  and  vertical 
electric  field  due  to  the  lightning  stepped  leader,  which  precedes 
the  first  return  stroke,  will  be  discussed  in  Appendix  3. 
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?.2  Experimental  Environment 


The  Mobile  Lightning  ^Laboratory  (MLL)  was  housed  inside  a  12  m 
long  trailer  which  was  located  near  Wimauma,  in  the  Tampa  Bay  area  of 
Florida  and  southeast  of  Tampa.   The  _Unenergized  Test  Line  (UTL)  was 
located  160  m  to  the  south  of  the  MLL.  A  network  of  five  remotely 
switchable  TV  cameras  was  located  over  an  area  of  roughly  60  km^ 
around  the  MLL.   One  of  the  five  TV  cameras,  a  still  camera,  and  a 
dynamic  thunder-sensing  microphone  were  located  at  the  MLL  site.  The 
geographical  locations  of  the  TV  cameras  and  their  fields  of  view, 
the  MLL,  and  the  UTL,  are  shown  in  Figure  2.1.   For  measurement 
accuracy  and  personal  safety  a  "true"  ground  plane,  60  m  in  diameter, 
was  built  from  a  wire  mesh  and  buried  in  the  ground  beneath  the 
MLL.   At  various  points,  3  m  ground  rods  were  sunk  vertically  into 
the  ground,  with  a  15  m  rod  at  the  center  of  the  mesh.  The  effective 
resistance  of  this  ground  plane  was  less  than  1  Ohm.   The  MLL  was 
bonded  to  this  ground  plane,  as  was  the  distribution  transformer 
which  supplied  power  to  the  MLL.  A  shallow  trench  was  dug  over  the 
distance  from  the  UTL  to  the  MLL,  and  all  measurement  cables  were 
buried  in  this  trench.  Figure  2.2  shows  a  sketch  of  the  MLL  and  the 
UTL  to  scale,  and  Table  2.1  lists  some  of  the  important  parameters  of 
the  line. 

2.3  Data  Acquisition  Systems 

The  purpose  of  the  research  reported  in  this  thesis  was  not  to 
collect  voluminous  data  on  induced  line  voltages  but  rather  to 
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Figure  2.2.    Detailed  sketch  of  the  experimental  test-line. 
Important  parameters  are  listed  in  Table  2.1. 
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correlate  the  voltage  measurements  with  the  electromagnetic  field 
produced  by  the  lightning  together  with  the  location  of  the  ground 
strike  point.  Thus,  there  are  three  pieces  of  information  which  make 
up  the  complete  data  set:  the  fields  produced  by  lightning,  the 
induced  line  voltage,  and  the  location  of  the  ground  strike  point. 

2.3.1  Field  Measurements 

Simultaneous  measurements  of  the  lightning  fields  and  the 
induced  line  voltages  were  made  on  two  separate  systems:  one 
continuous  and  one  triggered.  A  schematic  diagram  of  the  complete 
recording  system  is  shown  in  Figure  2.3. 

The  electric  field  sensing  was  done  by  three  circular  flat-plate 
antennas  one  having  an  area  of  0.5  m^  and  two  an  area  of  0.2  m^ 
each.  The  antenna  outputs,  which  give  the  derivative  of  the  electric 
field,  were  electronically  integrated  and  recorded  either  on  35  mm 
film  or  on  an  analog  7-channel  Ampex  FR-1900  instrumentation  Tape 
Recorder  (ITR),  at  a  speed  of  120  ips.  The  0.5  m^  antenna  and  one  of 
the  0.2  m^  antennas  were  followed  by  integrators  and  amplifiers  to 
provide  four  channels  of  electric  field  signals  with  a  3  db  bandwidth 
from  0.03  Hz  to  1  MHz  and  an  80  db  dynamic  range,  covering  electric 
field  strengths  from  4  V/m  to  40,000  V/m.  These  signals  were 
recorded  in  R1  mode  on  channels  1  through  4  of  the  ITR,  with  a 
record/reproduce  bandwidth  from  dc  to  500  kHz  (-6  db).  The  third 
antenna  (area  =  0.2  m^)  was  also  followed  by  an  integrator  and 
amplifier  and  provided  two  channels  of  electric  field  signals  with  a 
3  db  bandwidth  from  160  Hz  to  2  MHz  and  40  db  dynamic  range.  These 
two  channels  were  normally  recorded  in  direct  mode  on  channels  5  and 
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6  of  the  ITR,  with  a  3  db  bandwidth  from  400  Hz  to  1.5  MHz.  However, 
for  the  purpose  of  this  experiment  the  electric  field  signal  was 
recorded  only  on  channel  5.  The  input  to  direct  channel  5  of  the  ITR 
was  connected  in  parallel  to  the  input  of  the  Biomation  805  Digital 
W_aveforni  Recorder  (DWR).  The  DWR  was  triggered  off  the  initial 
portion  of  the  return  stroke  field  and  delayed  by  a  pre-set  amount  in 
order  to  display  the  complete  waveform  on  two  different  gain  settings 
of  the  top  traces  of  two  Textronix  555  dual -beam  oscilloscopes.  The 
lower  beams  were  chopped  and  used  to  display  the  W-E  component  of  the 
magnetic  field  and  the  UTL  voltage,  as  discussed  later.  The 
oscilloscopes  were  fitted  with  cameras  loaded  with  film  which  was 
moved  vertically  past  the  screens  at  a  linear  speed  of  about  5  cm/s, 
perpendicular  to  the  beam  sweep  direction,  so  that  a  permanent, 
spatially-separated  record  of  the  waveforms  of  each  individual  return 
stroke  is  obtained.  The  oscilloscope  time  base  was  generally  set 
to  200  ys  full-scale.  The  effective  3  db  bandwidth  of  these  records 
was  from  160  Hz  to  2  MHz. 

Inputs  to  R-l  channels  1  through  4  of  the  ITR  were  connected  in 
parallel  to  channels  1  through  4  of  a  Gould  Brush  Strip  Chart 
_Recorder  (SCR)  with  a  bandwidth  from  dc  to  100  Hz.  The  chart  speed 
was  either  5  mm/s  or  25  mm/s.  Timing  and  thunder  records  were  also 
displayed  on  the  SCR,  as  discussed  later. 

The  sensing  of  the  magnetic  fields  due  to  lightning  was 
accomplished  using  a  system  of  two  vertical  orthogonal  loops  mounted 
on  top  of  the  MLL  and  di recti onally  aligned  to  sense  the  S-M  and  W-E 
components  of  the  horizontal  magnetic  field.  The  signals  from  the 
loop  antennas  were  integrated  and  amplified  with  a  3  db  bandwidth 
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from  1  kHz  to  about  1.5  MHz  and  were  normally  fed  into  two  DWRs. 
However,  for  the  purpose  of  this  experiment,  only  the  signal  due  to 
the  W-E  component  of  the  magnetic  field  was  fed  into  a  DWR  to  be 
digitized  and  recorded  on  one  of  the  channels  of  the  chopped  lower 
beam  of  the  oscilloscopes,  the  other  being  used  to  display  the  UTL 
voltage,  as  we  shall  discuss.  The  magnetic  field  signals  were  not 
recorded  on  either  the  SCR  or  on  the  ITR. 

2.3.2  Voltage  Measurements 

With  the  exception  of  the  north  and  south  end  poles,  the  460  m, 
5-span  UTL  was  identical  to  standard  Tampa  Electric  Company  (TECO) 
single-phase  7,620/13,000  V  grounded-wye  primary  lines.  In  Figure 
2.2  the  location  of  the  UTL  and  the  MLL  were  shown  in  detail,  and 
important  parameters  of  the  UTL  were  listed  in  Table  2.1.  The 
impulse  insulation  level  of  the  end  poles  was  made  in  excess  of  250 
kV  by  using  25  kV-class  insulators,  unbonded  wood  crossarms  and 
fiberglass  guy  insulators.  Provision  was  made  for  the  connection  of 
a  10  kV  lightning  arrestor  on  the  north-end  pole,  but  none  was  ever 
used.  A  protective  gap  of  40  cm  was  also  installed  at  the  north  end 
to  limit  the  maximum  lightning  impulse  voltage  on  the  line  to  about 
250  kV.  At  the  south  end  a  driven  ground  rod  and  a  downlead  were 
provided  with  plans  for  disconnecting  the  neutral  and  grounding  the 
phase  conductor.  To  measure  the  voltage  induced  in  the  top  (phase) 
conductor,  a  downlead  was  provided  at  the  north  end,  closest  to  the 
MLL.  This  downlead  was  terminated  on  the  top  plate  of  a  specially 
designed  parallel  plate  capacitor.  The  bottom  plate  of  this 
capacitor  was   mounted   on  a   square   ground   plane  with  6  m  sides,   which 
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was  in  turn  bonded  to  the  true  ground  plane  established  for  the 
MLL.  The  capacitor  plate  separation  was  0.5  m.  The  electric  field 
signal  in  the  capacitor,  generated  by  the  voltage  on  the  line,  was 
sensed  by  a  flat-plate  antenna,  \jery  similar  to  the  ones  used  for 
sensing  the  vertical  field  due  to  lightning,  built  into  the  bottom 
capacitor  plate.  This  signal  was  transmitted  through  a  coaxial 
cable,  buried  in  the  service  trench,  into  the  MLL  to  be  integrated, 
amplified,  and  then  recorded. 

For  the  purpose  of  our  experiment,  the  input  to  direct  channel  6 
of  the  ITR  was  the  UTL  voltage  signal.  In  addition,  for  wery  close 
storms  the  electric  field  signal  was  disconnected  from  the  most 
sensitive  FA  channel  4,  and  the  UTL  voltage  recorded  instead. 
For  distant  storms,  the  UTL  signal  was  sometimes  recorded  on  FM 
channel  1 . 

As  noted  previously,  the  line  voltage  was  also  displayed  on  the 
oscilloscopes.  The  input  to  channel  6  of  the  ITR  was  connected  in 
parallel  to  the  input  of  a  Biomation  805  DWR.  The  output  of  the  DWR 
was  connected  to  one  of  the  channels  of  the  chopped  lower  beams, 
displayed  at  two  different  sensitivities  on  the  oscilloscope  screens, 
and  permanently  recorded  on  the  film,  together  with  vertical  electric 
field  and  the  W-E  magnetic  field. 

In  instances  of  very  close  lightning,  the  line  voltage  was  also 
connected  to  channel  4  of  the  SCR  in  place  of  the  electric  field. 

2.3.3  Lightning  Location 

As  discussed  in  Section  1.2.2  the  ground  strike  point  of  a 
lightning  return  stroke  is  determined  only  when  the  leader  comes 


32 


within  tens  of  meters  of  the  ground  and  initiates  the  attachment 
process.  In  order  to  correlate  the  remotely  measured  fields  with  the 
lightning  characteristics,  it  is  important  to  know  the  ground  strike 
point.  With  this  goal  in  mind,  a  network  of  five  TV  cameras  was 
employed  to  record  lightning  in  the  vicinity  of  the  UTL.  Continuous 
videotape  records  were  obtained  during  active  thunderstorm  periods. 
The  field  of  view  of  each  TV  was  calibrated  in  degrees.  The  MLL  was 
equipped  with  a  dynamic  microphone  which  sensed  the  thunder  due  to 
lightning.  This  signal  was  amplified  and  recorded  on  channel  5  of 
the  SCR.  The  time  difference  between  the  electric  field  record  and 
the  thunder  from  an  individual  lightning  was  used  to  determine  the 
distance  between  the  MLL  and  the  ground  strike  point. 

The  south-east  corner  of  the  roof  of  the  MLL  was  equipped  with  a 
plastic  dome.  During  a  thunderstorm  "run"  an  observer  was  positioned 
under  this  dome,  and  a  visual  identification  of  the  lightning  and  its 
location  was  made.  The  observer's  comments  were  recorded  on  the 
audio  channel  of  a  videotape  recorder  located  at  the  MLL  and 
independently  on  an  Audio  C_assette  Recorder  (ACR).  Comments  judged 
particularly  important  were  also  written  on  the  SCR  paper  in  real 
time.  The  ACR  proved  very  valuable,  since,  during  a  run,  changes  in 
equipment  settings  were  sometimes  not  recorded  in  the  data  book  and 
only  a  verbal  record  on  the  ACR  was  available. 

If  a  particular  flash  was  seen  on  at  least  two  TVs,  triangu- 
lation  was  used  to  determine  the  location  of  the  ground  strike 
point.  The  time-to-thunder  and  the  observer's  comments  were  then 
used  as  a  check  on  the  location.  If  the  flash  was  seen  only  on  one 
TV,  the  thunder  record  was  used  to  locate  the  ground  strike  point. 
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and  the  observer's  comments  were  used  as  a  check.  The  ground  strike 
points  of  all  flashes  located  for  this  study  were  seen  on  at  least 
one  TV. 

2.3.4  Relative  Timing 

In  order  to  locate  a  ground  flash,  accurate  relative  timing  had 
to  be  maintained  between  the  various  pieces  of  data.   The  time 
signals  recorded  on  the  field  and  voltage  records  were  derived  from  a 
Jjme    Code  Generator  (TCG)  in  the  MLL.   The  TCG  was  periodically 
synchronized  to  WWV  time.  The  TCG  output  was  recorded  in  the  m   mode 
on  channel  7  of  the  ITR.   It  was  also  used  to  drive  light-emitting 
diodes  (LEOs)  positioned  on  the  edge  of  the  oscilloscope  screens  to 
be  recorded  on  film  and  to  provide  the  time  display  for  the  TV 
monitor  and  videotape  at  the  MLL.   Hence,  there  were  no  problems  of 
relative  timing  associated  with  the  records  on  ITR,  SCR,  film,  and 
the  videotape  recorded  by  the  TV  located  at  the  MLL  site.   However, 
each  of  the  four  remotely  switchable  TVs  was  equipped  with  its  own 
TCG  located  at  each  remote  site.   All  four  TCGs  were  periodically 
synchronized  to  the  WWV  time  signal.  However,  failures  of  individual 
TCGs  occurred  due  to  excessive  heat  and  humidity,  line  voltage 
surges,  and  circuit-breaker  operations,  so  that  there  were  instances 
when  all  TCGs  were  not  synchronized.   In  spite  of  these  problems, 
however,  the  correlation  of  most  of  the  TV  records  was  possible. 
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2.4  The  Data  Base 

The  data  presented  in  this  thesis  were  taken  from  three  separate 
thunderstorms  in  the  vicinity  of  the  MLL  and  the  UTL  in  the  Tampa  Bay 
area  during  the  summer  of  1979.  In  Section  2.3  it  was  pointed  out 
that  a  complete  dataset  consists  of  three  essential  elements:  the 
electromagnetic  field  due  to  lightning,  the  line  voltage,  and  the 
location  of  the  ground  strike  point.  From  an  initial  review  of  the 
records,  it  was  found  that  on  three  days  of  the  1979  storm  season, 
the  field  and  voltages  were  simultaneously  recorded,  and  the  TV 
network  was  operational.  Various  elements  of  the  dataset  were 
subsequently  extracted  from  recordings  made  on  those  three  days.  We 
now  describe  in  detail  the  process  of  data  selection  for  any 
particular  day. 

2.4.1  Flash  Selection 


The  first  step  in  the  data  extraction  was  to  make  a  general 
survey  of  the  field  and  voltage  records  to  determine  if  there  were 
good  data  on  the  film  and  on  the  instrumentation  tape.  When  it  was 
determined  that  good  waveforms  could  be  obtained,  all  the  TV  tapes 
were  scanned,  and  lists  of  times  and  angles  were  made  for  each  flash 
seen  on  the  screen.  The  next  step  involved  looking  through  the 
strip-chart  records  to  see  if  there  were  identifiable  thunder  records 
and  observer  comments  for  the  flashes  seen  on  the  TVs.  Hard  copies 
were  made  from  the  TV  tapes  for  each  flash  and  a  copy  of  the  thunder 
record  was  made.  The  flash  was  identified  by  the  jJniversal  Time  (UT) 
of  its  occurrence. 
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The  audio  cassette  recorded  during  the  storm  was  played  back  to 
get  an  idea  about  the  storm  in  general  and  to  note  any  changes  in 
equipment  gain  settings  not  recorded  in  the  data  book.  At  the  end  of 
this  phase  it  was  possible  to  tell  whether  the  ground  strike  point  of 
a  flash  could  be  located,  either  by  tri angulation  if  records  were 
obtained  from  more  than  one  TV,  or  by  using  one  TV  and  the  thunder 
record.  As  noted  previously,  all  the  flashes  selected  for  this  study 
were  seen  on  at  least  one  TV. 

2.4.2  Field  and  Voltage  Records 

The  next  step  in  the  analysis  was  to  look  through  the  simulta- 
neously triggered  field  and  voltage  records  on  film  and  make  hard 
copies  on  a  microfilm  copier.  These  records  were  all  on  a  time  base 
of  200  us  full-scale.  The  final  step  to  completing  the  database  was 
to  obtain  the  field  and  voltage  waveforms  from  the  continuous 
instrumentation  tape  records  and  photograph  them  on  polaroid  film. 
The  time  base  for  these  could  be  varied  from  200  ms  to  40  us  full- 
scale  depending  on  the  chosen  sampling  rate.  During  this  stage  of 
data  extraction,  some  of  the  flashes  which  could  be  located  were 
eliminated  from  the  dataset  because  of  the  absence  of  field  and 
voltage  records.  The  final  list  of  flashes  which  emerged  at  the  end 
of  this  phase  were  those  which  could  be  located  and  for  which  field 
and  voltage  records  were  available. 

2.4.3  Location  and  Calibration 

If  a  flash  was  seen  on  two  TVs,  the  angles  of  the  flash  from 
each  TV  were  measured  from  the  calibrated  screen,  and  tri angulation 
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was  used  to  find  the  location  of  the  ground  strike  point.  If  the 
time  to  thunder  at  the  MLL  was  also  available,  the  distance  of  the 
flash  from  the  MLL  was  estimated  by  assuming  that  the  velocity  of  the 
sound  in  that  environment  was  345  m/s.  This  distance  estimate  and 
the  observer's  comments  were  used  to  check  on  the  strike  point 
obtained  from  the  TV  tri angulation.  However,  if  the  flash  was  seen 
on  only  one  TV,  the  angle  from  that  TV  was  used  together  with  the 
distance  from  the  MLL  obtained  by  thunder  ranging,  to  locate  the 
ground  strike  point.  In  that  case,  the  observer's  comments  were  used 
as  a  check. 

Let  us  first  discuss  the  calibration  of  the  triggered  records 
obtained  from  the  film.  As  can  be  seen  from  Figure  2.3,  the  vertical 
electric  field  and  the  line  voltage  signals  obtained  from  the 
respective  antennas  were  integrated  and  passed  through  the  delay  line 
which  is  essentially  as  a  2  MHz  low-pass  filter  and  then  connected  to 
the  Biomation  DWRs.  The  outputs  from  the  DWRs  were  used  to  display 
the  waveforms  on  the  oscilloscope.  The  electric  field  antennas  were 
mounted  on  top  of  the  MLL.  This  resulted  in  enhancing  the  electric 
field  at  the  actual  antenna  plate  relative  to  the  vertical  electric 
field  at  the  ground,  which  led  to  the  addition  of  an  "enhancement 
factor"  in  all  the  field  waveform  calibrations.  Thus,  we  have  the 
following  calibration  formulae  for  the  film  records: 


a  X  F  X  S  X  2.6  X  K^ 
V R '2.,) 


a"<i'         V^l_  =  a  X  F  X  S  X  2.6  X  Ky  X  d  (2.2) 
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where, 

a  =  number  of  divisions  on  the  grid 
F  =  scope  sensitivity  (V/division  on  grid) 
S  =  biomation  sensitivity  (V  full  scale/V) 
2.6  =  attenuation  due  to  delay  line 
'^E(V)  ~  combined  calibration  constant  of  antenna  and  integrator 
system  (V/m  per  V) 
R  =  enhancement  factor 
d  =  parallel -plate  capacitor  spacing  (m) 

The  instrumentation  tape  calibrations  for  the  electric  field 
records  obtained  from  FM  and  direct  channels  are  given  by 


a  X  F  X  S  X  K^H 
"gr        R  ^'^'^' 


a  X  F  X  S  X  K.. 
and.        Eg^  =  ^ IHI  (2.4) 

where, 

'^FM(dir)  "  combined  calibration  constant  for  antenna  integrator, 

and  recording  system  (V/m  per  V) 
a,F,S,R  =  as  defined  in  Equations  (2.1)  and  (2.2) 

The  calibration  formula  for  the  voltage  records  obtained  from 
the  direct  tape  channel  is 


V^^  =  a  X  F  X  S  X  2.6  X  K  X  d  (2.5) 

which  is  exactly  the  same  as  Equation  (2.2). 
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2.5  An  Example 


We  shall  illustrate  the  process  of  gathering  the  pertinent  data 
for  this  study  by  discussing  in  detail  one  fairly  typical  flash.  The 
flash  being  discussed  occurred  at  190610  jJniversal  Time  (UT)  on  July 
15,  1979,  during  designated  run  #79196TR27  (run  #27  at  the  MLL  which 
was  recorded  on  day  196  of  the  year  1979). 

2.5.1  Location 


From  an  initial  survey  of  the  strip-chart  records  it  was  found 
that  at  190610  UT  there  was  a  ground  flash  with  possibly  two  or  three 
strokes  and  about  23  seconds  to  thunder.  The  observer's  visual 
identification  of  the  flash  written  in  real-time  on  the  strip-chart 
identified  the  location  at  125°  with  respect  to  Morth  (0°).  A 
reproduction  of  the  strip-chart  record  with  the  comments  made  in 
real-time  is  shown  in  Figure  2.4. 

The  TV  records  for  run  #  79196TR27  indicated  that  the  TV  camera 
located  at  the  MLL  registered  a  flash  with  two  strokes  at  190610 
UT.  From  the  calibrated  TV  screen  the  location  of  the  ground  strike 
point  was  found  to  be  at  120°  from  true  North.  One  of  the  remote 
TVs,  designated  TV2S,  registered  a  three-stroke  flash  at  190610  UT, 
with  a  bearing  from  TV2S  to  the  ground  strike  point  of  141°  from  true 
North.  Using  the  information  obtained  from  the  two  TVs,  an  inter- 
section point  was  defined  as  shown  in  Figure  2.5.  The  ground  strike 
point  determined  from  the  TV  intersection  is  located  120°  from  true 
North  at  a  distance  of  8.8  km  from  the  MLL.  The  120°  location  angle 
agrees    very   well    with   the   rough    visual    angle   of   125°    obtained    from 


Figure   2.4.  Six-channel    strip-chart   record   for   the  190610   flash 

during     run     #79196TR27,     showing     four  channels     of 

electric    field   signals,    one  channel    of  thunder,    and 
one  of  time-code. 
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the  observer's  comments.  The  distance  obtained  by  thunder  ranging  is 
found  to  be  8  km,  again  in  reasonable  agreement  with  the  location 
distance  obtained  from  the  TV  intersection. 

2.5.2  Film  Records 

At  190610  UT,  film  records  were  available  from  both  the 
oscilloscopes.  The  trigger  threshold  was  such  that  the  film  recorded 
two  strokes.  Reproductions  of  the  200  ys  full-scale  records  obtained 
from  both  the  films  are  shown  in  Figures  2.6(a)  and  (b).  The 
magnetic  field  for  the  first  stroke  is  saturated  on  both  records. 
Scales  for  the  vertical  electric  field  and  the  test-line  voltage  are 
also  shown  in  the  figure. 

2.5.3  Instrumentation  Tape  Records 

The  vertical  electric  field  records  for  the  190610  flash 
obtained  from  one  R1  and  one  direct  channel  of  the  ITR  are  shown  in 
Figure  2.7  on  a  scale  of  20  ms/div.  Records  from  the  direct  channel, 
in  the  vicinity  of  the  fast  field  changes,  at  time-scales  ranging 
from  40  us/div  to  4  us/div,  indicated  that  each  of  the  three  fast 
changes  was  a  return  stroke.  For  comparison  with  the  film  records, 
we  show  the  waveforms  of  the  vertical  electric  field  and  the  test- 
line  voltage  for  all  three  strokes  on  a  time  scale  of  8  ps/div  in 
Figure  2.8.  The  calibration  for  the  electric  field  from  the  first 
stroke  has  been  adjusted  for  saturation  inside  the  isolation 
amplifier  through  which  the  electric  field  signal  is  recorded  on  the 
direct  channel . 
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Figure  2.8.  Vertical  electric  field  and  test-line  voltages  for 
all  three  return  strokes  of  the  flash  obtained  from 
the  direct  channels  of  the  instrumentation  tape 
recorder  for  the  190610  flash  during  run 
#79196TR27.  The  electric  field  and  voltage  scales 
are  shown  for  each  stroke.  All  records  are 
80  PS  full -scale.  The  relative  displacement  in  time 
between  the  electric  field  and  the  voltage  records 
is  due  to  the  positioning  of  the  record  and 
reproduce  heads  of  the  Instrumentation  Tape 
Recorder. 
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2.5.4  Conclusions 

Comparing  Figures  2.6  and  2.8,  we  note  that  the  third  return 
stroke  in  the  flash  is  missing  from  the  film  records  because  it  was 
below  the  trigger  threshold.  If  continuous  data  were  not  taken  on 
the  ITR,  this  flash  may  have  been  characterized  as  a  two-stroke  flash 
rather  than  a  three-stroke  one.  Another  point  to  note  is  the  fact 
that  the  voltage  waveform  measured  from  the  direct  tape  channel  falls 
to  zero  much  faster  than  that  measured  off  the  film.  This  is  due  to 
the  poor  low-frequency  response  of  the  direct  channel.  For 
comparison  purposes,  we  show  in  Figure  2.9  voltage  waveform  records 
on  two  tape  channels,  one  m  and  one  direct,  for  the  single  stroke 
flash  at  192238  UT  during  run  #79199TR31.  As  noted  in  Section  2.3.1 
the  m  channel  has  a  frequency  response  from  about  dc  to  500  kHz  (-6 
db).  As  can  be  seen,  the  direct  channel  record  is  less  noisy,  but 
has  a  much  faster  decay  characteristic;  the  peak  voltage  measured  on 
each  record  is  essentially  the  same. 

In  Table  2.2  the  initial  peak  values  obtained  from  the  various 
film  and  tape  records  for  the  example  being  discussed  are 
summarized.  The  values  obtained  from  the  measurements  differ  by 
about  +10%,  well  within  the  desired  accuracy.  This  difference  is 
attributed  to  several  factors.  As  noted  in  Section  2.3,  the 
oscilloscope  records,  the  m  and  the  direct  channels  of  the 
instrumentation  tape  recorder  all  have  different  bandwidths.  In 
addition  to  the  differences  in  waveform  associated  with  the  frequency 
response,  there  is  error  in  measuring  small  deflections  (see,  for 
example,  Figures  2.6(a)  and  (b)),  particularly  in  the  presence  of 
noise. 
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2.6  Data  Analysis 

As  described  in  the  earlier  sections  of  this  chapter,  the  data 
base  accumulated  for  this  thesis  was  put  together  from  thunderstorms 
on  three  separate  days.  In  Tables  2.3(a),  (b),  and  (c)  we  list  the 
flashes  used  in  this  study  for  each  of  these  three  days.  This  table 
gives  the  location  of  the  strokes  from  the  MLL  in  polar  coordinates 
as  illustrated  in  Figure  2.5  and  the  time  interval  between  strokes. 
The  values  derived  for  the  initial  peak  vertical  electric  field  and 
the  peak  line  voltage  from  the  oscillograms  and  instrumentation  tape 
records  as  outlined  in  the  previous  section,  are  also  listed.  In 
addition,  from  available  40  ys  full-scale  instrumentation  tape 
records,  risetimes  for  the  electric  field  and  the  voltage  were 
measured  and  have  also  been  tabulated. 

Analysis  of  the  vertical  electric  field  records  will  first  be 
used  to  show  that  the  sample  of  selected  lightning  is  not  unusual. 
Pertinent  parameters  will  be  plotted  and  compared  with  statistics  on 
lightning  in  Florida  and  in  the  other  areas  of  the  world.  The  data 
obtained  for  first  strokes  will  be  analyzed  separately  from  that 
obtained  for  subsequent  strokes.  The  next  step  will  be  to  plot 
various  statistics  related  to  the  voltage  waveforms  and  compare  them 
with  available  data  from  other  studies.  Finally,  we  shall  illustrate 
the  relationship  between  the  peak  line  voltage  normalized  to  the 
initial  peak  vertical  electric  field  as  a  function  of  the  angle  of 
the  ground  strike  point  around  the  line. 
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2.6.1  Analysis  of  Vertical  Electric  Field  Records 

In  Figures  2.10  and  2.11  we  plot  histograms  of  the  location  of 
the  various  strokes  used  in  the  analysis.  In  Figure  2.10  we  plot  the 
number  of  strokes  at  various  angles  from  the  MLL,  0°  being  true 
Morth.  In  Figure  2.11  we  plot  the  number  of  strokes  occurring  at 
various  ranges  from  the  MLL.  The  shaded  data  in  Figures  2.10  and 
2.11,  and  in  all  other  figures  involving  first  stroke  parameters, 
represent  subsequent  strokes  which  have  a  new  channel  to  ground  and 
are  classified  as  first  strokes,  as  discussed  in  Chapter  I.  Figure 
2.10  shows  that  there  are  strikes  in  all  four  quadrants  around  the 
test-line,  although  there  is  relatively  less  lightning  in  the  first 
quadrant  (0°  -  90°)  and  the  fourth  quadrant  (260°  -  360°).  This  is 
partly  because  of  the  positioning  of  the  TV  network  and  partly  due  to 
operation  problems  associated  with  the  remote  TVs.  Figure  2.11 
indicates  that  we  have  data  on  strokes  at  distances  from  about  1  km 
to  about  20  km  from  the  line.  The  range  at  which  the  TVs  can 
register  a  flash  limits  the  range  of  the  measured  data  because  the 
data  set  was  limited  to  those  flashes  which  could  be  seen  on  at  least 
one  TV  monitor. 

A  total  of  90  flashes  have  been  analyzed  in  this  thesis.  In 
Figure  2.12(a)  we  plot  a  histogram  of  the  number  of  strokes  in  a 
flash;  Figure  2.12(b)  shows  a  cumulative  distribution  function 
derived  from  the  histogram.  The  average  number  of  strokes  per  flash 
is  3.8  with  a  standard  deviation  of  3.0;  the  median  value  is  3,  and 
the  most  probable  value  is  1.  A  histogram  of  the  time  interval 
between  successive  strokes  in  a  multiple-stroke  flash  is  shown  in 
Figure  2.13(a).  The  associated  cummulative  frequency  distribution  is 


Figure  2.10.  Number  of  strokes  in  the  analyzed  data  plotted  as  a 
function  of  the  angle  from  the  Mobile  Lightning 
Laboratory.  The  shaded  data  indicate  subsequent 
strokes  which  have  been  classified  as  first  strokes. 
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Figure  2.11.  Number  of  strokes  in  the  analyzed  data  plotted  as  a 
function  of  the  distance  from  the  Mobile  Lightning 
Laboratory.  The  shaded  data  indicate  subsequent 
strokes  which  have  been  classified  as  first  strokes. 
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Figure  2.12(a).  Histogram  showing  the  number  of  strokes  per  flash 
for  the  data  analyzed.  The  mean  is  3.8  with  a 
standard  deviation  of  3.0.  The  median  is  3. 
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Figure  2.12(b).  Cumulative  distribution  function  for  the  number  of 
strokes  per  flash  showing  the  correlation  between 
data  obtained  in  this  thesis  and  other  studies. 
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Figure  2.13(a).     Histogram     of     time     intervals     between     successive 
return  strokes  in  a  flash. 


85 


INTERSTROKE  TIME  INTERVALS 
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given  in  Figure  2, 13(b).  The  mean  for  this  distribution  is  90  ms 
with  72  ms  as  the  standard  deviation;  the  median  is  70  ms,  with  the 
most  probable  values  between  40  and  60  ms.  In  Figures  2.12(b)  and 
2.13(b)  the  data  presented  here  are  compared  with  those  given  by 
Schonland  (1956)  and  Thomson  (1980a). 

The  parameters  plotted  in  Figures  2.12(a)  and  (b),  and  in 
Figures  2.13(a)  and  (b)  were  obtained  from  continuous  instrumentation 
tape  records  which  were  examined  at  various  time-scales  ranging  from 
400  ms  to  40  us  full-scale.  Thomson  (1980b)  has  presented  a  summary 
of  data  on  the  number  of  strokes  per  flash  and  the  interstroke  time 
interval  measured  by  various  researchers  all  over  the  world  using 
various  techniques  which  include  records  of  fields  on  a  slow  time- 
scale,  optical  data,  and  measurements  made  from  closed-circuit  video 
tape  records.  The  technique  used  here  is  less  prone  to  error  than 
any  technique  used  before.  From  Table  1  in  Thomson  (1980b)  it  can  be 
seen  that  there  is  considerable  spread  in  the  values  of  the  mean 
number  of  strokes  per  flash  and  mean  interstroke  time  interval  as 
reported  in  the  literature.  The  mean  number  of  strokes  per  flash 
ranges  from  a  low  value  of  2.1  (Brantley  et  al.,  1975)  to  a  high  of  6 
(Kitagawa  et  al . ,  1962).  The  mean  interstroke  interval  ranges  from 
38  ms  (McCann,  1944)  to  150  ms  (Takeuti,  1965).  The  statistics 
presented  in  Figures  2.12(a)  and  2.13(a)  give  values  which  are  within 
the  range  of  those  presented  in  the  literature.  In  another  paper, 
Thomson  (1980a)  has  given  a  cumulative  distribution  function  for  the 
number  of  strokes  per  flash  for  80  lightning  flashes  in  Papua,  New 
Guinea.  Schonland  (1956)  also  gives  such  a  distribution  due  to  1800 
flashes  over  a  long  period  of  time  in  South  Africa.   In  Figures 
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2.12(b)  and  2.13(b)  these  two  distributions  are  also  redrawn  and 
indicate  reasonable  agreement  with  the  data  presented  in  this 
thesis.  Anderson  and  Eriksson  (1980)  report  a  mean  of  2.5  strokes 
per  flash  computed  from  638  flashes.  However,  over  half  the  number 
of  flashes  in  their  data  have  a  single  stroke.  They  also  report  a 
median  value  of  50  ms  for  their  interstroke  time  interval  distri- 
bution. A  summary  of  the  flash  statistics  listed  in  Table  2.4 
indicates  that  the  sample  of  flashes  analyzed  in  this  thesis  is 
fairly  typical  of  worldwide  lightning. 

We  shall  now  discuss  various  statistical  parameters  of  interest 
related  to  the  return  stroke  field  change.  We  analyze  first  and 
subsequent  strokes  separately.  As  discussed  in  Chapter  I,  in  this 
thesis  we  define  a  "first"  stroke  as  one  which  is  associated  with  a 
new  channel  to  ground,  as  seen  on  the  TV  monitors.  However,  as 
mentioned  earlier,  the  data  points  associated  with  these  strokes  will 
be  shaded  in  the  histograms,  so  that  they  may  be  readily 
identified.  There  are  22  such  events  out  of  a  total  of  112  first 
strokes. 

The  initial  peak  in  the  vertical  electric  field  at  ground  is 
essentially  due  to  the  radiation  component  of  the  total  electric 
field  and  can  be  identified  readily  at  all  but  very  close  (under  500 
m)  distances  from  the  lightning  channel.  Figure  2.14(a)  shows  a  plot 
of  this  electric  peak  field  as  a  function  of  distance  from  the  MLL. 
This  initial  peak  may  be  normalized  to  a  distance  of  100  km  using  an 
inverse  distance  relationship  as  discussed  by  Lin  et  al .  (1979).  In 
Figure  2.14(b)  we  show  histograms  for  the  normalized  peak  for  112 
first  strokes:  the  mean  is  5.2  V/m  and  the  standard  deviation  is  3.4 
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Table  2.4.  Summary  of  flash  statistics  from  selected  studies. 

NUMBER     MEAN      S.D.     MEDIAN 

STROKES  PER  FLASH 

This  thesis  (Figure  2.12(a))  90 

Thomson  (1980a)  80 

Brantley  et  al .  (1975)  119 

86 

Kitagawa  et  al.  (1962)  87 

Schonland  (1956)  1800 


INTERSTROKE  TIME  INTERVAL  (ms) 

This  thesis  (Figure  2.13(a))   244 
Thomson  (1980a) 
Brantley  et  al .  (1975) 

Takeuti  (1965) 
Schonland  (1956) 
McCann  (1944) 
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Figure  2.14(b).  Peak  vertical  electric  field  of  the  return  stroke 
normalized  to  100  km.  For  112  first  strokes  the 
mean  is  5.2  V/m  with  a  standard  deviation  of  3.4 
V/m;  without  the  shaded  data  the  mean  for  90  first 
strokes  is  6.5  V/m  with  a  standard  deviation  of  3.5 
V/m.  For  237  subsequent  strokes,  the  mean  and 
standard  deviation  are  3.8  V/m  and  2.2  V/m, 
respectively. 
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V/m.   Eliminating  the  shaded  data  points  leads  to  a  mean  of  6.5  V/m 
with  a  standard  deviation  of  3.5  V/m  for  90  first  strokes.   For  237 
subsequent  strokes  the  mean  is  3.8  V/m  with  a  standard  deviation  of 
2.2  V/m.   Tiller  et  a1 .  (1976)  reported  a  mean  of  9.9  V/m  with  a 
standard  deviation  of  6.8  V/m  for  75  first  strokes;  for  163 
subsequent  strokes  the  mean  and  standard  deviation  were  reported  to 
be  5.7  V/m  and  4.5  V/m  respectively.  McDonald  et  al.  (1979)  studied 
two  night  storms  in  Florida,  one  over  the  Gulf,  the  other  over  the 
Atlantic.   The  Atlantic  storm  had  a  vertical  electric  field  peak 
distribution  which  gave  a  mean  and  standard  deviation  of  5.4  V/m  and 
2.1  V/m,  respectively,  for  54  first  strokes,  and  3.6  V/m  and  1.3  V/m, 
respectively,  for  119  subsequent  strokes.   Much  larger  means  were 
reported  for  the  Gulf  storm:  10.2  V/m  for  first  strokes  and  5.4  V/m 
for  subsequent  strokes.   Taylor  (1963)  reported  a  mean  value  of  4.8 
V/m  for  47  first  strokes  measured  in  Oklahoma.   Lin  et  al .  (1979) 
have  presented  the  normalized  values  of  the  initial  vertical  electric 
field  as  measured  simultaneously  at  two  stations,  one  within  15  km  of 
the  lightning,  and  the  other  either  50  km  or  200  km  away.   Their 
reported  mean  values  at  the  close  station  for  first  and  subsequent 
strokes  are  6.4  V/m  and  4.7  V/m,  respectively.   They  show  that  the 
average  normalized  value  measured  at  the  distant  station  for  the  same 
data  is  smaller  than  that  measured  at  the  close  station  as  a  result 
of  attenuation  due  to  propagation  over  lossy  earth.   To  summarize, 
the  values  for  the  vertical  electric  field  peak  normalized  to  100  km 
given  in  this  thesis  agree  very  well  with  the  extant  data. 

We  now  focus  on  the  various  parameters  associated  with  the  rise 
of  the  initial  vertical  electric  field  to  its  peak  value.   For  the 
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sake  of  consistency  all    risetime  measurements   reported  in  this  thesis 
have    been    made    from    the    direct    channel    records    obtained    from    the 
instrumentation    tape   on   a   time-scale   of  40   \s  full-scale.     When   such 
records  were   not  available,    no   risetime  measurements  were  made.     The 
error    involved    in    determining    the    zero    is    estimated    to    be    about 
±0.5   ys  in   most   cases.      A   further   error  in   measurement  is   estimated 
at  about  ±0.1   ys.      In   Figure  2.15  we  present  histograms  for  the  zero- 
to-peak    risetimes    of   the   vertical   electric    field.      The   distribution 
for   105   first  strokes  has  a  mean   of  4.4   ys  with  a   standard   deviation 
of  1.8   ys;  for      220  subsequent      strokes      the     mean      is      2.8   ys  with 
1.5   ys  as    the  standard  deviation.      Eliminating   the   shaded  data   gives 
a   mean   of  4.6   ys  with    1.8   ys  as    the    standard   deviation   for  84   first 
strokes.       In    a    Pennsylvania    study,    Fisher    and    Uman    (1972)    found    a 
mean        zero-to-peak        risetime       from       200   ys  full-scale        records 
of  3.6   ys  with    a    standard    deviation    of  1.8   ys  for    25    first   strokes; 
for     26     subsequent     strokes     the     mean     was  3.1    ys  with     a     standard 
deviation    of  1.9   ys.      All    these    strokes    were    within    25    km    of    the 
measuring    station.      For    storms    in   the   area   around   the   Kennedy   Space 
Center,    Florida,    Lin   and   Uman    (1973)    using    100   ys  full-scale   records 
reported     a     mean     zero-to-peak     risetime     of  4.0   ys  with     a     standard 
deviation   of  2.2   ys  for    12    first    strokes.      For   83    subsequent   return 
strokes,      the      mean      and      standard      deviation      were      reported      to 
be  1.2   ys  and  1.1    ys,     respectively.        Data     reported    for    two    close 
storms    around    Gainesville,     Florida,     by    Tiller    et    al .     (1976)    give 
a  mean       zero-to-peak      risetime      measured      from      200   ys  full-scale 
and  100   ys     full-scale    records    for    120    first    strokes   as    3.3    ys  with 
1.0   ys  as    the    standard    deviation;    for    163    subsequent    return    strokes 


Figure  2.15.  Distribution  of  the  zero-to-peak  risetime  of  the 
vertical  electric  field.  For  105  first  strokes  the 
mean  and  standard  deviation  are  4.4  us  and 
1.8  ys,  respectively;  for  220  subsequent  strokes, 
2.8  us  and  1.5  us,  respectively.  Without  the  shaded 
data  the  mean  for  84  first  strokes  is  4.6  us  with  a 
standard  deviation  of  1.8  us. 
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the  mean  zero-to-peak   risetime  was  reported  as  2.3   us  with  a  standard 
deviation     of  0.9  ys.       For     the     two-station     electric     field     data 
reported    by     Lin    et    a1 .     (1979)     and    recorded    on    time-scales    from 
40   ys  to  200   ys  full-scale,   the  mean  vertical    electric  field  zero-to- 
peak    risetime   as   measured   at  the   close   station   for  80   first  strokes 
was   found  to   be  2.5   ys  with  a   standard   deviation   of  1.3   ys.     For  142 
subsequent    strokes    the    mean    and    standard    deviation    were    found    to 
be  1.6   ys  and  0.8  ys,   respectively.        For    the    same    data,     the    mean 
risetimes    measured    at    the    far    station    are    much    larger    due    to    the 
degradation    in    the    front    due    to    propagation.        To    conclude,    the 
distributions    of    zero-to-peak    risetimes    for    first    and    subsequent 
return    strokes    analyzed    in    this    thesis    are   in    agreement  with    those 
reported   elsewhere    in    the    literature.       Differences    in    the    reported 
values  can  be  attributed  to  the  different  time  scales  of  measurement. 
It  is  difficult  to  define  the  exact  zero  field  value,   especially 
due   to   the   tape   recorder  noise.     The  10%-90X  risetime  is   therefore  a 
more      meaningful      parameter      to      measure      than      the      zero-to-peak 
risetime.        In    Figure    2.16    we    present    histograms    of    the    vertical 
electric     field     10^-90%     risetimes.         For    105     first    strokes,     the 
distribution    of    the    vertical     electric    field    10%-90X    risetime    was 
found      to      have      a      mean      of  2.5   ys  with      a      standard      deviation 
of  1.2   ys;   for      220      subsequent      strokes,      the     mean     and     standard 
deviation        were        found        to        be  1.5   ys  and  0.9   ys,   respectively. 
Eliminating  the   shaded  data,   we   find  a  mean  of  2.7    ys  with  a   standard 
deviation  of   1.2   ys  for  84  first  strokes.     There  appear  to  be  no  data 
on    the    10%-90%    risetime    of    the    vertical    electric    field   reported    in 
the    literature.       However,    the    values    shown    in    Figure    2.15    are    in 


Figure  2.15.  Distribution  of  the  10%-to-90%  risetime  of  the 
vertical  electric  field.  For  105  first  strokes  mean 
and  standard  deviation  are  2.6   us  and 

1.2  ys,  respectively;  for  220  subsequent  strokes, 
1.5  us  and  0.9  ys,  respectively.  Without  the  shaded 
data  the  mean  is  2.7  ts  with  a  standard  deviation  of 
1.2   us  for  84  first  strokes. 
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conformity  with  the  data  on  zero-to-peak  risetimes  reported  in  this 
thesis  and  elsewhere,  as  discussed  in  the  previous  paragraph. 

Weidman  and  Krider  (1978)  have  suggested  that  the  initial 
wavefront  of  the  vertical  electric  field  due  to  a  return  stroke  may 
be  decomposed  into  two  parts:  a  slow  initial  front,  followed  by  a 
fast  rise  to  peak.  The  fast  portion  of  the  wavefront  is  responsible 
for  the  maximum  rate-of-change  and  hence  it  is  a  parameter  of  prime 
importance  in  all  engineering  studies.  From  a  study  of  125  strokes 
propagating  about  50  km  over  sea  water,  Weidman  and  Krider  (1980) 
have  reported  a  mean  of  90  ns  and  a  standard  deviation  of  ^0  ns  for 
the  10%-90%  risetime  of  the  fast  portion  of  the  wavefront. 
Interestingly,  they  report  that  both  first  and  subsequent  return 
stroke  fields  have  about  the  same  maximum  rate-of-change.  However, 
the  fast  portion  of  the  wavefront  is  severely  degraded  as  a  result  of 
propagation  over  land.  For  example,  they  found  that  the  mean  10%-90% 
risetime  of  the  fast  transition  for  29  strokes  which  struck  sea  water 
over  a  distance  of  10-35  km  but  which  propagated  over  only  3  km  of 
land  was  201  ns,  more  than  twice  the  mean  obtained  when  there  was  no 
propagation  over  land.  In  Figure  2.17  we  plot  histograms  for  the 
10%-90%  risetime  of  the  fast  transitions  for  data  presented  in  this 
thesis.  The  mean  is  970  ns  and  the  standard  deviation  is  680  ns  for 
102  first  strokes.  Without  the  shaded  data  points  the  mean  is  1010 
ns  with  a  standard  deviation  of  710  ns  for  82  first  strokes.  For  217 
subsequent  strokes,  the  mean  and  standard  deviation  are  610  ns  and 
270  ns,  respectively.  All  propagation  in  our  data  is  over  land  with 
distances  ranging  from  about  1  km  to  20  km  as  shown  in  Figure  2.11. 


Figure  2.17.  Distribution  of  the  10%- to- 90%  risetime  of  the 
fast  transition  in  the  vertical  electric  field. 
For  102  first  strokes  the  mean  and  standard 
deviation  are  0.97  ps  and  0.68  ys,  respectively; 
for  217  subsequent  strokes,  0.61  us  and 
0.27  ys,  respectively.  Without  the  shaded  data  the 
mean  is  1  ys  with  a  standard  deviation 
of  0.70  ys  for  82  first  strokes. 
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In  Figure  2.18  we  plot  a  histogram  for  the  duration  of  the 
slow  initial  front  of  the  vertical  electric  field,  for  105  first 
strokes.    The  mean  and  standard  deviation  were  computed  to  be 
2.9  PS  and  1.3  us,  respectively.  Removing  the  shaded  data,  we  obtain 
a  mean  of  3.0  jjs  and  a  standard  deviation  of  1.3  ys  for  83  first 
strokes.    Weidman  and  Krider  (1978)  have  reported  a  mean  of 
about  4  MS  for  152  first  strokes,  which  is  close  to  the  mean  found  in 
this  thesis.   Another  parameter  of  interest  is  the  amplitude  of  the 
slow  front  expressed  as  a  percentage  of  the  associated  return  stroke 
peak  vertical  electric  field.  From  Figure  2.19  we  see  that  the  mean 
for  the  data  presented  in  this  thesis  is  28  percent  for  105  first 
strokes,  the  standard  deviation  being  15  percent.  Elimination  of  the 
shaded  data  leads  to  a  mean  of  27  percent  and  a  standard  deviation  of 
15  percent  for  83  first  strokes.  Weidman  and  Krider  (1978)  reported 
a  mean  of  about  44  percent  with  a  standard  deviation  of  20  percent 
for  34  first  strokes.   It  should  be  borne  in  mind  that  the  data 
presented  in  Figures  2.17,  2.18,  and  2.19  are  for  return  strokes 
which  propagate  over  ground  for  distances  ranging  from  1  km  and  20 
km.  The  data  presented  by  Weidman  and  Krider  (1978)  were  for  strokes 
propagating  over  30  to  50  km  of  sea  water.   Hence,  the  relatively 
large  values  of  the  risetime  associated  with  the  fast  field  change  of 
the  vertical  electric  fields  measured  in  this  thesis  can  easily  be 
understood. 

A  summary  of  the  return  stroke  vertical  electric  field 
statistics  is  given  in  Table  2.5.  Comparisons  with  statistics 
reported  in  the  literature  indicate  that  the  lightning  return  strokes 
studied  in  this  thesis  are  not  unusual. 


Figure  2.18.  Distribution  of  the  duration  of  the  slow  initial 
front  in  the  first  stroke  vertical  electric  field. 
For  105  first  strokes  the  mean  is  2.9  ps  with  a 
standard  deviation  of  1.3  ys;  without  the  shaded 
data  the  mean  and  standard  deviation  are  3.0  ys  and 
1.3  ys,  respectively  for  83,  first  strokes. 
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Figure  2.19.  Distribution  of  the  first  stroke  slow  initial  front 
amplitude  as  a  percentage  of  the  peak  vertical 
electric  field.  For  105  first  strokes  the  mean  is 
28%  with  a  standard  deviation  of  15%;  without  the 
shaded  data,  the  mean  is  27%  with  a  15%  standard 
deviation  for  83  first  strokes. 
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Table  2.5.     Summary  of  return  stroke  vertical   electric  field 
statistics  from  selected  studies. 


FIRST  STROKES 

SUBSEQUENT  STROKES 

N 

MEAN 

S.D. 

N 

MEAN 

S.D. 

Vertical  Electric  Field 

peak  normalized  to  100  km  (V/m) 

This  Thesis  (Figure  2.14(b)) 

112 

6.2 

3.4 

237 

3.8 

2.2 

Lin  et  al.  (1979) 

80 

6.4 

— 

142 

4.7 



McDonald  et  al .  (1979) 

54 

5.4 

2.1 

119 

3.6 

1.3 

52 

10.2 

3.5 

153 

5.4 

2.6 

Tiller  et  al.  (1975) 

75 

9.9 

6.8 

163 

5.7 

4.5 

Taylor  (1963) 

47 

4.8 



___ 

__. 

105 

4.4 

1.8 

220 

2.8 

1.5 

80 

2.5 

1.3 

142 

1.6 

0.8 

120 

3.3 

1.0 

163 

2.3 

0.9 

12 

4.0 

2.2 

83 

1.2 

1.1 

25 

3.6 

1.8 

26 

3.1 

1.9 

Vertical   Electric  Field  zero- 
to-peak  risetime   (ys) 

This  Thesis   (Figure  2.15) 

Lin  et  al.    (1979) 

Tiller  et  al .    (1976) 

Lin  and  Uman   (1973) 

Fisher  and  Uman   (1972) 

Vertical   Electric  Field 
10%-90%  risetime  (ys) 

This  Thesis   (Figure  2.16)  105       2.6       1.2         220       1.5       0.9 

Vertical   Electric  Field 
Fast  Transition 
10%-90%  risetime  (ns) 

This  Thesis   (Figure  2.17)  102       970       680         217       610       270 

Weidman  and  Krider  (1980)  125        90         40 

Vertical   Electric  Field 
slow  ramp  duration   (ps) 

This  Thesis   (Figure  2.18)  105       2.9       1.3 

Weidman  and  Krider  (1978)  152       4.0       1.7 


Table  2.5--continued. 
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FIRST  STROKES    SUBSEQUENT  STROKES 


N   MEAN  S.D.    N   MEAN  S.D. 


Vertical   Electric  Field  - 

Slow  Ramp 

amplitude  as  percent  of  peak 

This  Thesis  (Figure  2.19)    105   28    15 

Weidman  and  Krider  (1978)    34    40    20 
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2.6.2  Analysis  of  Voltage  Records 

Theories  of  the  induced  voltage  mechanism  as  presented  by  Wagner 
and  McCann  (1942),  Golde  (1954),  Lundholm  (1957),  Rusck  (1958). 
Chowdhuri  (1960),  and  Singarajah  (1971)  make  the  assumption  that  the 
line  voltage  is  induced  only  by  coupling  of  the  vertical  electric 
field  produced  by  nearby  lightning.  Thus,  induced  line  voltages 
produced  by  the  nearby  negative  return  strokes  to  ground  would  always 
be  of  positive  polarity,  as  discussed  in  Chapter  I.  Since  in 
previous  studies  the  inducing  electric  or  magnetic  fields  produced  by 
lightning  were  never  measured  in  the  vicinity  of  the  line  on  which 
the  voltage  measurements  were  made,  the  validity  of  the  above 
assumption  was  never  put  to  test.  We  will  now  present  some  of  the 
interesting  aspects  related  to  the  induced  voltages  measured  in  this 
thesis.  In  particular,  we  shall  show  that  the  induced  voltage  may  be 
either  positive  or  negative  in  polarity  for  the  same  value  of  the 
vertical  electric  field  peak,  the  polarity  and  amplitude  of  the 
voltage  being  a  strong  function  of  the  location  of  the  ground  strike 
point.  It  follows  that  the  coupling  cannot  be  due  to  the  vertical 
electric  field  only. 

In  Figure  2.20  we  plot  the  peak  measured  voltage  as  a  function 
of  distance  from  the  MLL  for  first  and  subsequent  strokes.  Peak 
voltages  induced  by  first  strokes  vary  from  +150  kV  to  -40  kV;  for 
subsequent  strokes,  the  values  range  from  +90  kV  to  -30  kV.  In 
Figure  2.21  we  show  the  distribution  of  the  absolute  value  of  the 
peak  induced  voltage  on  the  line.  For  112  first  strokes  we  find  a 
mean  of  22.6  kV  with  a  standard  deviation  of  22.4  kV.  For  237 
subsequent  strokes,   the  mean   and  standard  deviation  of  the 
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Figure  2.21.  Distribution  of  peak  induced  voltage  on  the  test 
line.  For  112  first  strokes  the  mean  and  standard 
deviation  are  22.6  kV  and  22.4  kV,  respectively;  for 
237  subsequent  strokes,  10.8  kV  and  9.0  kV, 
respectively. 


117 


PEAK  VOLTAGE  AMPLITUDE 


301 

or 

UJ  20 
DO 

:s 

3    10 
0 


E 


^ 


:2 


:^ 


2 


^ 


N=II2 

MEAN  =  22.6 
S.D=22.4 


Z2 


rr~i    I    I    1    i    ,    .    ,    ,    ,\\,    [— i    , 
0       10     20     30     40     50     60     70      80     90     100      150    160 

VOLTAGE  CKVD 


so- 
so - 

70 
60i 
LjJ  50- 

m 

^  40- 
30- 
20- 


10- 


N=237 

MEAN=I0.8 
S.D.=  9.0 


I    '     •    i^ — ' — 1 — I — I — I — I — I — I — R — r^  Vt — I — I — I 

0       10     20     30     40     50     60     70      80     90     100      150    160 


VOLTAGE  CKVD 


118 


distribution  were  10.8  kV  and  9.0  kV,  respectively.  Eliminating  the 
shaded  data,  we  find  a  mean  of  23.9  kV  with  a  standard  deviation  of 
24.4  kV  for  90  first  strokes.  Before  comparing  this  data  with  other 
studies,  we  take  one  more  look  at  Figure  2.20  and  note  that  most  of 
the  voltages  measured  in  this  thesis  were  induced  by  lightning  return 
strokes  which  were  between  about  4  and  12  km  distant.  Further,  the 
continuous  records  obtained  from  the  instrumentation  tape  recorder 
limit  the  lowest  measurable  voltage  to  the  channel  noise  level, 
estimated  to  be  about  1  kV. 

Measurements  of  overvoltages  on  the  specially  instrumented  9.9 
km  test  line  in  South  Africa,  reported  by  Eriksson  and  Meal  (1982) 
and  by  Eriksson  et  al .  (1982),  were  made  with  a  trigger  threshold  of 
about  12  kV.  The  most  common  type  of  surge  measured  in  that  study 
was  found  to  be  of  positive  polarity,  or  the  "classical  induced" 
type.  It  is  interesting  to  note  that  the  data  given  by  Eriksson  and 
Meal  (1982)  indicate  a  significant  number  of  positive  surges  (7  out 
of  21)  attributed  to  a  direct  strike  on  the  line.  Further,  a  small 
percentage  of  the  surges  exhibit  a  negative  polarity  and  are  not 
attributed  to  direct  strikes.  Out  of  a  total  of  281  surges  measured 
by  Eriksson  et  al .  (1982),  32  which  produced  "good  quality 
overvoltage  recordings"  were  located  accurately  by  tri angulation. 
Most  of  these  surges  were  produced  by  return  strokes  within  about  1 
km  from  the  line.  The  peak  voltage  for  these  32  records  has  a  mean 
of  45  kV  with  a  standard  deviation  of  35  kV,  and  there  does  not 
appear  to  be  any  strong  dependence  of  the  peak  voltage  on  the 
distance  from  the  line.  It  is  not  very  clear  whether  the  voltage 
recordings  are  always  due  to  the  first  return  stroke,  though  we  may 
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presume  that  to  be  the  case.  Voltages  induced  on  overhead  lines  due 
to  triggered  lightning  have  been  reported  by  Hamelin  et  al.  (1979). 
The  mean  peak  voltage  measured  on  the  2.1  km  telecommunication  line 
was  1  kV  for  lightning  1.4  km  from  the  line.  A  mean  of  74  kV  was 
reported  on  the  250  m  medium  voltage  line  for  lightning  50  m  away. 
We  conclude  that  the  peak  voltages  measured  on  our  test  line  are  much 
larger  than  the  values  reported  by  other  researchers. 

In  Appendix  D  we  present  a  histogram  of  measured  induced 
voltages  during  a  one  hour  and  thirty  minute  period  during  which  a 
storm  passed  very  close  to  the  line.  Since  the  TV  network  was  not 
operating  during  this  period,  no  flashes  were  located.  However,  some 
of  the  voltages  induced  on  the  line  during  this  period  were  due  to 
flashes  which  were  closer  to  the  line  than  those  presented  in  this 
thesis.  Figure  D.l  shows  several  voltages  greater  than  250  kV.  This 
does  not  appear  to  be  reasonable  since  the  line  was  designed  to 
spark-over  at  250  kV.  Thus,  either  the  spark-over  gap  was  improperly 
set  or  there  was  an  error  in  the  voltage  calibration  so  that  the 
actual  voltages  are  smaller  in  peak  amplitude  than  the  measured 
values  indicate. 

Figure  2.22  shows  the  distributions  of  the  zero-to-peak 
risetime  of  the  induced  voltage  waveforms  measured  in  this  thesis. 
All  risetime  measurements  have  been  made  on  a  time-scale 
of  40  MS  full-scale.  The  error  in  determining  the  zero  is 
estimated  at  ±0.5  us  in  most  cases.  A  further  error  in  measurement 
once  the  zero  is  defined  is  estimated  at  ±0.1  ms.  For  105  voltages 
induced  by  first  strokes,  the  mean  was  6.0  ms  with  a  standard 
deviation   of  3.8  us.    Without   the   shaded   data,   the   mean 


Figure  2.22.  Distribution  of  the  zero-to-peak  risetime  of  the 
induced  voltage  on  the  test-line.  For  105  first 
strokes  the  mean  and  standard  deviation  are 
6.0  MS  and  3.8  ys,  respectively;  for  218  subsequent 
strokes,  4.0   ys  and  2.3    ys,  respectively. 
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was  6.3  MS  with  4.0  ys  as  the  standard  deviation  for  84  first 
strokes.  For  218  voltages  induced  by  subsequent  strokes,  the 
mean  and  standard  deviation  were  found  to  be  4.0  vs  and 
2.3  us,  respectively.  Due  to  the  errors  involved  in  determining  the 
zero  and  the  peak  of  the  waveforms  in  the  presence  of  noise  it  is 
more  meaningful  to  talk  about  the  lO%-90%  risetimes  shown  in 
Figure  2.23.  For  105  first  stroke  induced  voltages  the  mean 
was  4.0  MS  with  3.2  ms  as  the  standard  deviation.  Without  the 
shaded  data  points,  the  mean  was  4.2  ms  and  the  standard 
deviation  3.4  ys  for  84  first  strokes.  For  218  subsequents,  the  mean 
and  standard  deviation  were  2.6  ms  and  1.7  ms,  respectively. 

Comparing  Figures  2.16  and  2.23,  we  can  conclude  that,  on  the 
average,  the  value  of  the  10%-90%  risetime  of  the  vertical  electric 
field  is  smaller  than  that  of  the  induced  voltage.  A  similar 
conclusion  may  be  drawn  with  regard  to  the  zero-to-peak  risetimes. 
However,  there  is  considerable  spread  in  the  ratio  of  the  10%-90% 
risetime  of  the  voltage  to  that  of  the  simultaneously  recorded 
vertical  electric  field.  This  distribution  is  shown  in  Figure 
2.24:  the  most  likely  value  of  the  ratio  is  between  one  and  two.  A 
summary  of  induced  voltage  statistics  presented  here  is  given  in 
Table  2.6. 

2.6.3  Correlation  Between  the  Voltage  and  Field  Records 

The  most  interesting  aspect  of  the  induced  voltages  on  the 
unenergized  test  line  measured  in  this  thesis  is  the  fact  that  the 
polarity  of  the  induced  voltage  was  found  to  be  either  positive  or 
negative  for  the  same  value  of  the  peak  vertical  electric  field. 


Figure  2.23.  Distribution  of  the  10%-to-90%  risetime  of  the 
induced  voltage  on  the  test-line.  For  105  first 
strokes  the  mean  and  standard  deviation  are 
4.0  iJS  and  3.2  us,  respectively;  for  218  subsequent 
strokes,  2.6  us  and  1.7  us,  respectively. 
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Figure  2.24.  Ratio  of  the  10%-to-90X  risetime  of  the  induced 
voltage  to  the  10%-to-90%  risetime  of  the 
simultaneously  recorded  vertical   electric  field. 


126 


10%  TO  90^  RISETIME  RATIO  OF  V/E 


a: 

UJ 

en 


20- 
15- 
10- 
5  - 
0 


FIRST 
N=105 
MEAN=I.8 
S.D=I.2 


M 


TlirfTipn 


f— ' 1 — 

7        8 


23456789 

RISETIME  RATIO 


10      II       12 


cr 

LU 

en 


30 

25  1 

20 

15- 

10- 

5  ■ 

0 


SUBSEQUENT 
N=2I7 
MEAN=2.2 
S.D.  =  I.9 


2       3       4 


Inn    nn.iiiiii        rg        n 
5       6        7       8       9       10      II 


RISETIME  RATIO 


.-a 

12 


127 


Table  2.6.  Summary  of  induced  voltage  statistics. 


FIRST  STROKES    SUBSEQUENT  STROKES 


N   MEAN  S.D.    N   MEAN  S.D. 


Test  Line  Voltage  peak 

amplitude   (kV)    (Figure  2.21)  112      22.6     22.4        237      10.8      9.0 

Test  Line  Voltage  zero-to-peak 

risetime  (ps)    (Figure  2.22)  105       6.0       3.8         218       4.0       2.3 

Test  Line  Voltage  10%-90% 

risetime   (ys)    (Figure  2.23)  105       4.0       3.2         218       2.6       1.7 

Ratio  of  10%-90%  risetime  of  V 

to  10%-90%  E  risetime  of  E 

(Figure  2.24)  105       1.8       1.2         217       2.2       1  9 
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This  IS  contradictory  to  the  various  theories  associated  with  the 
concept  of  the  classical  induced  surge.  In  Figure  2.25  this  point  is 
illustrated  with  the  help  of  three  subsequent  return  strokes  which 
are  at  approximately  the  same  distance  from  the  line.   The  located 
position  in  each  case  is  indicated  in  polar  coordinates.  Each  of  the 
three  return  strokes  has  about  the  same  value  of  the  peak  vertical 
electric  field  peak  as  measured  at  the  MLL,  the  measurement  of  line 
voltage  being  made  at  the  northern  end  of  the  line.  As  can  be  seen 
clearly  from  the  figure,  the  line  voltages  induced  by  essentially  the 
same  vertical  electric  field  are  remarkably  different  for  the  three 
return  strokes.   For  the  stroke  due  north  of  the  line,  the  induced 
voltage  waveform  is  of  negative  polarity.   For  the  stroke  due  south 
of  the  line,  the  voltage  is  of  positive  polarity  and  has  about  the 
same  peak  amplitude  as  that  induced  by  the  stroke  due  north  of  the 
line.    For  the  return  stroke  with  its  ground  strike  point 
perpendicular  to  the  line,  the  induced  voltage  is  buried  in  the  tape 
recorder  noise  level  and  is  very  small  in  peak  amplitude  compared 
with  the  other  two  induced  voltages.   It  therefore  appears  from  the 
data  that  the  coupling  of  the  lightning  electromagnetic  field  on  the 
line  cannot  be  satisfactorily  described  in  terms  of  the  vertical 
electric  field  alone.   In  fact,  from  Figure  2.25  there  are  strong 
indications  that  the  dominant  driving  function  for  the  induced 
voltage  is  probably  the  horizontal  component  of  the  electric  field 
and  not  its  vertical  component.   The  horizontal  electric  field 
component  for  a  perfectly  conducting  ground  plane  is  directly  related 
to  the  horizontal  magnetic  field,  as  discussed  in  the  next  Chapter. 
Similar  conclusions  have  been  drawn  by  Koga  et  al .  (1981)  who  report 


Figure  2.25.  Induced  voltages  on  the  test  line  and  vertical 
electric  field  due  to  lightning  for  return  strokes 
at  three  different  locations. 
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surges  of  opposite  polarities  measured  simultaneously  at  two  ends  of 
a  telephone  cable.  However,  they  did  not  have  the  locations  of 
lightning  ground  strike  points  in  their  study. 

To  investigate  further  the  relationship  between  the  induced 
voltage  and  the  location  of  the  ground  strike  point  the  following  was 
done.  The  peak  value  of  the  induced  voltage  was  first  normalized  to 
the  simultaneously  recorded  peak  vertical  electric  field.   This 
normalization  was  done  in  an  attempt  to  remove  the  distance  and 
return  stroke  current  dependence.    The  value  of  V/E,  or  the 
normalized  voltage,  was  plotted  as  a  function  of  the  angle  around  the 
test  line.   The  strokes  are  grouped  together  in  10°  intervals.  The 
mean  and  standard  deviation  of  the  values  of  V/E  within  the  angle 
range  were  computed.   In  Figure  2.26  these  values  are  plotted  for 
first  strokes;  Figure  2.27  shows  the  relationship  for  subsequent 
strokes.   The  circled  point  represents  the  mean,  the  vertical  bars 
represent  the  value  of  the  standard  deviation  on  each  side  of  the 
mean,  and  the  number  of  strokes  within  that  angle  range  is  shown  with 
each  point.   It  is  clear  from  these  figures  that  the  normalized 
voltage  V/E  has  a  strong  dependence  on  the  angle  around  the  line. 
The  values  for  first  strokes  are  not  significantly  dif-ferent  from 
those  for  subsequent  strokes.   This  is  not  surprising,  since  the 
voltage  has  been  normalized  with  respect  to  stroke  severity.   The 
curves  exhibit  a  positive  peak  at  around  180°,  go  through  zero  around 
90°  and  270°,  and  exhibit  a  negative  peak  around  0°. 

In  this  section,  we  have  examined  the  relationship  between  the 
lightning  return  stroke  field,  the  line  voltage,  and  the  location  of 
the  ground  strike  point.  We  have  found  that  the  induced  voltage  is 
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at  Us  maximum,  when  the  return  stroke  is  "end-on"  the  line,  the 
induced  effects  being  much  smaller  for  return  strokes  perpendicular 
to  the  line.   We  have  also  found  that  for  strokes  in  the  first  and 
fourth  quadrants  around  the  line,  the  induced  voltage  waveform  is  of 
negative  polarity;  strokes  in  the  second  and  third  quadrants  produce 
the  classical  induced  surge  of  positive  polarity.   Let  us  examine 
these  results.   It  appears  from  Figure  2.25  that  there  is  no  obvious 
relationship  between  the  induced  voltage  and  the  vertical  electric 
field,   since   essentially   the   same   vertical   electric   field 
illuminating  the  line  can  give  rise  to  widely  differing  voltages. 
Since,  for  a  negative  return  stroke,  the  vertical  electric  field 
always  points  into  the  ground,  the  horizontal  electric  field  produced 
by  the  finite  ground  conductivity  is  along  the  direction  of 
propagation  and  points  towards  the  lightning.   For  a  return  stroke 
due  north  of  the  line,  the  lightning  electromagnetic  field  would  be 
propagating  along  the  line  and  away  from  the  north  end  measuring 
point.    However,  for  a  stroke  due  south,  the  field  would  be 
propagating  along  the  line  and  toward  the  measuring  point.  For  both 
these  cases,  the  vertical  electric  field  along  the  line  would  be  in 
the  same  direction,  and  hence  if  the  vertical  electric  field  coupling 
were  dominant,  only  positive  polarity  voltages  would  be  produced. 
However,  the  component  of  the  horizontal  electric  field  along  the 
line  would  be  in  opposite  directions.   If  the  coupling  produced  by 
the  horizontal  electric  field  were  dominant,  it  would  explain  why 
voltages  of  either  polarity  are  induced.  This  would  also  explain  why 
the  induced  voltage  for  strokes  due  east  or  west  of  the  line  is 
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rather  small,  because  the  component  of  the  horizontal  electric  field 
along  the  line  would  be  very  small  in  that  case. 

Given  the  fact  that  negative  induced  surges  measured  on  the 
unenergized  test  line  for  this  thesis  are  almost  as  likely  to  occur 
as  positive  induced  surges,  why  is  it  that  measurements  of  induced 
surges  made  over  the  years    on    various  types  of  lines  indicate  a 
predominance  of  positive  waveforms?  To  answer  this  question,  we  note 
that  in  most  instances,  the  measurements  are  made  on  long  lines  which 
are  part  of  the  network  and  the  voltage  is  not  measured  at  the 
"end".   In  the  previous  paragraph  we  noted  that  the  induced  voltage 
would  have  a  negative  polarity  only  if  the  lightning  electromagnetic 
field  propagates  along  the  line  and  away  from  the  end  at  which  the 
voltage  is  measured.   If  the  measurement  is  made  in  the  middle  of  the 
line,  the  lightning  electromagnetic  field  is  always  propagating  along 
the  line  and  towards  the  measuring  station,  and  hence  the  induced 
surge  is  always  of  positive  polarity.  The  South  African  9.9  km  test 
line  was  equipped  with  one  measuring  station  at  one  end,  and  another 
measuring  station  in  the  middle.   However,  the  data  presented  by 
Eriksson  and  Meal  (1982)  and  by  Eriksson  et  al.  (1982)  are  drawn  only 
from  the  station  in  the  middle  of  the  line.   Hence,  it  is  not 
surprising  that  their  records  indicate  the  classical  induced  surge. 
The  experimental  setup  used  in  the  French  experiment  given  in  Hamelin 
et  al.  (1979)  is  such  that  only  positive  surges  would  be  measured. 
However,  results  from  the  recent  Japanese  study  on  a  telephone  cable 
shield  presented  by  Koga  et  al .  (1981)  and  discussed  earlier  show 
simultaneous  induced  surges  of  negative  and  positive  polarity 
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measured  at  the  two  ends  of  their  1  km  line.  This  is  in  keeping  with 
the  physical  arguments  presented  in  this  section. 

In  the  next  chapter  we  shall  mathematically  formulate  the 
problem  of  computing  induced  line  voltages  from  a  knowledge  of  the 
inducing  fields.  In  Chapter  IV  calculated  waveforms  will  be  compared 
with  all  available  data. 


CHAPTER  III 
THEORETICAL  ANALYSIS 

3.1  Introduction 

In  Chapter  I  the  various  theories  used  to  explain  induced  power- 
line  voltages  were  outlined,  e.g.,  Lundholm  (1957),  Rusck  (1958), 
Chowdhuri  (1956),  and  Singarajah  (1971).  These  theories  lead  to 
induced  surges  which  are  only  of  positive  polarity.  The  theory 
presented  by  Leteinturier  et  al.  (1930)  and  Djebari  et  al.  (1981), 
which  includes  the  effects  of  the  imperfectly  conducting  ground 
plane,  was  used  by  the  French  to  make  limited  comparisons  between 
their  theoretically  calculated  waveshapes  and  the  measured  waveform 
of  the  current  at  one  end  of  their  2.1  km  test  line.  It  was  pointed 
out  in  Chapter  I  that  the  major  drawback  of  their  theory  was  the  use 
of  a  very  simple  model  for  the  lightning  current.  Further,  their 
calcu-lations  were  performed  only  for  times  greater  than  2  us  when 
the  fields  do  not  vary  rapidly.  In  addition,  their  experimental 
setup  was  such  that  it  excluded  any  measurements  of  induced  surges  of 
negative  polarity.  No  general  calculations  were  done  which  predicted 
induced  line  surges  of  negative  polarity.  The  recent  Japanese  study 
reported  by  Koga  et  al.  (1981)  clearly  shows  measured  induced  surges 
of  negative  and  positive  polarity  on  their  test  line.  Although  some 
general  calculations  are  presented,  no  direct  comparison  between 
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their  theory  and  their  data  could  be  made  because  of  lack  of 
information  on  stroke  locations.  Furthermore,  they  do  not  present 
any  calculation  which  shows  negative  induced  surges  on  their  line. 

In  Chapter  II,  simultaneous  measurements  of  the  return  stroke 
vertical  electric  field  and  the  induced  voltage  on  a  specially 
constructed  unenergized  power  line  were  discussed.  It  is  clear  from 
the  data  that  the  "classical  induced"  surge  of  positive  polarity  is 
almost  as  likely  to  be  induced  on  our  460  m  test-line  as  a  surge  of 
negative  polarity.  Furthermore,  it  was  indicated  that  the  polarity 
of  the  surge  was  dependent  on  the  location  of  the  ground  strike 
point.  There  are  stronn  indications  from  the  data  that  the 
horizontal  electric  field  coupling  is  the  predominant  factor 
responsible  for  the  induced  voltage.  The  reasons  for  the  existence 
of  a  horizontal  electric  field  component  at  the  line  will  now  be 
investigated. 

(a)  Assuming  an  infinitely  conducting  ground  plane  and  a 
vertical  lightning  channel,  the  horizontal  electric  field  at  ground 
is  zero;  however  there  is  a  finite  horizontal  field  at  line  height, 
which  is  approximately  10  m  above  the  ground.  This  horizontal 
electric  field  may  be  predicted  using  the  return  stroke  model  of 
Master  et  al.  (1981).  However,  for  a  10  m  altitude  this  horizontal 
electric  field  is  too  small  to  produce  the  measured  voltages. 

(b)  More  often  than  not,  the  lightning  return  stroke  channel 
from  cloud  to  ground  is  not  perfectly  straight  and  vertical  as 
assumed  in  the  various  return  stroke  models,  but  exhibits  a  certain 
tortuosity  on  a  scale  from  about  1  m  to  over  1  km  (e.g.,  Evans  and 
Walker,  1963;  Hill,  1968).  A  dipole  with  an  arbitrary  orientation  in 
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space  would  result  in  different  electric  and  magnetic  fields  at  the 
line  height  compared  to  that  due  to  a  vertical  dipole  at  the  same 
point.  However,  data  analyzed  by  Hill  (1968)  indicate  that  the 
direction  of  sections  of  lightning  channels  are  randomly 
distributed.  For  segment  lengths  studied  (between  5  m  and  70  m)  and 
total  channel  lengths  analyzed  (between  1  km  and  4.3  km)  the  mean 
absolute  value  of  the  channel  direction  change  was  16°.  The  reported 
random  distribution  of  the  parameters  governing  channel  tortuosity 
indicate  that  the  net  effect  of  the  complete  lightning  channel  is 
such  that  it  can  be  modeled  as  vertical.  This  is  because  the  random 
distribution  of  tortuosity  cannot,  on  average,  provide  a  large  effect 
in  any  particular  direction.  It  is  in  fact  this  argument  which 
accounts  for  the  good  agreement  between  measured  data  and  theory  as 
presented  by  Lin  et  al.  (1979,  1980). 

Data  presented  on  the  measured  test-line  voltages  in  this  thesis 
for  lightning  at  some  fixed  distance  and  angle  from  the  MLL  do  not 
indicate  a  wide  range  of  variability.  This  further  indicates  that 
the  channel  shape  is  not  a  critical  factor  in  producing  the 
horizontal  electric  field.  In  fact,  the  averaged  data  on  the  peak 
line  voltage  normalized  to  the  vertical  electric  field  peak  as  shown 
in  Figures  2.26  and  2.27  indicate  very  little  spread,  and  hence  we 
argue  again  that  the  straight  vertical  channel  for  lightning  is  a 
good  approximation  to  an  actual  lightning  channel. 

(c)  Using  an  argument  that  the  ground  is  a  good,  but  not  a 
perfect  conductor,  a  horizontal  electric  field  may  be  computed  from 
the  vertical  component  at  ground,  based  on  the  wavetilt  of  the 
propagating  electromagnetic  wave.   This  method  gives  the  correct 
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waveshape  and  reasonable  amplitudes  for  the  theoretically  calculated 
voltages,  as  discussed  later. 

The  first  part  of  this  chapter  will  be  devoted  to  the 
investigation  of  the  horizontal  electric  field  produced  by  mechanisms 
(a)  and  (c)  outlined  above.  The  latter  half  of  the  chapter  will  be 
devoted  to  the  derivation  of  the  transmission-line  coupling 
equations,  and  a  discussion  of  their  numerical  solution.  It  will  be 
shown  that  the  total  voltage  at  any  point  on  an  overhead  line  may  be 
derived  from  a  knowledge  of  the  incident  horizontal  electric  field  at 
all  interior  points  and  the  incident  vertical  electric  field  at  the 
ends  of  the  line  and  at  the  measuring  point.  Listings  of  FORTRAN 
codes  used  to  obtain  the  various  solutions  are  given  in  Appendix  C. 

3.2  Electric  and  Magnetic  Fields  Illuminating  the  Test-Line 

3.2.1  The  Return  Stroke  Model 


The  word  "model"  as  used  in  this  thesis  may  be  defined  as  a 
mathematical  construct  of  the  properties  of  physical  phenomenon  which 
can  be  used  as  a  convenient  tool  to  further  the  understanding  of  that 
phenomenon.  The  lightning  return  stroke  has  been  subject  of  various 
models  for  over  four  decades.  An  important  objective  of  most  of 
these  models  is  to  be  able  to  describe  the  electric  and  magnetic 
fields  produced  by  the  return  stroke.  This  field  description  is 
arrived  at  from  a  knowledge  of  the  current  in  the  channel.  The 
modeling  of  the  current  itself  has  been  accomplished  to  varying 
degrees  of  sophistication,  e.g.,  Strawe  (1979),  Price  and  Pierce 
(1977),  Little  (1978),  Bruce  and  Golde  (1941),  Dennis  and  Pierce 
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(1964),  Uman  and  McLain  (1970),  Lin  et  al .  (1980),  and  Master  et  al. 
(1981). 

The  lightning  return  stroke  model  of  Lin  et  al .  (1980)  is  the 
only  available  model  consistent  with  measured  data  on  electric  and 
magnetic  fields.  The  model  currents  are  derived  from  data  recorded 
simultanously  at  two  recording  stations  on  the  ground,  one  within  15 
km  of  the  lightning  and  the  other  at  50  or  200  km.  The  model 
postulates  the  existence  of  the  following  three  current  components  in 
a  straight,  vertical  return  stroke  channel  above  a  perfectly 
conducting  earth: 

(a)  a  short  upward-propagating  pulse  of  current  associated  with  the 
electrical  breakdown  at  the  return  stroke  wavefront  and  responsible 
for  the  peak  current;  the  pulse  propagating  up  at  a  constant  velocity 
v; 

(b)  a  uniform  current,  which  may  already  be  flowing  in  the  channel 
as  leader  current,  or  may  begin  to  flow  soon  after  the  commencement 
of  the  return  stroke;  and 

(c)  a  "corona"  current,  caused  by  the  neutralization  of  the  charge 
stored  around  the  leader  channel,  and  discharged  by  the  passage  of 
the  return  stroke  wavefront. 

The  parameters  defining  these  three  current  components  are  derived 
from  the  measured  two-station  fields  and  are  discussed  in  detail  by 
Lin  et  al .  (1980). 

Although  the  model  of  Lin  et  al .  (1980)  was  primarily  applicable 
to  subsequent  return  strokes,  it  proved  reasonably  successful  at 
modeling  the  first  stroke  two-station  electric  and  magnetic  fields 
measured  at  ground  level  reported  by  Lin  (1978).  However,  Schonland 
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et  a1 .  (1935)  reported  that  the  peak  luminosity  of  the  first  return 
stroke  channel  decreased  with  height  above  the  ground.  Recent 
observations  reported  by  Jordan  and  Uman  (1980)  show  a  marked 
decrease  in  peak  luminosity  with  height  for  subsequent  strokes  as 
well.  These  observations  indicate  that  the  breakdown  pulse  current 
[(i)  in  the  previous  paragraph]  should  be  attenuated  as  it  propagates 
up  the  channel.  This  attenuation  would  not  substantially  alter  the 
electric  and  magnetic  fields  measured  at  ground  level;  however,  it 
would  significantly  alter  the  fields  measured  at  an  altitude.  This 
extension  of  the  model  of  Lin  et  a1.  (1980)  to  predict  fields  at 
altitudes  and  bring  the  model  currents  in  closer  conformity  with 
reported  observations  was  carried  out  as  an  integral  part  of  this 
thesis  and  was  reported  by  Master  et  a1 .  (1981).  This  model  can  be 
used  to  compute  the  electric  and  magnetic  fields  in  the  vicinity  of 
any  overhead  power  line.  The  fields  may  then  be  coupled  onto  the 
line  in  order  to  produce  the  induced  voltage. 

3.2.2  Electromagnetic  Field  Calculation 

In  Figure  3.1  we  indicate  the  geometry  of  the  return  stroke 
channel.  The  "image"  channel  will  be  used  to  simulate  the  effect  of 
the  perfectly  conducting  ground  plane  (see  Stratton,  1941,  pgs.  577- 
582).  Due  to  the  cylindrical  symmetry  of  the  problem  about  the 
return  stroke  channel,  the  solution  for  the  electric  and  magnetic 
fields  at  any  general  point  P  (r,  (*>,z)  will  be  obtained  using  a 
cylindrical  coordinate  system  with  the  origin  at  the  point 
at  which  the  return  stroke  makes  contact  with  the  earth.  The 
complete  channel  from  cloud  to  ground  is  split  into  small  channel 


Figure  3.1.     Geometry  for  field  computations  based  on  the  model 
of  Master  et  al.  (1981). 
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sections  or  current  dipoles.  The  electric  and  magnetic  fields  at 
point  P  (r,^,z)  produced  by  a  small  dipole  located  at  the  origin  and 
carrying  a  current  i(t)  have  been  derived  in  Appendix  A.  Equations 
(A. 35),  (A. 36)  and  (A. 37)  written  for  a  general  dipole  of  length  dz' 
and  at  height  z' ,  carrying  a  current  i(z',  t)  are  reproduced  below. 

pdz'        3i(z',t-^)   ^  _ 

B^(r,,,z,t)  =^[-^   ^T-^'-,      i(^'.t-§)]    (3.1) 


E^(r.^,z,t)  -  ^  l^^^^  f\u',  X  -'^-)dx  .  ^rUzlll  i(,.  ,t  -  i) 


3i(z',T-^) 


r(z-z')  °'^^  '  '  -  c 
:2r3 


'"^"-i^    3t-^^  (3.2) 


E,(r.*,z,t)=^[i(^I^:i!l!:!A-(z-,.-^)dx 


^0      R5 


R, 


?(z-zM2_r2  R      r2   9i(z'.t--) 

+  ii^A_Lj:Li(z',t  -h    --I -, ^] 

CR"^  ^     c2r3       3t 


(3.3) 


where,  y^  is  the  permeability  and  e^  the  permittivity  of  air  (or 
vacuum),  and  all  geometrical  factors  are  illustrated  in  Figure  3.1. 
The  current  i(z,t)  represents  the  currents  derived  from  the  models 
presented  by  Lin  et  al.  (1980)  and  Master  et  al.  (1981)  described  in 
Section  3.2.1.  The  first  terms  in  Equations  (3.2)  and  (3.3)  are 
generally  called  the  electrostatic  field;  the  second  terms  in 
Equations  (3.1),  (3.2)  and  (3.3)  are  known  as  the  induction  terms; 
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and  the  remaining  terms  which  involve  the  time-derivative  of  the 
current  are  called  the  radiation  terms. 

To  obtain  the  electric  and  magnetic  fields  at  P  (r,(j>,z)  due  to 
the  lightning  return  stroke,  we  have  to  sum  the  field  contributions 
due  to  each  individual  dipole  in  the  channel.  We  note  that  at  any 
time  the  complete  field  is  due  to  contributions  from  various  dipoles 
from  the  real  and  the  image  channels,  each  dipole  carrying  a 
different  current  which  depends  on  its  position,  on  the  channel,  and 
with  respect  to  the  field  point. 

The  computation  of  the  time-varying  electric  and  magnetic  field 
components  based  on  the  model  of  Master  et  al.  (1981)  is  performed  on 
a  digital  computer.  A  listing  of  the  FORTRAN  code  is  given  in 
Appendix  C.  As  input,  the  user  specifies  the  parameters  of  the 
return  stroke  model  current,  the  coordinates  of  the  field  point,  and 
the  total  time  for  which  the  field  description  is  desired.  The 
program  output  gives  the  electric  and  magnetic  field  components  and 
the  waveform  of  the  current  at  ground.  The  current  waveform  at 
ground  due  to  a  typical  subsequent  return  stroke  which  was  modeled  by 
Lin  et  al.  (19B0)  and  Master  et  al.  (1981)  is  shown  in  Figure  3.2. 
The  electric  and  magnetic  field  components  at  a  range  of  3  km  from 
the  channel  and  an  altitude  of  10  m  above  ground  produced  by  this 
current  are  shown  in  Figure  3.3. 

3.2.3  Wavetilt  Formulation 


If  the  earth  were  a  perfect  conductor,  a  propagating  electro- 
magnetic wave  would  have  no  tangential  electric  field  component 
at  the  earth's  surface.    However,  if  the  earth  has  a  finite 
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Figure  3.3.  Electric  and  magnetic  fields  produced  by  a  typical 
subsequent  return  stroke  at  a  distance  of  3  km  from 
the  channel  at  an  altitude  of  10  m. 
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conductivty,  a  propagating  electromagnetic  wave  has  a  non-zero  net 
tangential  electric  field  along  the  direction  of  propagation.  This 
is  because  the  tangential  component  of  the  propagating  magnetic  field 
at  the  ground  induces  a  current  density  in  the  ground,  which  for  a 
finite  conductivity  produces  a  tangential  electric  field  at  the 
surface.  The  classical  research  associated  with  the  understanding  of 
this  phenomenon  has  been  credited  mainly  to  Zenneck  (1907), 
Sommerfeld  (1909,  1926),  and  Norton  (1935,  1935,  1937).  An 
interesting  historical  perspective  is  given  in  Chapter  II  of  Wait 
(1962).  Special  mention  must  be  made  of  the  work  reported  by  Carson 
(1926)  which  was  concerned  with  propagation  in  overhead  wires  with 
ground  return.  For  a  moderately  good,  flat,  homogenous  ground,  the 
wavetilt  function  relating  the  Fourier  transform  of  the  vertical 
electric  field  to  that  of  the  horizontal  electric  field  at  ground  may 
be  adequately  approximated  by  (e.g.  Wait,  1962;  Vance,  1977) 

Eu(J'^)       1 

^(J'^^  =  r-rmv  =     ..  (3.4) 


where 


E|^  -  horizontal  electric  field 
Ey  -  vertical  electric  field 
W  -  wavetilt 


e  -  relative  permittivity  of  the  soil 

a  -  soil  conductivity 

e  -  permittivity  of  air  (or  vacuum) 


The  main  approximations  involved  in  arriving  at  the  above  formulation 
are  the  following:   (a)  the  incident  wave  is  a  radiation  field  with 
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its  electric  field  vector  in  the  plane  defined  by  the  direction  of 
propagation  and  the  normal  to  the  ground,  and  (b)  the  angle  of 
incidence  is  near  grazing  to  the  flat  earth's  surface.  Both  these 
approximations  are  reasonable  for  the  data  reported  in  this  thesis. 
Approximation  (a)  will  be  further  investigated  later. 

From  Equation  (3.4)  we  note  that  the  horizontal  electric  field 
may  be  obtained  from  the  vertical  electric  field  and  a  knowledge  of 
the  ground  parameters  e^  and  a.  For  a  radiation  field  in  air  (or 
vacuum),  the  vertical  electric  field  and  the  horizontal  magnetic 
field  are  related  by  the  free  space  impedance.  Due  to  this  fact. 
Equation  (3.4)  may  also  be  used  to  determine  the  horizontal  electric 
field  from  the  horizontal  magnetic  field. 

An  obvious  way  to  proceed  would  be  to  Fourier  Transform  the 
vertical  electric  field  and  then  obtain  the  Fourier  Transform  of  the 
horizontal  component  in  accordance  with  Equation  (3.4).  Inverse 
transforming  the  results  gives  E^(t).  For  any  real  waveform  this 
would  involve  using  the  Discrete  Fourier  Transform  (OFT)  and  the  Fast 
Fourier  Transform  (FFT)  technique.  A  more  direct  approach  will  be 
adopted  here.  Since  the  vertical  electric  field  may  be  obtained  in  a 
piecewise  linear  form,  the  ramp  response  of  W(s)  can  be  used  to 
obtain  the  horizontal  field  in  the  time-domain  without  resorting  to 
the  transform-inverse  transform  process. 

Let  us  consider  a  vertical  electric  field  signal  as  being 
represented  by  its  piecewise  linear  version  as  shown  in  Figure  3.4. 
Without  any  loss  of  generality  let  us  assume  that  there  are  only  four 
linear  segments  as  shown.  The  start  point  of  any  segment  j  is  given 
by  the  coordinates  (t.-,  Ew  )  and  the  end  point  by  (t^ , i ,  E,,   ). 

J      Vj  J+l      Vj+l 


Figure     3.4.  Example     of    a    piecewise    linear    vertical     electric 

field  due  to  a  return-stroke. 
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Thus,  an  m-segment  representation  will  be  defined  by  rm-l  points. 
Further,  it  will  be  assumed  that  outside  the  time  interval  given  by 
[t-| ,  t^+il  the  vertical  electric  field  is  zero.  Therefore,  from 
Figure  3.4,  an  analytical  expression  for  Ey(t)  may  be  written  down  as 
follows. 


Ev(t) 


^V2  ^2   ^3"  ^V2 


/t)  =  t^  •  t  .  U{t)  +  [-  i^  +   t,  -  tJ  •  ^^-^2)  •  U(t-t2) 


^  ^"   ta-  t2   ^   t.  -  t3^  •  (^-^3)  •  U(t-t3) 


^  ^-   t^-  t3   "   ts-  tJ  •  (t-t'*^  •  "(t-tO 


+  [-  ts-  t^  •  (^-^5)  -  Evs^  •  "(t-ts)        (3.5) 


where  u(t)  is  the  unit  step  function. 

It  can  be  seen  that  each  term  in  Equation  (3.5)  represents  a 
ramp,  except  for  the  last  term  which  is  a  step,  due  to  the  fact  that 
the  field  waveform  in  Figure  3.4  terminates  at  a  finite  non-zero 
value. 

At  this  point,  a  digression  is  necessary  in  order  to  determine 
the  ramp  and  step  response  of  W(s)  as  given  in  Equation  3.4.  For  a 
ramp  input,  i.e. 


we  have, 
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E,{s)   =  1 


/rr+~ojsT     s2 


Eh^s)     =  ^    •  -fTo ^ TTo  (3.7) 


/r       s3/2(s   +    a/e^e^)l/2 


Using    a    table    of    Inverse    Laplace    Transforms    (e.g.    Abramovitz    and 

Stegun,   1968),  we  have  from  Equation   (3.4) 

E^(t)   =-iL  .  t   .  e'P^  [iQ(pt)   +  I^(pt)]    .  u(t)               (3.8) 
^r 

where 

P            =    o/2e^e^ 

Iq,Ii    -    modified    Bessel's    Functions    of    the    first    kind    of 

the  zeroth  and  first  order  respectively. 

Similarly  for  a  step  input, 

Ey(s)   =  k  u(t)                                                                                (3.9) 

we  obtain 

E    ('■1    -                             ^                                                               ('<   ^^\\ 

"             /r  s^/^(s  +   a/e   £  )^/2 
r                      'or 

and  thus, 

E^(t)   =  A.e-P^  I^(pt)    .  u(t)                                              (3.11) 

r 
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To  summarize,  for  a  ramp  vertical  electric  field  given  by  Equation 
(3.6),  the  horizontal  electric  field  is  given  in  the  time-domain  by 
Equation  (3.8).  Similarly,  the  step  in  Equation  (3.9)  leads  to  a 
solution  given  by  Equation  (3.11).  Using  these  results  we  return  to 
the  waveform  shown  in  Figure  3.4  and  given  analytically  by  Equation 
(3.5).  For  that  vertical  electric  field,  the  horizontal  electric 
field  is  given  by. 


^2 

V^'>  =^  'rr  '^  '  ^'^^i^p^^i  +  ii(pt)  ]u(t) 

r 

1     ,    ^Vz      ^Va'  ^2^  -p(t-t2) 

r 

•     [lo(p(t-t2))    +    I^(p(t-t2))]u(t-t2) 
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^V3"^2      ^V.-^a  -p(t-t3) 
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1   ,  ^Vs  "  S,  -p{t-t5) 

.^[-^^^--].(t-t3)  .e 

r 


•  [I^ipit-ts))  +  I^(p(t-t5))]u(t-t5) 
1  ,  ,  -P(t-t5) 

+  -^  [-  E„  ]  •  e        .  I  (p(t-t5))  .  u(t-t5)    (3.12) 
r 

It  has  therefore  been  analytically  demonstrated  that  for  a  piecewise 
linear  vertical  electric  field  waveform  given  in  the  time-domain  by 
Equation  (3.5),  a  horizontal  field  waveform  may  be  derived,  also  in 
the  time-domain,  in  accordance  with  Equation  (3.12).  For  example, 
let  us  determine  the  horizontal  field  from  the  wavetilt  for  the 
vertical  electric  field  waveform  given  in  Figure  3.3.  In  Figure  3.5 
we  show  the  digitized  version  of  the  vertical  electric  field,  and  the 
derived  horizontal  electric  field  for  ground  parameters  given 
by  e^  =  10  and  a  =  10   mhos/m. 

One  of  the  assumptions  underlying  the  wavetilt  formulation  as 
given  by  Equation  (3.4)  is  its  restriction  to  radiation  fields.  It 
is  therefore  applicable  to  only  the  radiation  component  of  the  total 
electric  field.  However,  in  Figure  3.5  the  wavetilt  formulation  was 
applied  to  the  complete  vertical  electric  field.  We  now  provide 
justification  for  that  procedure.  In  Figure  3.6  we  show  the 
piecewise-linear  version  of  the  vertical  electric  radiation  field  at 
a  distance  of  3  km  from  the  channel.  This  field  is  obtained  by 
determining  the  complete  field  at  100  km  (which  is  essentially  all 
radiation)  and  extrapolating  it  to  3  km  using  an  inverse  distance 


Figure    3.5.  Piecewise-linear    version    of    the    vertical     electric 

field  given  in  Figure  3.3  and  the  derived  horizontal 

field  for   e^  =  10  and   a  =  10"-^  mhos/m. 
r 
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Figure    3.6.  Piecewise-linear    version    of    the    vertical     electric 

field  calculated  at  100  km  and  extrapolated  to  3  km 
and  the  derived  horizontal  field  for  e  =  10  and  a  = 
10"-^  mhos/m. 
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relationship,  in  accordance  with  standard  practice.  In  Figure  3.6 
the  horizontal  electric  field  derived  from  this  radiation  field  is 
shown  for  e  =  10  and  a  =  10  mhos/m.  A  comparison  between  Figures 
3.5  and  3.6  indicates  that  the  peak  horizontal  field  is  essentially 
the  same  in  both  instances.  This  result  is  not  very  surprising, 
since  from  Equation  (3.12)  it  is  clear  that  the  peak  horizontal  field 
occurs  at  a  time  which  corresponds  approximately  to  the  time  at  which 
the  vertical  electric  field  has  the  maximum  slope.  However,  the 
horizontal  electric  field  shown  in  Figure  3.6  decays  to  zero  at  a 
much  faster  rate  than  that  shown  in  Figure  3.5.  As  discussed 
earlier,  the  wavetilt  formulation  may  also  be  applied  to  the  magnetic 
field  of  the  radiated  wave.  The  resultant  horizontal  electric  field 
derived  from  the  radiation  magnetic  field  would  be  identical  to  that 
shown  in  Figure  3.6.  In  Figure  3.7  we  show  the  results  of  applying 
the  formulation  to  the  complete  magnetic  field.  Again,  the  peak 
horizontal  field  is  the  same.  However,  its  decay  back  to  zero  is  not 
as  fast  as  that  in  Figure  3.6  but  is  faster  than  that  of  the 
horizontal  electric  field  shown  in  Figure  3.5  because  of  the  absence 
of  the  electrostatic  component  in  the  magnetic  field.  Calculations 
of  the  horizontal  field  were  also  performed  using  the  vertical  fields 
at  200  m.  Again,  the  peak  value  of  the  horizontal  electric  field 
derived  from  the  complete  vertical  electric  field  is  essentially  the 
same  as  that  obtained  from  the  radiation  waveform  alone.  However, 
the  error  involved  in  the  horizontal  field  decay  back  to  zero  is 
larger  at  200  m  than  at  3  km  due  to  the  presence  of  a  larger 
electrostatic  component.  To  summarize,  we  conclude,  that  the 
wavetilt  formulation  may  be  applied  to  the  complete  vertical  electric 


Figure  3.7.  Piecewise  linear  version  of  the  magnetic  field  given 
in  Figure  3.3  and  the  derived  horizontal  electric 
field  for  e^  =  10  and  a  =   10"-^  mhos/m. 
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field  or  the  complete  magnetic  field  waveforms  without  much  error  in 
the  initial  peak  although  small  errors  may  occur  later.  However,  the 
only  measurements  of  the  magnetic  field  were  made  on  the  triggered 
system  on  a  time  scale  of  200  us  full-scale  as  discussed  in  Chapter 
II,  and  do  not  provide  sufficient  resolution  for  accurate  analysis, 
and  hence  the  complete  vertical  electric  field  waveform  will  be  used 
in  this  thesis. 

We  have  explained  how  the  vertical  and  horizontal  components  of 
the  electric  field  may  be  derived  in  the  vicinity  of  a  power  line 
using  the  model  of  Master  et  al .  (1981)  with  the  assumption  that  the 
ground  is  a  perfect  conductor.   Effects  of  an  imperfectly  conducting 
ground  plane  may  then  be  incorporated  as  a  perturbation  to  the  field 
solution  derived  for  a  perfectly  conducting  ground.   Either  the 
vertical  electric  field  predicted  by  the  model,  or  the  vertical  field 
directly  measured  near  the  line,  may  be  used  to  compute  the 
horizontal  field  due  to  the  ground  tilt.  In  this  thesis,  we  shall  do 
both:  we  shall  first  present  a  sequence  of  calculated  line  voltages 
for  various  ground  conditions  using  model  electric  fields.  We  shall 
then  use  the  measured  vertical  electric  field  near  the  test-line  to 
compute  a  voltage  waveform  to  be  compared  with  the  simultaneously 
measured  voltage.   For  ground  conductivities  lower  than  10"^  mhos/m, 
the  horizontal  electric  field  computed  from  the  wavetilt  is  much 
larger  than  that  obtained  from  the  model  of  Master  et  al.  (1981) 
for  a  line  height  of  10  m.   For  example,  from  Figures  3.5  and 
3.6,  a=   10"^  mhos/m  leads  to  a  horizontal  field  peak  of  9  V/m 
whereas  from  Figure  3.3  vve  find  an  initial  peak  of  only  0.4 
V/m.  Calculations   indicate   that   for   lightning   over   sea- 
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water  (a=  1  mhos/m;  e^  =  80)  the  peak  horizontal  electric  field 
calculated  from  the  wavetilt  is  of  the  same  order  of  magnitude  as 
that  calculated  from  the  model  of  Master  et  al.  (1981).  For  such 
high  conductivities,  both  horizontal  fields  would  have  to  be  taken 
into  account.  Even  for  10"^  mhos/m,  the  value  of  the  horizontal 
field  after  the  initial  peak  derived  from  the  wavetilt  is  of  the  same 
order  of  magnitude  as  the  value  derived  from  the  model  of  Master  et 
al.  (1981).  However,  in  this  thesis  we  shall  only  use  the  horizontal 
electric  field  derived  from  the  wavetilt  formula  and  the  vertical 
electric  field.  We  now  proceed  to  the  coupling  analysis  to  determine 
how  these  fields  induce  the  line  voltages. 

3.3  Induced  Line  Voltage 

In  this  section  we  discuss  how  the  electromagnetic  field 
generated  by  nearby  lightning  is  coupled  onto  an  overhead  power 
line.  We  shall  first  derive  the  coupled  differential  equations  for 
the  line  voltages  and  currents,  starting  with  Maxwell's  Equations  in 
the  time  domain.  These  differential  equations  will  then  be  converted 
into  a  set  of  difference  equations  which  can  be  solved  on  a  digital 
computer. 

3.3.1  Theoretical  Model 


The  heights  above  ground,  at  the  poles  and  midspans,  for  both 
the  conductors  of  our  unenergized  test-line  have  been  given  in  Table 
2.1.  From  these  values  we  find  average  heights  at  the  pole  and  at 
midspan  for  each  conductor,  and  hence  determine  an  average  value  of 
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the  sag  for  each  conductor.  The  average  height  of  the  line  is  then 
computed  using  the  following  engineering  approximation, 

h  =  hp  -  I  s 
where, 

h     =  effective  height  of  the  line 
hp  =  height  of  the  line  at  the  pole 
s     =  sag 

This    leads    to    the   following    values   which   will    be   used   in   this 
study, 

h^  =  9.8  m  —  height  of  phase  (top)   conductor 

h^  =  7.6  m  —  height  of  neutral    (bottom)   conductor. 

The  phase  and  neutral  conductors  were  of  the  same  type,  with 
diameters  equal  to  8  mm.  We  have  now  completely  defined  the  geometry 
of  the  unenergized  test-line.  We  have  to  determine  the  response  of 
this  line  when  it  is  excited  by  a  transient  electromagnetic  field  due 
to  nearby  lightning.  This  field  is,  in  general,  a  non-uniform 
electromagnetic  field. 

The  solution  to  this  problem  was  first  obtained  in  the  frequency 
domain  by  Taylor  et  al.  (1965),  and  a  more  convenient  form  of  this 
solution  was  presented  by  Smith  (1973).  Whitescarver  (1969)  started 
with  the  solution  given  by  Taylor  et  al.  (1965)  and  modified  it 
according  to  the  formulation  given  by  Putzer  (1968).  He  then  used 
these  results  to  make  frequency-domain  comparisons  between  the  theory 
and  his  experimental  results.  The  frequency-domain  solution  was 
transformed  by  Perala  (1974)  into  a  form  more  suitable  for  computer 
simulation.         This     solution     indicated     that    the     transmission     line 
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problem    could    be     solved     completely     in    terms    of    the    tangential 
electric   field  at   the   conductor   boundaries   and   the  vertical    electric 
field   at  the   end  points  of  the   line,   and  the  vertical    electric  field 
at    the    calculation    point.       A    general    time-domain    solution    to    the 
problem  was  presented  recently  by  Agrawal    et  al .    (1980).     The  various 
assumptions    underlying    the    solution    were    more    closely    examined.      A 
computational    scheme  was  outlined  using  finite-difference  techniques. 
For  the   sake   of  completeness,   we   shall    derive  the  Telegrapher's 
Equations    for   the   two-conductor   line    under    study,    using   a  perfectly 
conducting   ground  as   the   reference   conductor.     The  equations  will    be 
formulated    in    the    time-domain.       The    approach    adopted   will    parallel 
that    used    by    Agrawal    et    al.    (1980).       It    is    important   to    note    that 
most    of   the   earlier    studies    conducted   with    the    specific   purpose    of 
understanding   induced   voltages   on   power   lines,    i.e.    Lundholm    (1-957), 
Rusck    (1958),   Chowdhuri    (1956),    and  Singarajah    (1971),   did  not  follow 
the   approach    indicated   below,    but   assumed   that  the   line   voltage  was 
due  to  the  vertical   electric  field  coupling  alone. 

In  Figure  3.7  we  show  the  two-conductor  test  line  above  a  per- 
fectly conducting  ground.  The  effect  of  finite  ground  conductivity 
will  be  introduced  later.  Let  F(x,y,z,t)  and  ¥(x,y,z,t)  be  the  total 
electric  and  magnetic  field  in  the  vicinity  of  the  line.  However, 
due  to  the  orientation  of  the  axes,  y  =  0  for  all  points  on  the  line 
and  hence  we  omit  the  y  dependence;  the  time  dependence  will  be 
suppressed  for  brevity.  To  derive  the  first  Telegrapher's  Equation 
we  apply  Faraday's  Law  to  the  top  phase  conductor  over  the  region 
shown  by  dotted  lines  in  Figure  3.7.     Since, 
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cjl    .  dl    =  -]^  //F   .-Js,  (3.13) 


/  MEJx  +  Ax,z)   -  E   {x,z)}dz 


0  ^ 


x  +  Ax  x  +  Axh, 

-/  E(x.hi)dx=^     /  /  ^  B^{x.z)dz  dx     (3.14) 

X  ^  "^^     X  0         -^ 


where  we  have  made  use  of  the  boundary  condition 

E^lx.O)  =  0  (3.15) 

Dividing  each  term  in  Equation  (3.14),  by  Ax,  and  taking  the  limit 
as  Ax  -»■  0  we  have 


1^  /  ^  E^(x,z)dz  -  E  (x,hi)  =  -^  /  ^  B  (x,z)dz        (3.16) 
0  0    -^ 


At  this  point  it  is  convenient  to  split  the  complete  field  into  two 
components  (a)  "incident"  field  due  to  lightning,  which  would  be 
present  in  the  absence  of  conductors  and  (b)  "scattered"  field  due  to 
the  current  and  voltage  induced  in  the  conductor  by  the  incident 
field.  Thus, 


^z  =  ^;  ^  ^l 


h=^l     '^l  (3.17) 


I     =  B^  +  B^ 

y   y   y 
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where  the  superscripts  i  and  s  represent  the  incident  and  the 
scattered  components  of  the  total  field.  Substituting  these  into 
Equation  (3.16),  we  obtain 


1^  /  ^  E^  (x.z)dz  -  ^  /  ^  B^  (x,z)dz  -  E^(x,h^) 

=  -|^  /^  E^  (x,z)dz  +-|^  /^Bj(x,z)dz        (3.18) 
0  0    -^ 


This  equation  is  exact.  However,  we  shall  now  make  an  assumption 
that  the  scattered  fields  are  transverse  magnetic.  This  is  a  valid 
assumption  since  the  diameter  of  the  conductor  (~  1  cm)  is  very  small 
compared  to  the  line  height  (-10  m) .  This  ensures  that  the 
following  two  properties  hold  (Stratton,  1941,  pg.  352).  (a)  the 
scattered  voltage  is  single  .valued  in  the  transverse  plane,  and 
defined  by 


s        ^  . 
VJ  (x)  =  -  /  ^  E^  (x,z)dz  (3.19) 

0 


and  (b)  the  scattered  magnetic  flux  can  be  directly  related  to  the 
inductances  per  unit  length  of  the  line. 


/   B^  (x,z)dz  =  4ili(x)  +  42^^^^  ^2-2°^ 


where, 

L..  --  self-inductance  of  the  phase  conductor 
and      L^2  "-  ™tual  inductance  between  the  phase  and  the  neutral 
conductor 
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We  note  that  in  Equation  (3.20)  we  have  not  used  any  superscript  on 
the  currents  because  there  are  only  scattered  currents.  Substituting 
Equations  (3.19)  and  (3.20)  into  Equation  (3.18), 


3V=  (X)  ^  , 


ax 


-^   {4^l^(x)  +  l^^l^M]  +  E^(x,h^) 


^     /^  E^  (x,z)dz  -^  /^BNx,z)dz    (3.21) 
0  0    ■^ 


If  we  apply  Equation  (3.13)  to  the  incident  component  of  the  fields, 
we  obtain 

h  h 

-|^  /  1  E^  (x,z)dz  --^  /  ^  b]  (x,z)dz  =  eJ  (x,h^)   (3.22) 

where  we  have  again  used  the  boundary  condition  given  in  Equation 
(3.15).  The  total  electric  field  parallel  to  the  line,  E  (x,hj  in 
Equation  (3.21),  is  related  to  the  scattered  current,  and  produces  a 
voltage  drop  per  unit  length  along  the  line.  This  relation  may  be 
expressed  in  the  frequency  domain  using  an  impedance  per  unit 
length  of  the  line.  The  internal  inductance  of  the  conductor  is 
generally  negligible,  and  hence  we  shall  assume  that  the  relation 
between  E^(x,h^)  and  I(x)  is  purely  resistive. 

E^(x,h^)  =  Rjjl^(x)  +  Ri2^2^x)  (3.23) 

where  R^^  and  R^^  ai^e  the  self  and  mutual  resistance  per  unit  length 
for  conductor  1.  Substituting  Equations  (3.22)  and  (3.23)  into 
Equation  (3.21),  we  get  the  first  Telegrapher's  Equation  for  our 
line. 
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3V^    (X)        3 


8X 


--at    [LijHlj]-    [Rlj][lj]=E;;(x,h^)  (3.24) 


A  similar  analysis  for  the  neutral   conductor  leads  to 


^^2    (^^         3 


3X 


^It   [4jJ[^jJ  -^    t^2jJ[^j^=  ^x   ^^'^2^  (3.25) 


Equations  (3.24)  and  (3.25)  may  be  combined  together  into  a  matrix 
form, 

^  [V'(x)]+  [R][l(x)]  +  [L]-^  [I(x)]  =  [E^  (x,h)]  (3.26) 

To     derive     the     second     Telegrapher's    Equation    we    start    with 
Maxwell's  formulatiorv  of  Ampere's  Law 


VxH  =  J+i^  (3.27) 


and  integrate  over  a  closed  cylindrical  surface  of  length  Ax  just 
outside  the  phase  conductor.  Using  the  Gauss  Divergence  Theorem, 

(jl^    .^  +  ^  ^TJ^  .^  =  0  (3.28) 

We  have  already  made  the  assumption  that  the  current  in  the  wire 
flows  parallel  to  its  axis,  and  hence 


/  T  •  "33"  =  I  (x  +  ax)  -  L(x)  (3.29] 

ends        ^  ^ 


If  we  assume  that  the  air  conductivity  is  negligible, 
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/  T  •  "?s  =  0  (3.30) 

cyl 


We  substitute  Equations  (3.29)  and  (3.30)  into  Equation  (3.28), 
divide  by  Ax,  and  take  the  limit  as  Ax  ->■  0, 


8Ii(x)   3Q.(x) 

+  -4;r-=  0  (3.31) 


8x     at 


where  Q.  is  the  net  free  charge  per  unit  length  on  the  phase 
conductor,  obtained  from  Gauss's  Law, 

(/U^  ."a^  =  Q^  (3.32) 

Q^  may  be  related  to 'the  capacitances  per  unit  length  of  the  line  and 
the  scattered  voltages. 


\  -  'lA  '  hz'l  (3.33) 


where, 

C^^   --  capacitance  between  phase  conductor  and  ground, 
C^P  ~  capacitance  between  phase  and  neutral   conductors. 

As  shown  in  Agrawal  et  al .  (1980),  Equation  (3.33)  rests  on  the 
assumption  that  the  scattered  axial  electric  field  E^  does  not  vary 
significantly  along  the  x-direction.  Substituting  Equation  (3.33) 
into  Equation   (3.31)  we  have  the  second  Telegrapher's  Equation, 


^         +If   [CiJ[v;(x)]  =  0  (3.34) 


3X  3t    L^lj  JL'j 
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Similarly  for  the  neutral   conductor 


-ir-'lt    [C23][Vj'(x)]=0  (3.35) 


Combining  Equations  (3.34)  and  (3.35)  into  a  matrix  formulation  leads 
to 

4-  tKx)]  +    [C]4   tv'(x)]  =  0  (3.36) 

Equations  (3.26)  and  (3.36)  are  the  coupled  set  of  first-order 
partial  differential  equations  which  may  be  solved  for  the 
variables  [l(x)  ]  and  [v^(x)  ].  Note  that  the  source  term  is  present 
in  only  one  equation.  However,  the  boundary  conditions  at  the  ends 
of  the  line  also  involve  a  knowledge  of  the  vertical  electric  field, 
as  we  shall  see  later.  The  current  [l(x)  ]  represents  the  actual 
current  which  is  flowing  in  the  line.  However  [V^(x)  ]  represents 
only  the  scattered  voltage.  Thus,  instead  of  solving  the 
Telegraphers'  Equations  for  the  currents  and  actual  voltages  on  the 
line,  we  solve  only  for  the  currents  and  the  scattered  voltages.  The 
advantage  of  this  formulation  will  be  discussed  below.  After 
obtaining  the  solution  for  the  scattered  voltage,  the  actual  voltage 
at  any  point  x  on  the  line  may  be  determined  in  accordance  with 
Equations  (3.17)  and  (3.19), 


[V(x)]  =  [V^x)]  +  [v'(x) 


h  . 

=  [V^(x)]  -  [/  E^  (x,z)dz]  (3.37) 


We  now  transform  Equations  (3.26)  and  (3.36)  into  another  form  where 
the  variables  are  the  actual  currents  on  the  line  and  the  total 
voltage.     Let  us  substitute  Equation   (3.37)   into  Equation   (3.26), 


.  h   . 

ir  {[V(x)]+    [fEl   (x,z)dz]}+    [R][l(x) 

0 


+    [L]^  [Kx)]  =    [EJ   (x,h: 


Hence, 


l^  [V(x)]+    [R][l(x)]+    [L]^  [Kx)]  =    [eJ  (x,h)] 

3  ^    ■ 

"1^  [/  E^   (x,z)dz]  (3.38) 

0 


From  Equation  (3.22)  it  is  clear  that  the  forcing  terms  in  Equation 
(3.38)  can  be  rewritten  in  terms  of  the  incident  magnetic  field,  so 
that 


4  [V(x)]+   [R][l(x)]+    [l]^  [IM]=       9     [/V   (x,z)dz]       (3.39) 

0     ^ 


Similarly  Equation   (3.35)  may  be  transformed  into 

4  tl('<H+    [C]^  [V(x)]  —    [C]^  [!\l   (x,z)dz]     (3.40) 

0 

We  see  that  the  forcing  terms  in  Equations  (3.39)  and  (3.40)  are  much 
more  involved  than  the  single  term  in  the  set  of  Equations  (3.26)  and 
(3.35).   Furthermore,  the  forcing  terms  in  Equations  (3.39)  and 
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(3.40)  involve  taking  derivatives  of  the  incident  field,  and  hence, 
in  general,  are  not  efficient  for  computer  simulation  studies. 

In  this  thesis  the  set  of  coupled  differential  equations  given 
by  Equations  (3.26)  and  (3.36)  will  be  used.  The  effect  of  the 
imperfectly  conducting  earth  will  be  introduced  in  the  formulation  by 
using  the  horizontal  electric  field  determined  from  the  wavetilt. 
Furthermore,  a  resistance  per  unit  length  will  be  introduced  to 
account  for  attenuation  due  to  the  ground.  The  calculation  of  the 
horizontal  electric  field  from  the  vertical  component  was  illustrated 
in  Figures  3.5  and  3.5.  We  now  present  the  theory  to  account  for  the 
per  unit  line  resistance. 

For  a  single  conductor  at  height  h  located  above  earth  which  has 
finite  conductivity,  the  impedance  per  unit  length  of  the  ground 
plane  given  by  Vance  (1977)  is 


where 


Y   =  {j%(  a  +  Jwe)  } 

a,  £  -  conductivity  and  permittivity  of  the  ground. 


Simplified  versions  of  Equation  (3.41)  may  be  written  as  follows.  If 
the  earth  is  a  good  conductor  to  the  extent  that  the  displacement 
current  term  is  negligible,  then 


_  1  +  J 
^q  ~  JW^  (3.42) 
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If  the  earth  is  a  relatively  poor  conductor  such  that  the  skin  depth 
is  large  compared  to  the  height  h. 


2g  =-r"+  J'^2^^°9  T= (3.43) 

^  ^  Y^h 


where 


5  =  l/(TTfM^a) 
Yg  =  1.7811 


Consider  a  ground  plane  with  a=  10"  mhos/m  and  e  =  3.   Then,  at 
1  MHz, 

6  =  1/  (tt  X  10^  X  4ti  X  lO'-^x  10"^  J  ^  =  50 


For  h  =  10  m  and  using  Equation  (3.43)  we  obtain 


Zg  =  0.99  +  ja)6.01  x  10'^  (3.44) 


Thus  we  see  that  the  impedance  in  this  case  is  almost  purely 
resistive.  The  resistance  per  meter  of  the  copper  conductor  is 
negligible  compared  to  the  1  Ohm/m  calculated  in  Equation  (3.44).  On 
the  other  hand  the  inductance  per  unit  length  of  a  conductor  10  m 
above  ground  is  about  two  orders  of  magnitude  larger  than  that  due  to 
the  earth  conductivity. 

3.3.2  Computer  Solution 

The  coupled  differential  equations  given  by  Equations  (3.26)  and 
(3.36)  will  now  be  converted  into  a  set  of  difference  equations  using 
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the  point-centered  finite-difference  technique.  We  shall  not  dwell 
on  the  existence  and  uniqueness  of  the  solutions  to  these  equations; 
the  interested  reader  is  referred  to  Lax  (1958).  The  difference 
formulation  of  Equation   (3.26)  centered  about  the  point  (x,t)   is, 

^V(x.Ax/2,t)^V(x-   Ax/2,t)j^    [R][l(x.t)] 

.    [L][^(^»t^^V2)^-   Kx.t-    At/2)  J  __    [,^(^^^^j  (3^45^ 

where  the  source  term  for  conductor  j  is  given  by 


E5j(x,t)  =  E^(x,hj,t)  (3.46) 


Similarly,  Equation   (3.36)  may  be  written  as 

rl(x  +  Ax/2,t)   -   I(x  -  Ax/2.t)  1 
^  AX  J 

^   ^^^(JU,t^  At/2)^-  V(x,t  -   At/2)  J  ^  Q  (3^4^^ 

Note    that   in   Equations    (3.45)    and   (3.47)   we   have   dropped   the    super- 
script on  the  scattered  voltage  V^  for  notational   brevity. 
Let  us  introduce  the  following  notation: 


VJJ  =    [vKk  -  1]ax,    (n  -  l)At 


ij  =   [I{(k  -i)Ax,    [n  -l]At 
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V"     =    [Vj[k  4)Ax,    (n  -^)At}] 


where, 

n  =  1,2,3,    "max"'"^ '    (t"""!^  count) 

k  =  1,2,3,    '^max'*'^ '   (position  count) 

and, 

^^.^^   =  L  =  total  line  length: 
max  ^  ' 

"m,^/'^"t  =  T  i  total  time  of  interest, 
max 

Hence,  Equation  (3.45)  becomes 

wn+1   yn+1      ^n+1  ,  ,n      ,n+l   ,n   v""^^  +  v" 
!w_\_,[,]iL^.[L]^L_lik=i_i     ,3.48, 

and  Equation  (3.47)  gives 


From  Equations  (3.48)  and  (3.49)  we  obtain 


\  Ut      2^     lUt      2J^k  ^  2 AX ^  ^^•^°' 


and 


-1   l"      -   l" 
yn+1   =  v"  +    [Ll         "^-^         "^  (3.51) 

k  k         At-'  M. 


The    solution    scheme   presented    in   Equations    (3.50)    and    (3.51)    may   be 
understood  in  terms  of  Figure  3.8  and  the  following  algorithm: 


o 

• 

( > 

o 

(U 

x: 

4-> 

-tJ 

(0 

3 

^_ 

CT 

o 

(1) 

<u.> 


186 


lV'(l.-u)=uA^ 
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1.  Initialize  the  current  array. 

2.  Define  voltage  at  all  nodes  at  n  =  1  as  zero. 

3.  Compute  l[J  for  k  =  1,  2,  . . . .  k^^^. 

4.  Compute  vj"*"^  for  k  =  2,  3,  ....  k^^^. 

5.  Compute  V^  ,  v|^  ^.   from  boundary  conditions. 

max 

6.  Store  information  for  plots. 

7.  Do  steps  3  through  6  f or  n  =  1 ,  n^^^. 

The  solution  proceeds  sequentially,  within  each  step  of  the  do- 
loop  we  solve  for  one  row  of  current  and  the  next  row  of  voltage,  as 
shown  in  above  diagram. 

A  close  look  at  Equation  (3.51)  indicates  that  it  cannot  be  used 
for  points  k  =  1  and  k  =  k^a^"*"^'  ^ '^-  ^^^  ^"^  points.  That  is  not 
very  surprising  since  these  two  points  define  the  boundaries  of 
the  transmission  line.  Solutions  to  the  voltages  at  these  end  points 
are  determined  by  the  terminations  between  various  conductors  and 
between  the  conductors  and  ground.  With  reference  to  Figure  3.8  we 
note  that  the  current  is  not  known  at  the  boundary,  but  is  known  at 
points  Ax/2  and  3 Ax/2  inside  the  boundary.  To  determine  the  current 
at  the  end  we  assume  that  currents  at  these  two  closest  nodes  and  the 
current  at  the  end  itself  can  be  approximated  by  a  parabolic  curve 
at  each  time  t.  Using  the  coefficients  of  that  fit  the  value 
of  81/ 3x  may  be  obtained  and  then  plugged  directly  into  Equation 
(3.36)  to  determine  the  end  voltage.  We  illustrate  this  procedure  by 
determining  the  voltage  at  x  =  0,  for  some  time  step.  We  have, 


I  =  ax2  +  bx  +  c  (3.52) 

where  a,  b,  c  will  be  determined. 

If  I^  and  I^  define  the  currents   at  points  x  =  ax/2  and 
X  =  3AX/2,  and  if  l"  denotes  current  at  x  =  0,  then 

=  ml"l-'2-8lS}  (3.53) 


x=0 


Let  us  assume  that  the  termination  network  at  x  =  0  is  a 
parallel  R-C  network.  This  more  general  termination  as  opposed  to  a 
purely  resistive  termination  and  will  be  needed  to  simulate  the 
effect  of  the  measuring  capacitor  at  the  north  end  of  the  unenergized 
test  line.  It  is  important  to  note  that  the  scattered  voltage  used 
in  our  formulation  in  Equation  (3.35)  is  not  a  physical  real 
voltage.  The  boundary  condition  at  the  end  of  the  line  will 
therefore  be  in  terms  of  the  physical  total  voltage  denoted  by  vJ, 

T      dvj 
^0  =  -  K^'l-    K^-^  (3.54) 

where  [G^]  and  [C^]  are  the  conductance  and  capacitance  terminating 
matrices  at  x  =  0  and 

v[  =  V^  +  vj  (3.55) 

and  V^  represents    the    voltage   due    to   the    incident   vertical    electric 
field,  and  is  given  by 
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..  h    . 

V}  =    [-   /E^(0.z)dz]  (3.56) 

0 

Substituting     Equations     (3.53),      (3.54),     and     (3.55)     into     Equation 
(3.25)  and  writing  it  in  a  difference  form  we  obtain 


f""^""   =    [c  +^  G    +    ^     r    1     irr       4AX  ^     ^    8      „    1  .,n 
^1  ^^  ^  3AX  \  "^  lA^^tJ     ^LC  -  3^Gl  +^^  CJ  V^ 

_^_  r     n         r).       4  At    r      .    f-  i"""^         i" 


fl  . n+1         .  n 

-wtSHvJ         -V}]}  (3.57) 

A  similar  analysis  may  be  carried  out  to  determine  the  voltage  at 
X  =  L.  In  Figure  3.10  we  show  an  equivalent  circuit  for  a  single- 
conductor  line  above  ground.  Note  that  at  the  ends  of  the  line,  the 
vertical  electric  field  is  needed;  at  all  other  points  only  the 
horizontal  electric  field  is  required  for  a  solution  of  the 
Telegraphers'  Equations  for  the  scattered  voltage  and  current. 

To  summarize,  the  solution  to  the  currents  and  scattered 
voltages  on  the  line  may  be  obtained  from  the  knowledge  of  the 
horizontal  electric  field  along  the  conductors  and  the  vertical 
electric  fields  at  the  two  ends.  However,  to  determine  the  actual 
voltage  physically  present  at  any  point  along  the  line,  knowledge  of 
the  vertical  electric  field  at  that  point  is  also  required  as  in 
Equations  (3.55)  and  (3.56). 


(O    (/)    O)  1 —    o 

Qj     O   -r-     (O     S- 


<T3     =  ^ 


(O  ^ 


-^     »/)  "^   -r-     .-C 


C    r-  C 


r—     'I' 


^   -^   -P  TD 
C   ,— 

O     •!- 

cut- 


<J    CDM-   +-> 


191 


CHAPTER  IV 
RESULTS 

4.1  Introduction 

In  this  chapter  we  present  calculations  based  on  the  theory 
presented  in  Chapter  III.  First,  we  derive  typical  voltages  which 
would  be  measured  on  the  short  test  line  for  various  soil  conditions, 
using  typical  waveforms  for  the  return  stroke  vertical  electric  field 
derived  from  the  model  of  Lin  et  al .  (1980)  and  Master  et  al. 
(1981).  Second,  we  perform  similar  calculations  for  a  5  km  long 
line.  Finally,  we  use  some  of  the  data  presented  in  Chapter  II  on 
simultaneously  measured  vertical  electric  fields  and  induced  voltages 
on  the  460  m  test  line  to  make  comparisons  between  the  theory  and  the 
data  on  a  case-by-case  basis. 

4.2  Short  Line 

In  this  section,  we  calculate  induced  voltage  waveforms  on  a 
short  line  produced  by  a  typical  return  stroke  at  various  locations 
for  a  range  of  earth  parameters  and  line  terminations.  The  overhead 
line  is  chosen  to  be  a  500  m  single-phase  13  kV  distribution  line, 
identical     in     construction     to     the    actual     test    line    described    in 
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Chapter  II.  The  line  is  assumed  to  have  a  North-South  orientation. 
The  bottom  neutral  wire  is  grounded  at  the  North  end  in  a  20  Ohm 
resistance,  and  at  the  South  end  in  a  60  Ohm  resistance.  The  voltage 
is  calculated  at  the  North  end  of  the  phase  wire.  Calculations 
indicate  that  the  presence  of  the  terminated  neutral  wire  reduces  the 
value  of  the  voltage  on  the  phase  wire  by  about  20%,  without  making 
any  significant  changes  in  its  waveshape.  However,  when  the  neutral 
wire  is  left  open-circuited,  there  is  no  reduction  in  the  voltage  on 
the  phase  wire.  To  illustrate  the  range  of  possible  voltage 
waveforms  and  their  sensitivity  to  stroke  location,  earth  parameters, 
and  line  terminations,  calculations  will  be  performed: 

(i)  for  the  case  where  the  phase  wire  is  open-circuited 
(10"  Ohms)  at  both  ends  and  for  the  case  where  both  ends  of  the  phase 
wire  are  terminated  in  500  Ohm  resistances; 

(ii)  for  lightning  at  three  distances  from  the  line:  200  m,  1 
km,  and  5  km.  For  each  of  these  distances,  we  choose  five  locations 
for  the  ground  strike  point  at  various  bearings  from  the  line;  and, 

(iii)  for  four  soil  conditions:  a  =  10'^  mhos/m  with  e  =  15;  a  = 
10"-^  mhos/m  with  e^  =  8;  a=  10"'^  mhos/m  with  e  =  3;  and,  a=  10"^ 
mhos/m  with  e  =3. 

It  was  shown  in  Chapter  III  that  the  induced  voltage  at  either 
end  of  the  line  may  be  determined  from  a  knowledge  of  the  incident 
vertical  electric  field  at  the  ends  and  the  incident  horizontal 
electric  field  along  the  line.  The  first  step,  therefore,  is  to 
determine  the  vertical  electric  fields  at  the  ends  of  the  line.  In 
this  study,  we  shall  use  the  typical  subsequent  return  stroke,  with  a 
18  kA  peak  current,  modeled  by  Lin  et  al .  (1980)  and  Master  et  al . 
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(1981).  Since  the  waveform  of  the  vertical  electric  field  measured 
at  two  close  distances,  typically  under  3  km,  varies  considerably, 
the  vertical  electric  fields  to  be  used  at  the  two  ends  of  the  line 
have  been  computed  directly  from  the  model.  For  each  of  these  two 
vertical  electric  fields,  a  horizontal  electric  field  component  is 
calculated  using  the  wavetilt  formulation  described  in  Chapter  III. 
An  average  horizontal  field  waveshape  is  then  computed.  The  com- 
ponent of  this  average  horizontal  field  along  the  line  is  used  with  a 
1/D  extrapolation  as  the  source  function  along  the  line,  the  calcu- 
lated vertical  electric  fields  being  used  as  the  source  functions  at 
the  two  ends.  This  procedure  will  be  adopted  for  studies  with 
lightning  at  200  m  and  at  1  km  from  the  line.  For  lightning  at  5  km, 
the  source  functions  will  be  determined  by  using  a  1/D  extrapolation 
for  the  calculated  vertical  and  horizontal  electric  fields  at  5  km. 
As  discussed  in  Chapter  III,  it  will  be  assumed  here  that  the 
horizontal  electric  field  produced  by  the  wavetilt  is  the  only 
horizontal  source  driving  the  line. 

For  a  conductivity  of  10"^  mhos/m,  the  ground  is  a  reasonably 
good  conductor,  with  a  skin  depth  at  1  MHz  equal  to  5  m.  Hence,  the 
resistance  per  unit  length  introduced  into  the  transmission  line 
formulation,  as  discussed  in  Chapter  III,  will  be  calculated  from 
Equation  (3.42)  and  is  about  0.2  Ohms/m  at  1  MHz.  For  conductivities 
lower  than  10"^  mhos/m.  Equation  (3.43)  is  used,  giving  a  resistance 
per  unit  length  of  about  1  Ohm/m  at  1  MHz.  From  Equations  (3.42)  and 
(3.43)  it  is  clear  that  the  resistance  per  unit  length  is  frequency 
dependent.  However,  we  shall  use  fixed  values,  determined  at  1  MHz, 
in  the  time-domain  calculations. 
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The  calculations  of  the  induced  voltage  waveforms  for  various 
conditions  are  shown  in  Figures  4.1  through  4.12.  In  each  figure  the 
soil  parameters  a  and  e^  are  indicated  and  so  is  the  distance  of  the 
lightning  from  the  line.  In  all  cases,  voltages  are  calculated  at 
the  North  end  of  the  line  shown  by  a  'V  and  the  ground  strike  point 
is  indicated  by  an  asterisk  (*).  The  waveforms  drawn  with  a  solid 
line  are  for  the  case  in  which  the  line  is  open-circuited  at  both 
ends;  the  waveforms  drawn  with  dashed  lines  are  obtained  when  the 
line  is  terminated  in  500  Ohms  at  both  ends.  The  factor  'x  a'  used 
on  some  of  the  waveforms  indicates  that  the  waveform  shown  is  'a' 
times  the  actual  waveform;  that  is,  the  actual  waveform  is  'a'  times 
smaller  than  that  shown  on  the  figure.  The  distance  in  km  and  the 
bearing  of  the  ground  strike  points  from  the  North  end  of  the  line 
are  as  follows: 

(i)   0.2   km         A:    0.20|0^       ;   B:    0.20190^     ;   C:    0.32|143° 
D:   0.541158°    ;   E:   0.701 180° 

(ii)   1  km  A:    l.OOlO^      ;  B:    1.00|90^    ;  C:    1.03|104° 

D:    1.121117°   ;  E:    1.50|180° 

(iii)   5  km  A:   5.00|0^      ;  B:   5.00|45^    ;  C:   5.01|93° 

D:   5.111136°   ;   E:   5.501 180° 

Several    features    of    the   waveforms    shown    in   Figures   4.1    through   4.12 
are  discussed  below. 

(a)  In  Figures  4.1  through  4.12  time  t  =  0  corresponds  to  the 
time  at  which  the  propagating  lightning  fields  illuminate  the  first 
section  of  the  line.  Thus  some  of  the  waveforms  exhibit  a  time 
delay.       For   example,    in    Figure    4.1.    the    voltages    for    lightning    at 


Figure  4.1.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  l^ne  for  lightning  200  m  from  the 
line  for  a  =  10"^  mhos/m  and  ey  =  15.  The  solid 
lines  represent  the  calculated  waveforms  when  the 
line  is  open-circuited  at  both  ends.  The  dotted 
lines  represent  the  calculated  waveforms  for  the 
line  terminated  in  500  a  at  each  end.  The  factor 
'x  a'  on  any  waveform  indicates  that  the  waveform 
shown  is  'a'  times  the  actual  waveform.  The  ground 
strike  points  are  indicated  by  asterisks  (*). 
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Figure  4.2.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  lijie  for  lightning  200  m  from  the 
line  for  a  =  10"^  mhos/m  and  e^  =  8.  The  solid 
lines  represent  the  calculated  waveforms  when  the 
line  is  open-circuited  at  both  ends.  The  dotted 
lines  represent  the  calculated  waveforms  for  the 
line  terminated  in  500  fi  at  each  end.  The  factor 
'x  a'  on  any  waveform  indicates  that  the  waveform 
shown  is  'a'  times  the  actual  waveform.  The  ground 
strike  points  are  indicated  by  asterisks  (*). 
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Figure  4.3. 


Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  lyie  for  lightning  200  m  from  the 
line  for  a  =  10"^  mhos/m  and  ey  =  3.  The  solid 
lines  represent  the  calculated  waveforms  when  the 
line  is  open-circuited  at  both  ends.  The  dotted 
lines  represent  the  calculated  waveforms  for  the 
line  terminated  in  500  Q  at  each  end.  The  factor 
'x  a'  on  any  waveform  indicates  that  the  waveform 
shown  is  'a'  times  the  actual  waveform.  The  ground 
strike  points  are  indicated  by  asterisks  (*). 
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Figure  4.4.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  lirie  for  lightning  200  m  from  the 
line  for  a  =  10"^  mhos/m  and  ey  =  3.  The  solid 
lines  represent  the  calculated  waveforms  when  the 
line  is  open-circuited  at  both  ends.  The  dotted 
lines  represent  the  calculated  waveforms  for  the 
line  terminated  in  500  ^  at  each  end.  The  factor 
'x  a'  on  any  waveform  indicates  that  the  waveform 
shown  is  'a'  times  the  actual  waveform.  The  ground 
strike  points  are  indicated  by  asterisks  (*). 
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Figure  4.5.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  line  for  lightning  1  km  the  line  for 
a  =  lO'*^  mhos/m  and  ey  =  15.  The  solid  lines 
represent  the  calculated  waveforms  when  the  line  is 
open-circuited  at  both  ends.  The  dotted  lines 
represent  the  calculated  waveforms  for  the  line 
terminated  in  500  fi  at  each  end.  The  factor  'x  a' 
on  any  waveform  indicates  that  the  waveform  shown  is 
'a'  times  the  actual  waveform.  The  ground  strike 
points  are  indicated  by  asterisks  (*). 
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Figure  4.6.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  line  for  lightning  1  km  from  the  line 
for  a  =  10"-^  mhos/m  and  ey  =  8.  The  solid  lines 
represent  the  calculated  waveforms  when  the  line  is 
open-circuited  at  both  ends.  The  dotted  lines 
represent  the  calculated  waveforms  for  the  line 
terminated  in  500  Q  at  each  end.  The  factor  'x  a' 
on  any  waveform  indicates  that  the  waveform  shown  is 
'a'  times  the  actual  waveform.  The  ground  strike 
points  are  indicated  by  asterisks  (*). 


207 


DISTANCE  =  1km 

€    =8 

r 

a=  ^Q-^  mhos/m 


x5 


'oH' 


Wv-v- — ' 


.ll  x5 


W-v:^*->x_:i-:.-- 


500m 


*B 


*C 


1 
5kV 
T 


^10/is(-<- 


xO.5 


.^         xO.5 


0: 


*E 


Figure  4.7.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  line  for  lightning  1  km  from  the  line 
for  a  =  10"^  mhos/m  and  sy  =  3.  The  solid  lines 
represent  the  calculated  waveforms  when  the  line  is 
open-circuited  at  both  ends.  The  dotted  lines 
represent  the  calculated  waveforms  for  the  line 
terminated  in  500  n  at  each  end.  The  factor  'x  a' 
on  any  waveform  indicates  that  the  waveform  shown  is 
'a'  times  the  actual  waveform.  The  ground  strike 
points  are  indicated  by  asterisks  (*). 
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Figure  4.8.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  line  for  lightning  1  km  from  the  line 
for  a  =  10"  mhos/m  and  e^  =  3.  The  solid  lines 
represent  the  calculated  waveforms  when  the  line  is 
open-circuited  at  both  ends.  The  dotted  lines 
represent  the  calculated  waveforms  for  the  line 
terminated  in  500  n  at  each  end.  The  factor  'x  a' 
on  any  waveform  indicates  that  the  waveform  shown  is 
'a'  times  the  actual  waveform.  The  ground  strike 
points  are  indicated  by  asterisks  (*). 
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Figure  4.9.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  line  for  lightning  5  km  from  the  line 
for  a  =  10"^  mhos/m  and  ey  =  15.  The  solid  lines 
represent  the  calculated  waveforms  when  the  line  is 
open-circuited  at  both  ends.  The  dotted  lines 
represent  the  calculated  waveforms  for  the  line 
terminated  in  500  ^  at  each  end.  The  factor  'x  a' 
on  any  waveform  indicates  that  the  waveform  shown  is 
'a'  times  the  actual  waveform.  The  ground  strike 
points  are  indicated  by  asterisks  (*). 
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Figure  4.10.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  line  for  lightning  5  km  from  the  line 
for  a  =  10"-^  mhos/m  and  e^  =  8.  The  solid  lines 
represent  the  calculated  waveforms  when  the  line  is 
open-circuited  at  both  ends.  The  dotted  lines 
represent  the  calculated  waveforms  for  the  line 
terminated  in  500  n  at  each  end.  The  factor  'x  a' 
on  any  waveform  indicates  that  the  waveform  shown  is 
'a'  times  the  actual  waveform.  The  ground  strike 
points  are  indicated  by  asterisks  (*). 


215 


*A 


± 

500  V 

T 
Hio^sh- 


DISTANCE  =5km 

er=8 

(7=  10~^  mhos/m 


•6i\ 


/ 


V   ± 

500  m 

T 


*c 


xO.5 


xO.25 


o: 


K         xO 


.5 


.|\  xO.25 


0: 


*E 


Figure  4.11.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  line  for  lightning  5  km  from  the  line 
for  a  =  10"^  mhos/m  and  e^  =  3.  The  solid  lines 
represent  the  calculated  waveforms  when  the  line  is 
open-circuited  at  both  ends.  The  dotted  lines 
represent  the  calculated  waveforms  for  the  line 
terminated  in  500  n  at  each  end.  The  factor  'x  a' 
on  any  waveform  indicates  that  the  waveform  shown  is 
'a'  times  the  actual  waveform.  The  ground  strike 
points  are  indicated  by  asterisks  {*). 
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Figure  4.12.  Calculated  voltage  waveforms  at  the  North  end  of  the 
500  m  overhead  line  for  lightning  5  km  from  the  line 
for  a  =  lO"'^  mhos/m  and  e^  =  3.  The  solid  lines 
represent  the  calculated  waveforms  when  the  line  is 
open-circuited  at  both  ends.  The  dotted  lines 
represent  the  calculated  waveforms  for  the  line 
terminated  in  500  n  at  each  end.  The  factor  'x  a' 
on  any  waveform  indicates  that  the  waveform  shown  is 
'a'  times  the  actual  waveform.  The  ground  strike 
points  are  indicated  by  asterisks  (*). 
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point  A  commence  at  t  =  0;  for  lightning  at  point  E,  they  commence 
after  a  time  lag  of  1.7  ms  which  corresponds  to  the  travel  time  down 
the  line. 

(b)  When  the  ground  strike  point  is  equidistant  from  the  two 
ends  of  the  line  (i.e.,  at  point  C),  the  voltage  waveform  calculated 
for  the  open-circuited  case  is  similar  to  the  vertical  electric  field 
waveform.  The  waveforms  for  the  terminated  case,  however,  look  more 
like  the  derivative  of  the  vertical  electric  field.  This  behavior  is 
exhibited  for  all  distances  and  soil  conditions,  and  is  due  to  the 
fact  that  for  point  C,  the  voltage  coupling  is  essentially  due  to  the 
vertical  electric  field,  and  the  line  behaves  like  a  classical 
vertical  electric  field  antenna.  However,  there  is  some  initial 
ripple  in  the  open-circuit  waveforms  due  to  propagation  and 
reflections  along  the  finite  line  length.  The  ripple  is  stronger  for 
lower  conductivities,  indicating  that  it  is  due  to  the  horizontal 
electric  field  coupling. 

(c)  When  the  ground  conductivity  is  high  {a=  10"^  mhos/m) ,  the 
voltage  waveforms  for  the  open-circuited  case  are  similar  for  all 
ground  strike  points  and  look  like  the  vertical  electric  field. 
However,  the  initial  portion  of  the  waveforms  due  to  lightning  at 
points  A  and  E  (i.e.  parallel  to  the  line)  show  significant 
oscillations  due  to  the  reflections  produced  by  the  coupling  of  the 
horizontal  electric  field.  The  oscillations  appear  to  be  very  small 
at  200  m  range  because  they  are  imbedded  within  the  large  electro- 
static component  of  the  vertical  electric  field.  At  1  km  and  5  km 
range,  when  the  electrostatic  field  is  weaker,  these  oscillations 
become  more  apparent. 
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(d)  When  soil  conductivity  is  reduced  to  from  10"^  mhos/m 
to  10'  mhos/m,  the  horizontal  field  coupling  becomes  comparable  to 
the  vertical  field  coupling  and  the  voltage  waveforms  for  the  open- 
circuit  case  for  lightning  at  point  A  become  bipolar.  As  the 
conductivity  is  reduced  further,  to  a  =  10"^  mhos/m  and  a=  10"^ 
mhos/m,  the  horizontal  field  coupling  becomes  dominant  and  the 
voltage  waveforms  exhibit  negative  polarity  for  lightning  at  points  A 
or  B  and  positive  polarity  for  lightning  at  points  D  or  E.  For 
lightning  at  point  C,  the  induced  voltage  waveforms  were  discussed  in 
(a). 

(e)  We  note  that  in  general  the  induced  voltage  on  the  line  is  a 
function  of  the  vertical  electric  field  coupling  and  the  horizontal 
electric  field  coupling.  The  vertical  electric  field  coupling  tends 
to  produce  a  voltage  of  positive  polarity  irrespective  of  the  ground 
strike  point,  in  conformity  with  the  concept  of  a  positive  induced 
classical  surge.  However,  the  voltage  produced  by  the  horizontal 
electric  field  coupling  is  of  negative  polarity  when  the  ground 
strike  point  is  located  such  that  the  incident  lightning  electro- 
magnetic field  is  propagating  along  the  line  and  away  from  the 
calculation  point  (e.g.,  points  A  or  B).  For  ground  strike  points 
located  such  that  the  incident  lightning  electromagnetic  field  is 
propagating  along  the  line  and  toward  the  point  at  which  the 
calculation  is  done,  the  horizontal  electric  field  coupling  produces 
a  positive  surge  (e.g.,  points  D  or  E).  In  addition  to  the  location 
of  the  ground  strike  point,  the  horizontal  coupling  is  determined  by 
the  soil  parameters,  the  coupling  being  stronger  for  soils  with  lower 
conductivity.  Thus,  it  can  be  seen  that  a  negative  induced  surge  may 
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be  produced  at  one  end  of  the  line  for  soils  for  which  the  horizontal 
coupling  is  the  dominant  voltage  producing  mechanism. 

(f)  The  magnitude  of  the  peak  induced  voltage  is  also  a  function 
of  the  soil  parameters,  the  peak  amplitude  being  larger  for  poorer 
soils.  For  example,  at  a  distance  of  1  km  and  the  ground  strike 
point  as  E,  the  peak  induced  voltage  for  the  open-circuit  case  is 
about  5  kV  for  a  =  10"^  mhos/m  and  about  45  kV  for  a  =  10"^  mhos/m. 
For  the  ground  strike  point  at  A  the  peak  voltage  is  about  5  kV 
for  a=  10"^  mhos/m,  and  about  40  kV  at  a  =  10"^  mhos/m,  in  addition 
to  reversing  its  polarity.  This  indicates  that  in  areas  where  the 
soil  has  a  low  conductivity,  the  induced  line  voltages  would  be 
larger  in  peak  amplitude.  This  is  of  special  significance  to  the 
power  industry  because  under  poor  soil  conditions  grounds  on 
arrestors  and  transformers  often  have  relatively  high  grounding 
resistances.  Higher  equipment  failure  has  often  been  attributed  to 
these  higher  grounding  resistances.  However,  calculations  presented 
here  indicate  that  poor  soil  conductivities  also  lead  to  higher 
induced  voltage  amplitudes  via  horizontal  field  coupling,  so  that 
both  effects  may  be  responsible  for  arrestor  and  transformer 
failures. 

Making  comparisons  between  the  calculations  presented  in  this 
section,  as  shown  in  Figures  4.1  through  4.12,  and  the  data  on 
positive  and  negative  induced  voltages,  as  presented  in  Figure  2.25 
of  Chapter  II,  we  conclude  that  the  effective  soil  conductivity  under 
which  the  data  were  recorded  is  in  the  range  10""^  to  10"^  mhos/m.  If 
the  soil  conductivity  were  high  ( a  ~  10"^  mhos/m)  the  angular 
variation  in  the  voltage  polarity  illustrated  in  Figures  2.26  and  2.7 
would  not  have  been  observed  on  our  450  m  test  line. 
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4.3  Long  Line 

In  this  section,  we  present  calculations  of  induced  voltage 
waveforms  on  a  5  km  long  overhead  line.  The  vertical  electric  field 
illuminating  the  line  is  obtained  from  the  return  stroke  model  of  Lin 
et  al.  (1980)  and  Master  et  al.  (1981)  as  discussed  in  the  previous 
section,  and  the  horizontal  electric  field  is  calculated  from  the 
wavetilt  formulation.  The  overhead  line  is  a  single  conductor 
identical  to  the  short  line  discussed  in  the  previous  section,  except 
for  the  missing  neutral  conductor.  Again,  the  line  is  assumed  to 
have  a  North-South  orientation,  with  the  North  end  short-circuited 
(1  Ohm),  and  the  South  end  open-circuited  (10^  Ohm).  Voltage 
waveforms  are  calculated  at  two  points  along  the  line:  one  at  a 
distance  of  2.2  km  from  the  short-circuited  North  end,  and  the  other 
at  the  open  circuited  South  end.  This  seemingly  peculiar  configur- 
ation is  chosen  because  it  is  representative  of  the  9.9  km  South 
African  test  line  (Eriksson  and  Meal,  1980)  discussed  previously. 
Calculated  waveforms  will  be  shown  for  two  soil  conditions:  a=  10"^ 
mhos/m  with  e  =  15;  and,  a=  10"'^  mhos/m  with  e  =  3.  The  ground 
strike  point  will  be  assumed  to  be  200  m  or  1  km  from  the  line. 
Three  ground  strike  points  are  selected  at  three  different  bearings 
in  each  instance.  The  computation  strategy  is  identical  to  that  used 
for  the  short  line,  except  for  the  fact  that  in  this  case,  we  also 
need  the  vertical  electric  field  at  the  calculation  point  in  the 
middle  of  the  line. 

Figures  4.13  through  4.16  show  calculated  voltage  waveforms  for 
various  cases.   In  each  figure,  the  soil  parameters  a  and  e  are 


Figure  4.13.  Calculated  voltage  waveforms  for  a  5  km  line  for 
lightning  200  m  from  the  line  for  a  =  10"'^  mhos/m 
and  ey  =  15.  The  solid  lines  represent  the  voltage 
waveform  calculated  for  point  V^;  the  dotted  lines 
represent  the  voltage  waveform  calculated  for  point 
V2.  Waveforms  marked  'x  a'  indicate  that  the 
waveform  shown  is  'a'  times  the  actual  waveform. 
The  location  of  the  ground  strike  points  are 
indicated  by  asterisks  (*). 
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Figure  4.14.  Calculated  voltage  waveforms  for  a  5  km  line  for 
lightning  200  m  from  the  line  for  a  =  10"'^  mhos/m 
and  Ej,  =  3.  The  solid  lines  represent  the  voltage 
waveform  calculated  for  point  V-j^;  the  dotted  lines 
represent  the  voltage  waveform  calculated  for  point 
Vp.  Waveforms  marked  'x  a'  indicate  that  the 
waveform  shown  is  'a'  times  the  actual  waveform. 
The  location  of  the  ground  strike  points  are 
indicated  by  asterisks  (*). 
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Figure  4.15.  Calculated  voltage  waveforms  for  a  5  km  line  for 
lightning  1  km  from  the  line  for  a=  10"^  mhos/m  and 
e  =  15.  The  solid  lines  represent  the  voltage 
waveform  calculated  for  point  V,;  the  dotted  lines 
represent  the  voltage  waveform  calculated  for  point 
V2.  Waveforms  marked  'x  a'  indicate  that  the 
waveform  shown  is  'a'  times  the  actual  waveform. 
The  location  of  the  ground  strike  points  are 
indicated  by  asterisks  (*). 
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Figure  4.16.  Calculated  voltage  waveforms  for  a  5  km  line  for 
lightning  1  km  from  the  line  for  a=  10"  mhos/m  and 
ey  =  13.  The  solid  lines  represent  the  voltage 
waveform  calculated  for  point  V^;  the  dotted  lines 
represent  the  voltage  waveform  calculated  for  point 
V2.  Waveforms  marked  'x  a'  indicate  that  the 
waveform  shown  is  'a'  times  the  actual  waveform. 
The  location  of  the  ground  strike  points  are 
indicated  by  asterisks  (*). 
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shown,  the  voltage  points  are  indicated  'Vi'  and  ^M^\  and  the  ground 
strike  points  are  shown  by  asterisks  {*).  The  wavefor.ns  shown  with 
bold  lines  and  labeled  '1'  represent  the  induced  voltages  calculated 
at  point  V^;  the  waveforms  shown  with  dashed  lines  and  labeled  '2' 
are  calculated  at  point  V2.  The  factor  'x  a'  used  with  some  of  the 
waveforms  indicates  that  the  waveform  shown  is  'a'  times  the  actual 
waveform.  The  distance  in  km  and  the  bearing  of  the  ground  strike 
points  from  the  short-circuited  North  end  are: 

(i)     0.2  km       A:    1.02|168.7°   ;  B:   4.00|177.1°   ;  C:   5.2|180° 
(I'i)     1  km  A:    1.4112351      ;   B:   4.1212651      ;  C:   6.0|180° 

We  now  note  some  of  the  interesting  features  of  the  voltage 
waveforms  shown  in  Figures  4.13  through  4.16. 

(a)  Time  t  =  0  in  Figures  4.13  through  4.16  correspond  to  the 
time  at  which  the  incident  electromagnetic  field  due  to  lightning 
illuminates  the  closest  section  of  the  line.  For  example,  for 
lightning  at  point  C.  there  is  no  time  delay  associated  with  the 
calculation  at  point  V^.  Point  V^  being  2.8  km  away  corresponds  to  a 
time  delay  of  about  9  seconds. 

(b)  For  the  good  ground  [0=  lO'^  mhos/m,  e^  =  15),  the  reflec- 
tions in  the  voltage  waveforms  from  the  short-circuited  and  open- 
circuited  ends  of  the  line  may  be  readily  identified.  Recall  that  a 
voltage  pulse  reflected  from  a  short-circuited  end  reverses  polarity, 
whereas  it  retains  its  polarity  when  reflected  from  an  open-circuited 
end.  For  low  ground  conductivity  (a=  10""^  mhos/m,  e^  =  3),  these 
reflections  are  essentially  damped  out  because  of  the  resistance  per 
unit  length  introduced  into  the  transmission  line  formulation. 
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(c)  An  voltage  waveforms  calculated  at  V^^  for  different  soil 
conductivities  and  different  ground  strike  points  are  of  positive 
polarity.  This  is  because  V-|^  being  in  the  middle  of  the  line,  the 
lightning  electromagnetic  field  is  always  travelling  along  the  line 
towards  V^  as  discussed  under  (d)  in  the  previous  section. 

(d)  For  the  low  conductivity  case  {a  =  10"^  mhos/m,  e  =  3)  the 
calculated  voltage  at  V2,  for  lightning  at  point  C,  is  of  negative 
polarity.  This  is  in  keeping  with  arguments  presented  in  (c)  above, 
and  (d)  in  the  previous  section. 

(e)  The  peak  amplitude  of  the  voltage  surges  is  greater  for  poor 
soil  conditions,  because  of  the  fact  that  there  is  more  horizontal 
electric  field  coupling.  The  problems  associated  with  this  larger 
induced  surge  were  discussed  under  (e)  in  the  previous  section. 

Thirty-two  good  quality  recordings  reported  by  Eriksson  et  al . 
(1982)  on  the  South  African  9.9  km  test  line,  as  discussed 
previously,  give  an  average  induced  voltage  peak  of  45  kY  with  a 
standard  deviation  of  35  kV.  The  configuration  and  termination  of 
their  line  was  similar  to  the  5  km  line  discussed  in  this  section. 
Furthermore,  the  location  of  the  calculation  point  V,,  with  respect 
to  the  line  length,  as  shown  in  Figures  4.13  through  4.16,  is  also 
close  to  the  point  at  which  measurements  reported  by  Eriksson  et  al. 
(1982)  were  made.  The  locations  of  most  of  the  ground  strike  points 
for  their  reported  data  (Figure  5  of  Eriksson  et  al . ,  1982)  may  be 
represented  by  locations  A  and  B  in  Figures  4.13  through  4.16. 
Though  no  direct  comparisons  can  be  made  between  the  theory  presented 
here  and  the  South  African  data  due  to  lack  of  electric  or  magnetic 
field  measurements,  there  appears  to  be  a  general  agreement  between 


234 


the  calculated  peak  voltage  shown  in  Figures  4.13  through  4.16  and 
the  data  reported  by  Eriksson  et  al.  (1932)  for  an  effective  ground 
conductivity  of  the  order  of  lO'^  mhos/m.  No  conductivity  measure- 
ments have  been  reported  by  Eriksson  et  al.  (1982);  however  a  value 
of  10"  mhos/m  is  not  unreasonable. 

Similar  calculations  were  also  performed  for  the  2.1  km  French 
telecommunication  line  reported  by  Hamelin  et  al.  (1979)  which  was 
terminated  in  500  Ohms  at  each  end.  Again,  we  find  reasonably  good 
agreement  between  measured  and  calculated  peak  voltages  for  an 
effective  ground  conductivity  of  the  order  of  lO'^  mhos/m,  a  value 
which  has  been  used  by  Djebari  et  al.  (1931)  as  representative  of 
their  "wet  ground"  conditions. 

4.4  Comparison  of  Measured  Voltages  with  Theory 

Pertinent  data  for  simultaneously  measured  vertical  electric 
field  and  induced  test-line  voltages  due  to  nearby  lightning  are 
listed  in  Tables  2.3(a),  (b),  and  (c).  In  this  section,  we  use  some 
of  that  data  to  make  direct  comparisons  between  the  voltage 
measurements  reported  in  this  thesis  and  calculated  voltage 
waveforms.  We  select  six  return  strokes  at  various  distances  and 
angles  from  the  North  end  of  the  line  where  the  voltage  measurement 
was  made.  The  simultaneously  measured  vertical  electric  field  is 
digitized  and  the  horizontal  electric  field  is  determined  from  the 
wavetilt  formulation.  From  the  discussion  in  Section  4.1  it  is  clear 
that  the  effective  ground  conductivity  in  the  vicinity  of  the  test 
line  must  have  been  relatively  low.    Since  no  conductivity 
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measurements  were  made  in  the  area  at  the  time  the  data  were 
recorded,  we  assume  soil  parameters  to  be  a  =  10"^  mhos/m  and  e  =3 
for  the  calculations.   Though  these  values  represent  an  unusually 
poor  soil,  the  experience  of  the  field  engineers  and  line  crews  in 
the  Tampa  Bay  area  indicates  that  ground  conditions  are  poor.   The 
experience  of  the  field  crew  also  indicate  vast  differences  in  soil 
parameters  for  regions  separated  by  a  mere  100  m.  During  the  summer 
of  1982,  three  years  after  the  experiment,  indirect  conductivity 
measurements  were  made  in  the  area  in  which  the  test-line  voltage 
measurements  were  made.  Measurements  of  the  steady-state  resistance 
to  ground  of  a  1.25  cm  diameter  vertically  driven  ground  rod  were 
made  at  1.5  m  depth  intervals  between  1.5  m    and  6.0  m.   Using  a 
formula  given  by  Sunde  (1958)  we  translated  these  ground  rod 
resistance  values  to  equivalent  ground  conductivities.  Measurements 
indicate  that  the  ground  conductivity  decreases  with  depth,  the 
lowest  value  being  about  10"-^  mhos/m  at  a  depth  of  about  6  m.  This 
value  is  in  agreement  with  the  average  effective  conductivity  derived 
from  radio  frequency  propagation  of  lightning  signals  over  200  km  of 
Florida  soil  by  Serhan  et  al.  (1980).  However,  a  value  between  10"'* 
and  10"^  mhos/m  has  to  be  used  to  reproduce  the  positive  and  negative 
polarity  measured  voltages.   Since  there  are  large  local  variations 
in  the  ground  conductivity  and  the  effective  ground  conductivity 
increases  with  depth  in  the  area,  as  discussed  above,  these  low 
values  of  conductivity  may  not  be  totally  unreasonable. 

The  following  strokes  have  been  selected  for  this  study  and 
represent  a  selection  of  ground  strike  points  all  around  the  line. 
The  range  and  bearing  of  the  ground  strike  point  from  the  North  end 
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of  the  line,  at  which  the  voltage  measurement  was  made,  and  the 
stroke  identification  are  given. 

(a)  Range  =  10.7  km;  angle  =  8°.   Stroke  #1  for  the  170455  UT 
flash,  during  run  #79199TR31. 

(b)  Range  =  13.6  km;  angle  =  53°.   Stroke  #1  for  the  171407  UT 
flash,  during  run  #79199TR31. 

(c)  Range  =  1.2  km;  angle  =  82°.  Stroke  #1  of  the  single  stroke 
flash  at  184600  UT,  during  run  #79199TR31. 

(d)  Range  =  5.5  km;  angle  =  92°.   Stroke  #2  for  the  190004  UT 
flash,  during  run  #79196TR27. 

(e)  Range  =  5.8  km;  angle  =  104°.   Stroke  #1  for  the  185842  UT 
flash,  during  run  #79196TR27. 

(f)  Range  =  10.2  km;  angle  =  175°.  Stroke  #2  for  the  220319  UT 
flash,  during  run  #79208TR44. 

In  Figures  4.17  through  4.22  we  present  direct  comparisons 
between  the  data  and  theory.  In  each  case,  the  vertical  electric 
field  obtained  on  40  us  full-scale  from  the  m  channel  of  the 
Instrumentation  Tape  Recorder  is  digitized.  The  wavetilt  computer 
program  listed  in  Appendix  C  is  used  to  determine  the  horizontal 
electric  field  with  the  soil  parameters  assumed  to  be  a=  10"^  mhos/m 
and  e^  =  3.  These  incident  fields  are  then  used  in  the  coupling 
program  (also  listed  in  Appendix  C)  to  determine  the  calculated 
voltages  under  various  terminal  conditions. 

In  Figures  4.17  through  4.22  we  first  show  a  piecewise-linear 
version  of  the  measured  vertical  electric  field,  and  the  calculated 
horizontal  electric  field.  The  next  waveform  is  that  of  the 
simultaneously  measured  voltage  at  the  Morth  end  of  the  line  obtained 


Figure  4.17.  Comparisons  between  data  and  theory  for  the  first 
stroke  in  the  170455  UT  flash  during  run  #79199TR31 
located  at  a  range  of  10.7  km  and  bearing  of  8°  from 
the  North  end  of  the  line  at  which  the  voltage  is 
measured.  The  top  two  traces  are  the  measured 
vertical  electric  field  and  the  calculated 
horizontal  electric  field  for  a  =  10"^  mhos/m  and 
ty  =  3.  The  remaining  traces  show  the  measured 
voltage  waveform  (V^^gaj)  and  the  calculated  voltages 
when  the  South  end  is  open-circuited  (Vqq),  short 
circuited  to  the  60  Ohm  ground  rod  CVsq),  and 
terminated  in  500  Ohms  (V|^,|). 
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Figure  4.18.  Comparisons  between  data  and  theory  for  the  first 
stroke  in  the  171407  UT  flash  during  run  #79199TR31 
located  at  a  range  of  13.6  km  and  bearing  of  53° 
from  the  North  end  of  the  line  at  which  the  voltage 
is  measured.  The  top  two  traces  are  the  measured 
vertical  electric  field  and  the  calculated 
horizontal  electric  field  for  o  =  10"  mhos/m  and 
ey  =  3.  The  remaining  traces  show  the  measured 
voltage  waveform  (Vj^g-jg)  and  the  calculated  voltages 
when  the  South  end  is  open-circuited  (Vqq),  short 
circuited  to  the  60  Ohm  ground  rod  (v^q),  and 
terminated  in  500  Ohms  (Vf^). 
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Figure  4.19.  Comparisons  between  data  and  theory  for  the  first 
stroke  in  the  184600  UT  flash  during  run  #79199TR31 
located  at  a  range  of  1.2  km  and  bearing  of  82°  from 
the  North  end  of  the  line  at  which  the  voltage  is 
measured.  The  top  two  traces  are  the  measured 
vertical  electric  field  and  the  calculated 
horizontal  electric  field  for  a  =  10"  mhos/m  and 
e^,  =  3.  The  remaining  traces  show  the  measured 
voltage  waveform  (V^ggg)  and  the  calculated  voltages 
when  the  South  end  is  open-circuited  (Vqq),  short 
circuited  to  the  60  Ohm  ground  rod  IV5Q),  and 
terminated  in  500  Ohms  (V[,;|). 
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Figure  4.20.  Comparisons  between  data  and  theory  for  the  second 
stroke  in  the  190004  UT  flash  during  run  #79196TR27 
located  at  a  range  of  5.5  km  and  bearing  of  92°  from 
the  North  end  of  the  line  at  which  the  voltage  is 
measured.  The  top  two  traces  are  the  measured 
vertical  electric  field  and  the  calculated 
horizontal  electric  field  for  a  =  10"^  mhos/m  and 
Ey,  =  3.  The  remaining  traces  show  the  measured 
voltage  waveform  (Vj^gaj)  and  the  calculated  voltages 
when  the  South  end  is  open-circuited  (Vqq),  short 
circuited  to  the  60  Ohm  ground  rod  {V5Q),  and 
terminated  in  500  Ohms  (Vf^). 
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Figure  4.21.  Comparisons  between  data  and  theory  for  the  first 
stroke  in  the  185842  UT  flash  during  run  #79196TR27 
located  at  a  range  of  5.8  km  and  bearing  of  104° 
from  the  North  end  of  the  line  at  which  the  voltage 
is  measured.  The  top  two  traces  are  the  measured 
vertical  electric  field  and  the  calculated 
horizontal  electric  field  for  a  =  10"^  mhos/m  and 
Ey,  =  3.  The  remaining  traces  show  the  measured 
voltage  waveform  (V^g^g)  and  the  calculated  voltages 
when  the  South  end  is  open-circuited  (Vqq),  short 
circuited  to  the  60  Ohm  ground  rod  (Vsq),  and 
terminated  in  500  Ohms  (V[^). 
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Figure  4.22.  Comparisons  between  data  and  theory  for  the  second 
stroke  in  the  220319  UT  flash  during  run  #79203TR44 
located  at  a  range  of  10.2  km  and  bearing  of  175° 
from  the  North  end  of  the  line  at  which  the  voltage 
is  measured.  The  top  two  traces  are  the  measured 
vertical  electric  field  and  the  calculated 
horizontal  electric  field  for  o  =  10"  mhos/m  and 
e^  =  3.  The  remaining  traces  show  the  measured 
voltage  waveform  (V^^ggj)  and  the  calculated  voltages 
when  the  South  end  is  open-circuited  (Vqq).  short 
circuited  to  the  60  Ohm  ground  rod  (V^q),  and 
terminated  in  500  Ohms  (Vj^). 
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from  the  electric  field  measurement  inside  the  specially  designed 
capacitor.  This  experimentally  obtained  waveform  is  compared  with 
three  calculated  waveforms.  The  first  calculated  waveform  labeled 
Vqq  is  obtained  when  the  South  end  is  left  open-circuited  (10^ 
Ohms).  From  all  available  records  of  the  experiment,  it  appears  that 
this  was  the  condition  under  which  all  measurements  reported  in  this 
thesis  were  made.  Calculated  waveforms  are  also  presented  for  the 
case  in  which  the  South  end  is  clamped  to  the  grounded  60  Ohm 
downlead  (V^q)  and  for  the  case  when  the  South  end  is  terminated  in  a 
500  Ohm  resistance  (Vf^). 

In  Figures  4.17  through  4.22,  with  the  exception  of  Figure  4.19, 
the  calculated  waveform  for  the  open-circuited  case  looks 
consistently  reasonable  compared  to  the  measured  waveform.  However, 
there  is  a  consistent  discrepancy  in  peak  amplitude:  the  calculated 
values  are  lower  than  the  measured  values  by  about  a  factor  of  4. 
The  waveforms  for  the  matched  case  (South  end  terminated  in  500  Ohms) 
also  look  very  similar  to  the  measured  waveforms,  except  for  the 
amplitude  of  the  calculated  waveform  which  is  again  lower  than  the 
measured  amplitude  by  about  a  factor  of  4.  The  waveforms  calculated 
for  the  short-circuited  case  (South  end  terminated  in  60  Ohms)  still 
have  features  similar  to  the  measured  waveforms,  however  the 
agreement  between  theory  and  data  is  not  as  good  as  that  obtained  for 
the  open-circuit  or  matched  case,  especially  in  Figures  4.17  and 
4.18.  In  Figure  4.19  the  agreement  between  theory  and  data  is  not 
very  good  for  all  three  cases  of  the  calculated  waveform.  The 
measured  waveform  is  bipolar  in  nature  with  an  initial  positive 
excursion  followed  by  an  approximately  equal  negative  excursion.  The 


250 


calculated  waveforms  exhibit  a  very  small  positive  excursion  for  the 
first  2  ys,  after  which  the  voltage  is  essentially  negative.  The 
amplitude  of  the  negative  peak  in  the  calculated  waveform  for  the 
open-circuited  case  is  again  about  4  times  smaller  than  the  negative 
peak  in  the  data.  The  differences  between  the  theory  and  data  in 
this  case  may  be  attributed  to  return  stroke  channel  tilt  or 
tortuosity. 

The  consistent  disagreement  in  the  voltage  amplitudes  between 
the  measured  and  calculated  waveforms  may  be  due  to  two  principal 
reasons.  It  is  conceivable  that  the  calibration  for  the  voltage  is 
in  error,  especially  the  calibration  constant  associated  with  the 
measuring  capacitor.  This  possibility  was  raised  in  Chapter  II  and 
in  Appendix  D.  However,  after  a  thorough  check  of  the  available 
records  of  the  calibration,  we  can  find  no  such  error.  Another 
possible  source  of  error  could  be  in  the  theory.  In  Chapter  III  the 
wavetilt  formulation  was  used  to  compute  the  horizontal  electric 
field  from  the  vertical  component.  The  formulation  was  applied  to 
the  complete  vertical  electric  field  instead  of  applying  it  to  only 
the  radiation  component  of  the  vertical  electric  field.  This 
procedure  was  justified  by  the  fact  that  essentially  the  same  initial 
horizontal  electric  field  would  be  obtained  in  either  case.  This, 
however,  leads  to  the  conclusion  that  the  induction  and  electrostatic 
components  in  the  vertical  electric  field  do  not  make  significant 
contributions  to  the  horizontal  electric  field.  It  is  possible  that 
the  induction  and  electrostatic  components  in  the  vertical  electric 
field  do  give  rise  to  a  horizontal  component  which  is  determined  in  a 
manner  different  from  the  wavetilt  formulation  outlined  in  Chapter 
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III.  If  that  were  the  case,  the  additional  horizontal  electric  field 
would  produce  a  greater  coupling  and  hence  lead  to  higher  voltage 
amplitudes. 


CHAPTER  V 
CONCLUSIOMS 

In  this  thesis  we  have  presented  data  on  the  induced  voltages 
measured  on  a  specially  constructed  test-line  and  simultaneously 
recorded  vertical  electric  fields  due  to  return  strokes,  correlated 
with  the  location  of  the  lightning  ground  strike  point.  For  the  460 
m  test-line,  the  measured  voltage  waveforms  exhibit  a  wide  range  of 
variability,  depending  on  the  location  of  the  lightning  ground  strike 
point  around  the  line.  In  particular,  we  presented  data  to  show  that 
when  the  ground  strike  point  is  along  the  imaginary  line  representing 
an  extension  of  the  test-line,  the  induced  voltage  is  maximum. 
Furthermore,  we  indicated  that  similar  lightning  return  strokes  along 
this  imaginary  line  could  produce  voltages  of  either  polarity 
depending  on  the  location  of  the  ground  strike  point.  For  lightning 
along  the  perpendicular  bisector  to  the  line,  the  induced  voltages 
were  found  to  be  relatively  small  in  peak  amplitude. 

These  data  appear  to  be  in  conflict  with  most  of  the  previously 
measured  induced  overvoltages  which  are  of  positive  polarity,  the  so- 
called  "classical  induced"  surge.  Classical  theories  used  to  explain 
the  induced  surges  on  power  lines  couple  only  the  vertical  electric 
field  due  to  lightning  and  produce  only  positive  surges.  It  is 
obvious  that  the  mechanism  of  induced  voltage,  at  least  on  our  test- 
line,  is  different  because  of  the  fact  that  measured  induced  surges 
are  of  either  polarity. 
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We  have  therefore  presented  a  new  theory  which  can  explain  both 
positive  and  negative  surges.  Some  of  the  elements  of  the  theory 
have  been  developed  over  the  years  to  determine  coupling  on  wires 
produced  by  a  nuclear  electromagnetic  pulse  (NEMP)  caused  by  a  high- 
altitude  nuclear  explosion.  However,  the  theory  has  never  been 
applied  to  specifically  study  lightning  induced  power-line  voltages. 

Based  on  the  new  theory,  we  have  calculated  a  wide  range  of 
voltage  waveforms  which  would  be  measured  on  short  and  long  lines 
under  a  broad  range  of  conditions.  We  noted  that  as  the  ground 
conductivity  decreased,  the  horizontal  electric  field  coupling 
dominated,  giving  rise  to  induced  voltages  of  negative  or  positive 
polarity  depending  on  the  location  of  the  ground  strike  point. 

Comparisons  were  made  for  individual  voltage  data  and 
theoretically  calculated  waveforms  with  reasonable  good  results. 
However,  we  find  a  discrepancy  in  peak  amplitude  between  the  data  and 
theory  which  is  attributed  to  errors  in  measurement  and/or  flaws  in 
the  theory,  as  discussed  below.  Some  general  comparisons  are  also 
made  between  our  theory  and  measurements  made  on  a  9.9  km  South 
African  power  line  (Eriksson  et  al . ,  1982)  and  measurements  made  on  a 
2.1  km  French  telecommunication  line  (Hamelin  et  al.,  1979).  In  both 
cases  we  find  good  agreement  between  theoretically  obtained  values 
and  the  reported  data. 

Some  data  on  induced  voltages  produced  by  nearby  stepped  leaders 
together  with  simultaneous  vertical  electric  field  measurement  are 
also  presented.  It  is  concluded  that  voltages  induced  on  our  test- 
line  by  nearby  stepped  leaders  are  about  an  order  of  magnitude  lower 
than  that  induced  by  the  following  return  stroke. 
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One  of  the  basic  problems  underlying  the  theory  presented  in 
this  thesis  is  the  fact  that  even  for  very  low  values  of  ground 
conductivity,  10"^  mhos/m,  the  magnitude  of  the  horizontal  electric 
field  is  not  large  enough  to  produce  the  measured  peak  voltage 
magnitudes.  Conductivities  of  the  order  of  10"-^  to  10"^  mhos/m  are 
generally  considered  the  smallest  reasonable  values  and,  as  noted  in 
Chapter  IV,  10"-^  mhos/m  was  measured  at  the  test-line  three  years 
after  the  experiment.  On  the  other  hand,  the  average  measured  peak 
voltages,  23  kV  for  first  strokes  and  11  kV  for  subsequent  strokes, 
are  an  order  of  magnitude  larger  than  those  previously  reported  if  we 
take  into  account  that  the  average  distance  of  lightning  from  the  460 
m  test-line  is  about  8  km,  so  there  are  some  doubts  about  the 
validity  of  the  voltage  calibrations.  Additional  questions  about  the 
voltage  measurements  are  raised  in  Appendix  D. 

Another  apparent  difficulty  is  associated  with  the  absence  of 
reflections  in  the  experimental  data.  If  the  experimental  test-line 
were  indeed  open-circuited  at  the  south  end,  why  are  no  reflections 
observed  in  the  measured  waveform?  The  answer  to  this  question 
probably  lies  in  the  fact  that  there  was  sufficient  damping  of  surges 
along  the  line  such  that  the  reflections  were  attenuated.  We  pointed 
out  in  Chapter  III  that  the  damping  resistance  is  really  frequency 
dependent  with  larger  damping  at  high  frequencies.  However,  in  our 
time-domain  calculations  we  determined  a  fixed  value  for  the  damping 
resistance  at  10^  'Hz  and  used  it  for  the  complete  time-domain  return 
stroke  waveform  which  has  its  frequency  spectrum  concentrated  between 
10  and  10°  Hz.  Thus  we  introduced  additional  damping  at  frequencies 
below  10^  Hz  and  adequate  damping  for  frequencies  around  10^  Hz.   In 
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spite  of  this,  calculated  waveforms  indicate  identifiable  reflections 
for  normal  values  of  conductivity,  e.g.,  10"-^  mhos/m.  This  is  again 
another  fact  which  points  to  the  validity  of  the  low  value  of 
conductivity  used  in  this  thesis.  As  an  alternative  to  the  use  of 
attenuation  to  eliminate  reflections,  if  the  south  end  of  the  test- 
line  were  somehow  terminated  in  about  500  Ohms,  the  correlation 
between  the  measured  and  calculated  waveforms  would  be  good  even  at 
10"  mhos/m,  except  that  the  calculated  peak  amplitudes  would  ba 
about  an  order  of  magnitude  lower  than  those  measured. 

Obviously,  a  direct  knowledge  of  the  horizontal  electric  field 
at  the  surface  of  the  earth  is  the  missing  link  in  our  understanding 
of  the  experimental  data.  Since  no  horizontal  electric  fields  due  to 
lightning  have  been  measured  to  date,  we  have  had  to  resort  to  using 
the  wavetilt  formulation  to  compute  a  horizontal  electric  field  from 
the  measured  vertical  electric  field.  This  theory  may  exhibit  some 
errors,  as  discussed  in  Chapter  III. 

An  experiment  should  be  performed  to  measure  the  complete 
lightning  electromagnetic  field.  The  simultaneous  measurement  of  the 
horizontal  electric  field,  vertical  electric  field,  and  the 
horizontal  magnetic  field  should  be  used  to  examine  the  assumptions 
made  in  this  thesis  and  introduce  correction  terms  wherever  needed. 
The  experimental  set-up  should  be  used  to  measure  the  complete 
lightning  field  in  the  presence  of  different  soil  conditions. 
Simultaneous  measurement  of  the  voltage  on  a  test-line  could  also  be 
recorded.  Different  terminations  could  be  used  at  different  times 
during  a  recording  season  in  order  to  see  any  differences  in  the 
measured  voltages.   The  coupling  analysis  presented  in  this  thesis 
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could  be  directly  applied  to  the  test-line  using  the  measured 
electric  fields,  thus  providing  a  check  of  that  portion  of  the 
theory. 


APPENDIX  A 

ELECTRIC  AND  MAGNETIC  FIELDS  PRODUCED  BY 
A  VERTICAL  DIPOLE  ABOVE  GROUND 


The  electric  and  magnetic  fields  at  any  arbitrary  point  in  space 
due  to  a  small  isolated  vertical  dipole  antenna  may  be  derived  in  the 
time  domain  from  Maxwell's  Equations.  For  a  homogenous,  isotropic 
medium, 

V  X  B  =   pj  +  yejl  (A.l) 

-VxE  =  -|  (A.2) 

e V    .  E   =   p  (A. 3) 

y   V   •  H  =  0   .  (A. 4) 

We  will    obtain   the   solution   to  Maxwell's   equation   in  terms   of  scalar 
and  vector  potentials.     Since   V«B  =0,  we  may  write 

B  =  V  X  A  .  (A. 5) 

Substituting  Equation  (A. 5)  into  Equation  (A.2)  we  obtain, 

-V  X  (E  *  -|)  =  0; 
which  leads  to 
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i.e.  I   =-V^  -  ^  .  (A. 6) 

Note  that  A  and  <^  are  arbitrary  functions  at  this  stage:  they  are  to 

be  determined. 

Also,  we  have  (using  Equation  (A.3)), 


v.(-v*--|)=| 


I.e.  V  (t,  +  -^  (V  .A)  =  -  I  (A. 7) 

From  Equations  (A.l),  (A. 5)  and  (A. 6)  we  may  write 

which  when  simplified  gives 

2 


V  (v-A)  _7^A  =  pO-yeV^-  ue — 

at 


at    .^2 


(A. 8) 


We  note  that  Equations  (A. 5)  through  (A. 8)  contain  the  same 
information  as  Equations  (A.l)  through  (A. 4).  What  we  have  done,  in 
essence,  is  to  transform  the  variables  in  the  equations.  To  uniquely 
define  the  vector  A  we  have  to  specify  its  curl  and  divergence.  In 
Equation  (A. 5)  we  specified  the  curl  of  A.  We  now  specify  its 
divergence  as,  '^ 


V  •  A  =  -y£^  ,  (A. 9) 
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since  it  leads  to  decoupling  of  the  equations  for  (>  and  A  .   Using 
Equation  (A. 9)  into  Equations  (A. 7)  and  (A. 8)  we  obtain 


V  <!.-  ytA-*=  --B  (A.IO) 

3t2    ^ 

2        A 

and      V  A  -  pe^-^  =  -  uJ.  (A. 11) 

at  2 


Thus,  we  need  to  solve  for  A  and  ()>  from  the  above  equations  and 
then  use  Equations  (A. 5)  and  (A. 6)  to  obtain  the  required  expressions 
for  the  fields. 

However,  we  do  not  really  need  to  solve  explicitly  for  (j)  since 
Equation  (A. 9)  relates  ij)  to  A  directly, 

1   ^ 
^=  -  -^     /(V-A)  dt  .  (A. 12) 

The  solution  to  Equation  (A. 11)  is 

J  (R'.  t  -  "^  :  ^'^) 

A(R,  t)  =^   /   ? dV  (A. 13) 

^"^   V       |R  -  R'l 

2     1 

where  c  =  — r  is  the  velocity  of  light  in  the  medium,  and  the 
integral  is  taken  over  the  source  volume.  The  solution  to  Equation 
(A. 11)  obtained  in  Equation  (A. 13)  may  be  veri/ied  by  direct 
substitution. 

Once  A  has  been  determined,  the  expressions  for  the  fields  may 
be  written  down  using  Equations  (A. 5),  (A. 6),  and  (A. 12), 

B(R,t)  =  V  X  A  ,  (A. 14) 
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F(R,t)  =  -^+  c^  /  V(V  .  A)  dt  .  (A. 15) 


Equations  (A. 13),  (A. 14),  and  (A. 15)  will  now  be  applied  to  the 
problem  of  a  small  dipole  antenna  carrying  a  current  i(t).  Consider 
the  point  P(r,4>,z)  in  Figure  A.l  at  which  we  need  to  determine  the 
fields.  Since  the  problem  has  cylindrical  symmetry,  we  shall  compute 
the  electric  and  magnetic  fields  using  a  cylindrical  coordinate 
system.  Uman  et  al.  (1975)  have  done  the  same  problem  in  spherical 
coordinates.  The  first  step  is  to  determine  the  vector 
potential  A  of  a  dipole  located  at  the  origin  of  the  coordinate 
system. 


M-  i(o,t  -  ^) 
^(R.t)  =  {  ^ ^ ^~  L  }  a"^  =  A^{R,t)a"^  (A. 16) 


Using  Equation  (A. 16)  we  obtain, 

5(r,(t),z,t)  =  V  X  A 

3A  / 

u„L  ,    i (t  -  ^) 
B(r.*,z.t)=-  [^l,     I^T^l  ]a;.  (A.17) 

We  have  the  following  relationships, 

2     2     2 

r  =  X  +  y  ; 

2     2     2 

and,   R  =  r  +  z  .  (A. 18) 


Figure  A.l.  Definition  of  the  geometrical  factors  used  in  the 
derivation  of  the  electric  and  magnetic  fields  at 
any  arbitrary  point  due  to  a  current  carrying  dipole 
at  the  origin. 
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Pir,<P,z) 
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Hence, 


and, 


^lilt-|)  )  =  -^.,-(t-|)  ; 


4int-l)   }=r(t-i)  ; 


where  primes  denote  derivative  with  respect  to  the  whole  argument, 
t  -  R/c.  Hence, 


Also, 

IF  '-r'  =  -  ;7  (A.20) 


Going  back  to  Equation  (A. 17),  we  have 


y  L     8i(t-  h  R   a  1    > 

i(r,M,t)  =  -^[^— 3^.i(t-§)^(^)  ]a- 


Using  Equations  (A. 19)  and  (A.20),  we  find 

\'    , 

y  L  3i{t-  -)  o 

5(r,,,z,t)   =^[^^—3^.-1.   i(t-^)    ]a^^  (A.21) 

cR  ^ 

We  now  determine  the  expression  for  E. 

9A,       y_L     3      i(t-^) 


3z  411      9z    ^  R  J 

R3 
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Since  there  is  no  (^variation,  the  gradient  operator  becomes, 

Thus, 

R 

uL    -   ,   8i(t-^) 

T^  I-  8z  ^  TT    3z     ^  "  9Z 

Expanding  Equation   (A. 23)  we  have, 

v(v.S)   _    r  1     ''^'^^-^^  --     '^"^^-^^ 


+  ^  [  ^  i  i         ,       ^      1  -     -^  I  ^  i(t-  ^)    }   ]  a"        (A. 23) 

r3  C  Z 


UqL      "   ^  R        3r3z        ^  ^-  ^3^         8r~ 


3r  R  z    9^'(t-F^        . 


,       3i(t-|)                         3i(t-^) 
+    [i    -^+    (--^)   ^ ^i(t-^) 


3Z 


2  n3  9Z  1,3  C 


3z  R  z    ^■'(^-F^        . 


We  have  from  Equation  (A. 19) 

3  R  r    9i(t-;^) 

Similarly, 

3  R  7      9i(t--^) 
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Using  these  relations,  we  have 


''"'-"'  3,         z'^'^-l' 


3^2  9Z     '        cR  3t 


cR         9t  C     ^  [^3    ^  at  cR         8z3t 


3i(t-^)  2  9,-(t    R)  2  a'iCt-^) 

+ -^ 5r-^  +  — ^ ^         (A. 261 


^^  ^^  CR3  ^  c2r2         ^2 

Similarly,   it  may  be  shown  that 


32i(t-^)       ^^      3i(t-^)  ^^    a2i{t-^) 


cR3  ''^  c2r2  at2 

Substituting     Equations     (A. 19),      (A. 25),      (A. 26),     and     (A. 27)     into 
Equation   (A. 24),  we  find, 

r5  ^       c'       j^3  ^      cR       at  ^   i\ 

3i(t-^)        ^2    a2i(t-f) 


Ml 

(  f-   '  +  '' 

-^' 

z     ^'-'t-F' 

cR        at 

^  c2r2  3t 


-i(t-^)  +  3l!    i(t    R) 
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8i(t-^) 


Z     r         _Z  C 

3    ^    '   CR  3t  J     J   ''z 


]  a 

R 


r  3rz  R         3rz    ^^^^-  F^         r7      ^i^t--^)       , 

R5  ^         cR'*  ^^  c2r3        3t2  1^ 

R5        R3  C  ^p,        ^j^2  3t 


,2      32i(t-^)         ^ 

'      ]a",  (A.28) 


c2r3  8t2  ^      Z 


Hence  we  can  write. 


t 

=  c2     /  v(  V  .  A)  dt 


cj 

0     ■     ■    R5      o'  "^  CR"^  ^'  c2r3  '^^  ^    "^ 


L     r   ;  3rz      r\-(.    Rn   .,.  3rz  ,,.     R,    ,     rz     ^i^*--) 
t 


R5         R3  0  ^  cR'^       cR2  ^ 

,2      9i(t-§)         ^ 

+  -^ ^  }a"     J  (A. 291 

c2r3  9t2  '      Z    ' 


Also,   from  Equation   (A. 16),  we  have 


0 


In  Equations  (A. 29)  and  (A. 30)  we  have  used  the  relation, 
.2_1 


[A. 301 


(A. 31) 
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Finally,  we  are  in  a  position  to  write  down  the  complete  expression 
for  E  from  Equation  (A. 15).  For  the  sake  of  convenience  we  write 

^  =  ^r^""  ^z^z  (A. 32) 

where  the  components  in  Equation  (A. 32)  may  be  obtained  by 
substituting  Equations  (A. 29)  and  (A. 30)  into  Equation  (A. 15). 


E^{r,*.z,t)  =4^[^   A-(T-^)  ...^nt-l) 


rz   9i(t-f) 
+  ^ 5F-^  ]  (A. 33) 


:2r3 


at 


and 


E,(r.z,,,t)=^[  (21^]  A(.-f)d..  i^^^]Ht-l) 
^  ^  %  R5    0      ^        cR"^        ^ 


,2      3i(t-f) 


;2r3    9t 


(A. 34) 


Equations  (A. 21),  (A. 33)  and  (A. 34)  represent  the  field  solution 
due  to  a  small  dipole  at  the  origin.  For  our  purpose,  the  dipole 
will  lie  along  the  z-axis  at  an  arbitrary  source  coordinate  z' .  In 
that  case,   the  solution  is  modified  and  may  be  written  as, 

y  L                ai(z',t  -  ■^)         .  n 

S^'-'*'^'^^  -^^-~ at -'f,  -n^'.t-l)]         (A.35) 

CK  K 
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E^(r,*,z,t)   = 


L      r  3r(z-z') 


0  R^'        0 


/i(z',x-  J)dx-.MAlz:ii(,.,t-i) 


cR' 


r(z-z')      9^*(z'.^-c^ 


c2r3 


at 


;a.36) 


flr^7f\-       L      r  2(z-z')^-r2      }■   ...  R.. 

t^(r,(j),z,t)   -  4^j^   L J     i(z   ,T--)dT 


?f7_7M2_r2  D  ^2         ai(z',t--) 


CR'*  ^'         c2r3  3t 


APPENDIX  B 
TEST  LINE  VOLTAGES  INDUCED  BY  NEARBY  STEPPED  LEADERS 

In  this  Appendix,  data  on  voltages  induced  by  nearby  stepped 
leaders  are  presented.  Simultaneously  measured  vertical  electric 
field  pulses  due  to  the  stepped  leader  will  also  be  shown. 

Figure  B.l  shows  the  vertical  electric  field  change  produced  by 
a  single  stroke  flash  which  occurred  at  185337  UT  during  run 
#79196TR27  on  a  time  scale  of  10  ms/div:  both  direct  and  PI  channels 
are  given.  In  Figure  B.2  the  location  of  the  ground  strike  point  is 
indicated. 

Instrumentation  tape  records  of  the  return  stroke  field  change 
and  the  corresponding  voltage  on  the  unenergized  test  line  are  shown 
in  Figure  B.3  on  a  time  scale  of  200  \s  full-scale.  The  peak  voltage 
is  about  150  kV  and  is  of  positive  polarity.  Sampling  the  continuous 
voltage  and  field  records  at  a  faster  sampling  rate  and  at  increased 
sensitivity,  we  can  look  at  the  waveforms  in  detail,  for 
about  50  MS  preceeding  the  return  stroke.  In  Figure  B.4  this  detail 
is  shown  on  two  expanded  time  scales:  the  top  photograph 
at  8  us/div  and  the  bottom  at  4  lis/div.  The  stepped  leader  pulses 
can  be  easily  identified  on  the  vertical  electric  field  records. 
Corresponding  line  voltages  can  also  be  identified  over  the  noise 
level  of  the  tape  channel. 
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The  peak  vertical  electric  field  due  to  the  stepped  leader 
pulses  is  an  order  of  magnitude  lower  than  the  return  stroke 
radiation  peak.  This  is  in  conformity  with  stepped  leader  electric 
field  measurements  reported  by  Krider  et  al.  (1977).  It  is  an 
interesting  fact  that  the  peak  induced  voltage  on  the  UTL  is  also  an 
order  of  magnitude  lower  i.e.,  about  13  kV,  than  that  due  to  the 
following  return  stroke. 

This  method  of  trying  to  extract  the  stepped  leader  induced 
voltages  was  tried  for  several  other  strokes.  However,  in  all  cases, 
the  induced  line  voltage  was  found  to  be  buried  below  the  channel 
noise  level.  Thus,  it  may  be  concluded  that  the  peak  voltage  induced 
on  the  UTL  by  stepped  leaders  is  an  order  of  magnitude  lower  than 
that  induced  by  the  following  return  stroke.  These  data  represent 
the  first  simultaneous  recording  of  the  voltages  on  an  overhead  line 
and  the  vertical  electric  field  due  to  stepped  leaders,  and  are  in 
contradiction  with  calculated  values  of  the  order  of  1  MV  derived  by 
Owa  (1964),  based  on  the  pre-discharge  theory  put  forth  by  Griscom 
(1958).  This  is  not  unexpected  since  the  calculation  by  Owa  (1964) 
assumed  a  stepped  leader  current  of  100  kA,  whereas  modern  measure- 
ments reported  by  Krider  and  Radda  (1975)  and  by  Krider  et  al .  (1977) 
indicate  a  peak  stepped  leader  current  of  the  order  of  1  kA. 


APPENDIX  C 
COMPUTER  PROGRAMS 

In  this  Appendix  we  give  FORTRAN  listings  of  various  computer 
codes  which  have  been  written  as  part  of  this  dissertation.  A  brief 
note  of  their  usage  is  also  included. 

C.l  Return  Stroke  Program 

The  user  specifies  the  various  parameters  associated  with  the 
three  component  currents  derived  from  two  station  field  measurements 
as  given  by  Lin  et  a1.  (1980).  The  field  point  is  specified,  and  the 
printer  and  plotter  format  is  also  specified.  As  output  the  program 
gives  the  electric  and  magnetic  field  components  at  the  field  point 
and  the  return  stroke  current  at  ground  for  each  time  step.  As  an 
option,  the  electrostatic,  induction,  and  radiation  components  of  the 
fields  may  also  be  printed.  A  FORTRAN  listing  follows. 
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C.2  Parameters  of  Lossless  Transmission  Lines 

The  following  is  a  FORTRAN  listing  of  a  computer  program  which 
can  be  used  to  determine  the  inductance  and  capacitance  matrices  for 
any  overhead  transmission  line.  The  program  allows  lines  with 
multiple  phases,  multiple  conductors  per  phase  (bundling)  and 
multiple  skywires. 
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C.3  Wavetilt  Program 


The  following  FORTRAN  listing  is  for  a  program  written  to 
compute  the  horizontal  electric  field  from  the  vertical  electric 
field  using  the  wavetilt  formulation  described  in  Chapter  III.  The 
user  specifies  a  piecewise  linear  version  of  the  vertical  electric 
field,  and  the  soil  parameters  a  and  e  .  The  computation  time  step 
and  the  total  time  of  interest  are  also  specified.  The  vertical  and 
horizontal  fields  are  plotted  on  the  Gould  plotter  as  output.  The 
vertical  and  horizontal  electric  fields  as  a  function  of  time  are 
written  onto  a  disk  file  to  be  used  in  the  coupling  program  listed  in 
the  next  section. 
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C.4  Coupling  Program 

The  following  is  a  FORTRAN  listing  of  the  program  used  to 
determine  the  voltages  and  currents  induced  on  overhead  lines  by 
nearby  lightning.  The  inductance  and  capacitance  parameters  are 
specified  as  obtained  from  the  program  listed  in  Section  C.2.  The 
program  allows  parallel  R-C  termination  at  either  end.  The 
coordinates  of  the  line,  the  location  of  the  ground  strike  point,  the 
time  step,  the  number  of  segments  into  which  the  line  is  split,  and 
the  total  time  for  which  calculations  are  desired  are  all  specified 
as  input.  The  incident  fields  are  read  off  the  disk  file  generated 
by  the  program  listed  in  Section  C.3.  The  output  presents  the 
voltage  plotted  on  the  Gould  plotter  as  a  function  of  time  for  any 
point  on  the  line. 
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APPENDIX  D 
DISTRIBUTION  OF  PEAK  VOLTAGES 

In  this  Appendix  we  present  a  distribution  of  the  peak  voltages 
induced  on  the  test-line  by  first  return  strokes  during  a  one  hour 
and  thirty  minute  period  of  run  #79220TR57.  The  storm  was  primarily 
to  the  S-SW  of  the  Mobile  Lightning  Laboratory  (MLL)  and  in  the  end 
moved  east.   Due  to  the  location  of  the  storm,  most  of  the  induced 
line  voltages  were  of  positive  polarity.   Unfortunately,  the  TV 
network  was  not  operating  during  this  period  and  hence  no  flashes 
were  located.  However,  from  the  records  kept  in  the  data  book  it  was 
clear  that  during  the  last  half  hour  the  storm  moved  very  close  to 
the  MLL.  Hence,  relatively  large  voltages  were  induced  on  the  test- 
line  at  that  time  due  to  flashes  which  were  much  closer  to  the  test- 
line  than  any  presented  in  the  body  of  this  thesis.  In  Figure  D.l  we 
present  the  distribution  of  the  measured  voltages.   We  note  that 
there  are  induced  voltages  in  excess  of  250  kV.  Voltages  this  large 
do  not  appear  to  be  reasonable  because  the  line  was  designed  to 
spark-over  at  250  kV.   Thus,  either  the  line  spark-over  gap  was 
improperly  set  or  there  was  an  error  in  the  voltage  calibration  so 
that  the  actual  voltages  were  smaller  than  the  measured  values 
indicated. 
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